Field trial evidence of non-transgenic and transgenic Bt. rice genotypes in north of Iran

Background Field-testing genetically modified crops provides scientists with an opportunity to collect information on environmental interactions and agronomic performance, which is critical to a full environmental safety assessment as required by regulatory authorities. As a result, the goal of this research was field trial of transgenic and conventional rice genotypes. The experiment was carried out in a randomized complete blocks design (RCBD) with four replications and seven genotypes in three isolated regions under the biosafety standard protocol in north of Iran in 2016. In this study, four transgenic lines with an active cry1Ab gene in the vegetative stage and three conventional genotypes (control) as treatment were assessed. Results The findings demonstrated that in all three regions, transgenic lines derived from Khazar cultivar, were similar to their parent(s) in terms of growth phenology, agronomical traits, grain amylose content, gel consistency, and gelatinization temperature. In all the three regions, the highest number of panicle per m2, number of filled spikelet per panicle, and filled spikelet percentage per panicle were obtained for transgenic lines. Khazar cultivar compared to the transgenic lines showed lower paddy yield. In all the three regions, transgenic lines had lower yield loss than their non-transgenic parent. The lowest number of white heads belonged to transgenic Tarom Molaii and transgenic lines. The most positive direct effect on paddy yield was related to the number of filled spikelet per panicle. Thus, this can be a good trait to achieve higher yield derived from reducing the negative indirect effect of dead heart and white heads. Conclusion It can be concluded that by producing transgenic rice, which is resistant to stem borer with an active promoter in the reproductive stage, farmers can reduce significant part of yield loss resulting from dead heart and white heads directly correlating with the number of filled spikelet per panicle and paddy yield.


Background
Rice (Oryza sativa L.) is the staple food of more than half of the world's population and has an obvious effect on feeding, income, and job creation of people in the world, especially in Iran [27]. According to the official statistics released by FAO, the area cultivated with rice in the world during the past years was from 145 million to over 160 million hectares [8]. The last global statistics showed that paddy yield and white rice productions were 742 and 492.2 million tons respectively in 2014, respectively [8]. The same amount was predicted for 2016. Iran has 550,000 ha of paddy field and two million tons of white rice production and has a 0.4% share in rice production and cultivation area in the world, most of which (about 75%) is located in the northern strip, i.e., the provinces of Guilan, Mazandaran, and Golestan, and the remaining (25%) of the paddy field is located in other 13 provinces with different climates [23]. Rice production in the world has been facing serious challenges since 2000 making the international community concerned about food supply for people in the world. As shown by investigations, the most important method to fight rice stem borer (Chilo suppressalis Walker) in Iran and also most regions in the world has been the chemical inputs [23]. Therefore, making transgenic rice cultivars by the gene cry1Ab of the bacterium Bt. is the best method to give plants resistance against the invasion of pests and to reduce the environmental damage by chemical pesticides [9].
Plant biotechnology has traditionally encompassed the application of cell and tissue culture for crop improvement. Since the mid-1980s, the development and application of transgenic plants has been the dominant research activity associated with plant biotechnology. More recently, research has expanded to include the application of genomics technologies [2]. Field-testing genetically modified crops provides scientists with an opportunity to collect information on environmental interactions and agronomic performance, critical to a full environmental safety assessment as required by regulatory authorities.
To estimate the yield benefit and cost of insectresistant GM rice, Xia et al. [39] compared field productivity of three insect-resistant rice lines with their non-transgenic parental cultivars. Their findings revealed that great benefits obtained by the yield-related traits were detected in the transgenic lines when high insect pressure was recorded, but a cryptic yield loss was observed when the level of insect pressure was extremely low. In fact, yield loss produced in the transgenic lines was under low insect pressure. A strategic field extension should be required when insect-resistant lines are commercialized to circumvent the unnecessary yield losses; this probably applies to other insect-resistant transgenic crops [39].
Among the agronomic traits of rice, yield is the most important one. It possesses three main components: number of tillers (panicles) per plant, number of grains per panicle, and 1000-grain weight. Yield is a complex quantitative trait; in addition to the three main components, there are also some other traits relevant to rice yield, such as plant height, plant architecture, biomass, and seed-setting rate [18,24,35]. Therefore, in the process of improving rice yield trait, we should take multiple agronomic traits into consideration [6]. Klumper and Qaim [21] considered a meta-analysis of the impacts of genetically modified crops. They carried out a meta-analysis of the agronomic and economic impacts of transgenic crops to consolidate the evidence. On average, GM technology adoption has reduced chemical pesticide use by 37%, increased crop yields by 22%, and increased farmer profits by 68%. Yield gains and pesticide reductions are larger for insect-resistant crops than for herbicide-tolerant crops. Yield and profit gains are higher in developing countries than in developed countries [21]. Perhaps great economic and impressive approach to protecting rice against insect attack is to cultivate insect-resistant rice cultivars. Overall, the objectives of this experiment were (1) to examine the quantitative and qualitative parameters of nontransgenic and transgenic rice cultivars, (2) to estimate the pest loss and stem borer resistance parameters including dead heart and whit heads, and (3) to identify genetic correlation and regression model by path analysis.

Description of the regions
Experiments were done in isolated paddy fields in three sites, under the control of Iran Rice Research Institute between Alborz Mountains range and the Caspian Sea in Mazandaran and Guilan provinces in north of Iran (Table 1) in 2016. The geographical coordinates and soil properties of the three sites are shown in Table 1. Solar radiation was estimated using sunshine hours and extraterrestrial radiation [31][32][33]. For calculating solar radiation, the Srad_calc program was used. This program uses the sunshine hour data to calculate solar radiation. For counting day length, the PP_calc program was applied. Srad_calc and PP_calc programs can also be downloaded from "https://sites.google.com/site/Crop-Modeling."

Description of the experiment
The experiment was carried out in a randomized complete blocks design (RCBD) with four replications. In this experiment, there were three transgenic rice lines including KHT 2 , KHT 3, and KHT 4 (driven from the backcross of Khazar cultivars with the transgenic line of Tarom Molaii) containing the gene cry1Ab from the bacterium Bt. (Bacillus thuringiensis) with the transgenic line of Tarom Molaii (containing the gene for resistance against striped stem borer as the non-restorer parent) along with the non-transgenic cultivars. Three nontransgenic rice cultivars were Tarom Hashemi, Nemat, and Khazar. In fact, in this experiment breeding and transgenic of cultivars were considered to use genetic potential of Tarom Molaii cultivar. Therefore, in backcross, 12.5% of Tarom Molaii background was maintained and in F 2 . In addition to transgenic selection, selection was made based on superior rice characteristics and morphological traits. The diagram of breeding scheme of backcross lines is shown in Fig. 1. The profile of the genotypes is shown in Table 1.
The isolation conditions were met according to biosafety principles (standard of biosafety protocol) consisting of cultivation with distance from other farmers' Table 1 Description of the geographical coordinate, soil properties (0-30 cm) and climatic parameters of three rice production sites. As well as description of name, origin, and other characteristics of rice genotypes  paddy fields. This allowed for putting enough margin space around the field, harvesting the margin separately, and eliminating the plants' residue after the harvesting. The experiment's time and space distance were also considered adhering to biosafety principles. Considering the climates of Mazandaran and Guilan provinces, seedlings were transplanted in 3-4 leaf stages. Considering the type of the cultivar, the transplanting operation was done in all the three regions with similar situations. The size of each plot was decided to be 4 × 7 m 2 , and planting density was 16 plants per square meter by 25 × 25 cm 2 planting arrangement. Nitrogen, phosphorus, and potassium fertilizers were used in each region according to the suggestions of Iran Rice Research Institute and by considering the result of soil analysis ( Table 1). All phosphorus amount and one-third of nitrogen fertilizer were used as basal in paddy field preparation stage. Two-thirds of nitrogen fertilizer was used as top-dressing in panicle initiation and full heading stages. Sixty percent of potassium fertilizer was used as basal and the remaining amount was used as topdressing in tillering and panicle initiation stages (splitting to 20% in each stage). Depth of irrigation was set at 5 cm according to agricultural principles. In order to control weeds' growth and the mixing of water and fertilizer in paddy fields, nylon plastic cover was put at the borders to the depth of 30 cm.
Crop protection practices, irrigation, weeding, and fertilization, were done in the isolated experiment paddy field in each region. But, no pesticide was used so that the lines' resistance could be accurately observed. Nits of stem borer were taken from infected fields and distributed in the experimental fields ultra-optimally. Other crop management practices were done according to the standards of biosafety protocol and Standard Evaluation System (SES) of International Rice Research Institute (IRRI).

Agronomical traits
During the growing period, after the removal of marginal effect in each site, all agronomical traits (phenological traits, morphological traits, and yield components) were randomly measured according to Standard Evaluation System (SES) of International Rice Research Institute (IRRI). Paddy yield, straw yield, and biological yield were measured by harvesting hills from 4 m 2 in the middle part of each plot based on 12% of moisture. From the relationship between paddy yield and biological yield, the harvest index was calculated and expressed in percentage.
The grain filling period was determined based on the time between pollination and physiological maturity (50% of the panicles turned yellowish green). The average of grain filling rate in the unit of area was obtained by dividing the paddy yield based on the length of grain filling period according to g m −2 d −1 . According to the difference in growing pattern and cultivars' maturity, the thermal index of GDD was calculated for grain filling period from pollination to physiological maturity stages according to the formula below: Where T max is maximum daily temperature, T min is minimum daily temperature, and T b is based temperature. The temperatures below 10 and above 30°C were considered 10 and 30, respectively. In order to calculate leaf area index the during pollination stage, the maximum length and width for all leaves in ten hills per plant were measured in centimeters, and it was estimated by the formula LA = 0.75 × LW, where LA is leaf area, and L and W are the leaf's length and width, respectively.

Determination of grain quality traits
Measuring baking qualitative traits, such as grain amylose content, gelatinization temperature, gel consistency, grain length, grain width, grain elongation, baked grain length, baked grain width, baked grain elongation, conversion efficiency, and healthy and broken grain percentage was done by Standard Evaluation System (SES) of International Rice Research Institute (IRRI).

Estimation of yield loss
For estimating yield loss, Islam and Karim's [15] method was used. In this method, in the white head stage, 16 plants were chosen per experimental plot. In every plant, the measured number of healthy panicle (N), the number of fairly healthy panicle (N´), and the number of blank panicle (N´´) were measured. The number of healthy spikelet per panicle (P), the number of fairly healthy spikelet per panicle (P´), and total paddy yield (Y) were measured randomly by 16 plants per plot.

Evaluation of pest loss
Sorting of infected percentage of plants for determining sensitive and resistance lines was evaluated by Heinrichs [14] method. In order to determine the dead heart, 30, 50, and 65 days after transplanting, and in order to determine the white heads, sampling was done a week before the harvest. In every sampling from each plot, 10 hills per plant were randomly selected, and the number of plants with dead hearts (DH) related to the damage from the insect's first generation, and the number of white heads (WH) related to the damage from the

Statistical analysis
After normalization, data analyzed by the SAS statistical software and averages' comparison were calculated by LSD tests in a 5% probability level. To determine those traits that had the most effect on yield, after applying regression to eliminate alignment, path analysis was performed by genotype correlation with the Path software.

Combined mean square analysis of investigated traits
Findings of the compound analysis variance in Table 2 showed that all investigated traits including phenological traits, morphological traits, yield components, quantitative yield and harvest index, yield loss, dead heart, and white heads were significant in 1% probability level on genotype treatment. DT, DF, DFM, SL, FS, SY, BY, HI, and YL were significant in 5% probability level under the effect of regions; also DFM, PL, TH, FTH, ITH, PM, FSP, BS, BSP, TGW, PY, DH, and WH were significant in 1% probability level on region. At interaction of region × genotype, DFM, FLL, FTP, ITH, ITP, and HI were significant statistically in 1% probability level. FLA, LAI, FTH, PM, YL, and WH were significant in 5% probability level under the interaction of region and genotypes (Table 2).

Mean comparison analysis of investigated traits
According to the findings, the vegetative and reproductive periods in Rasht region were higher than two other sites ( Table 2). The lowest vegetative and reproductive growth periods of cultivars were observed in Sari region. Nemat cultivar and transgenic lines derived from Khazar cultivar were long-term in terms of vegetative and reproductive growth period, and Khazar cultivar was in the next rank. Two local cultivars, Tarom Hashemi and Tarom Molaii, stood in third ranks. Attending to interaction of location and genotypes, we can find out in the three sites that the highest DFM was observed for nontransgenic Nemat cultivar and transgenic lines including KHT 2 , KHT 3, and KHT 4 . In the three sites, the least growth period was obtained for local cultivars, including Tarom Hashemi and Tarom Molaii ( Table 3). The highest and lowest FP was observed for all cultivars in Rasht and Sari. GDD in pollination stage varies in different regions, indicating that temperature variations vary the overall rice growth period. Mean comparison of genotypes showed that the maximum grain FP was related to the two KHT 2 and KHT 3 lines, Nemat, Khazar and Tarom Hashemi cultivars. Nemat and Khazar cultivars along with KHT 2 and KHT 3 lines had the highest GDD in pollination stage. The lowest grain FP and GDD in pollination stage was obtained for Tarom Molaii cultivar, but the least grain FR was related to Khazar cultivar (Table 2). Genetic factors (cultivar) somewhat determine the grain FR, and environmental factors (temperature) determine the extent of grain FP. It seems that grain FR of most cultivars was higher due to PM. When number of panicle per of unit area is higher, it causes faster source depletion (leaf and stem).
Mean comparison of treatments showed that among genotypes, Tarom Hashemi and Tarom Molaii had the highest PL and SL. Among semi-dwarf cultivars, the lowest PL and SL belonged to Nemat, Khazar and transgenic lines derived from Khazar. The highest LAI belonged to Khazar cultivar, the main cause of which was the higher FLA of this cultivar ( Table 2). Mean comparison by slice interaction showed that in all the three regions, the highest FLL belonged to Khazar, Tarom Hashemi, and Tarom Molaii cultivars. In all the three regions, transgenic lines derived from Khazar cultivar were similar to their parents in terms of FLL and FLA. In all the three areas, the highest LAI belonged to Khazar cultivar. Nemat, Tarom Hashemi, and Tarom Molaii cultivars with a lower difference stood in later ranks. In all the three regions, the lowest LAI belonged to Khazar transgenic lines (Table 3). Mean comparison shows differences between the three regions in terms of these traits due to the climatic conditions and soil properties of these regions. Genotypes mean comparison showed that the highest PM belonged to KHT 3 and KHT 4 lines. The low PM in Tarom Molaii cultivar has its genetic potential, which has lower degree of tiller production than modified cultivars ( Table 2).
According to mean comparison, using slice interaction transgenic lines had higher FTH, fertile tiller percentage, lower ITH, and ITP compared to non-transgenic cultivars (Table 3). In all the three regions, mean comparison by slice interaction method showed that the highest PM was obtained for transgenic lines derived from Khazar cultivar. Also, in Rasht region, Nemat cultivar along with Khazar transgenic lines showed the highest PM. Local cultivars, Tarom Hashemi and Tarom Molaii, showed the lowest PM (Table 3).
Genotype mean comparison showed that the highest paddy yield belonged to Nemat cultivar and transgenic lines derived from Khazar cultivar. Khazar cultivar compared to the transgenic lines, KHT 2 , KHT 3 , and KHT 4 , showed 25.49, 25.88, and 19.26% lower paddy yield, respectively ( Table 2). The lowest paddy yield was attributed to two local cultivars, Tarom Hashemi (4520 kg ha −1 ) and Tarom Molaii (3579 kg ha −1 ). In terms of straw yield and biological yield, cultivars had a high variation (Table 2). Mean comparison by slice interaction method showed that in all the three regions, transgenic lines derived from Khazar had higher harvest index than their non-transgenic parent. Also, in all the three regions, the lowest harvest index was observed for Tarom Molaii cultivar (Table 3).

Qualitative traits
According to findings, there were no significant differences between transgenic lines and their parents in terms of technological traits and grain baking related traits. In terms of grain length before and after baking, the transgenic lines KHT 2 , KHT 3 , and KHT 4 were similar to Khazar cultivar. In fact, the coherence in the shape and paddy yield is the first factor in the quality corroboration. Farmers consider these factors for the cultivation of new cultivars ( Table 4). The grain appearance is important to retain its quality and market-friendliness. Grain length before baking is one of the qualitative factors of rice grain. This trait was not influenced by transgenic lines. Transgenic lines and their conventional parent have no significance in terms of grain length after baking and grain elongation (Table 4). In terms of grain conversion efficiency and healthy grain percentage, transgenic lines have a better condition compared to non-transgenic parents. Rice grains increase after baking, and if this increase in length is without increase in thickness, it can be a positive factor for cultivars. There were no significant differences between transgenic lines and their parents in terms of grain amylose content, gel consistency, and gelatinization temperature (Table 4).   Yield loss (YL), dead heart (DH), and white heads (WH) Based on the results, the highest YL was observed for Khazar cultivar and Nemat cultivar. Yield loss of transgenic lines derived from Khazar cultivar was significantly lower than their non-transgenic parent. Also, the lowest YL was obtained for Tarom Molalii (Table 2). Mean comparison by slice interaction method showed that in all the three regions, transgenic lines derived from Khazar had lower YL than their non-transgenic parent. In all the three regions, the lowest YL was observed for Tarom Molaii cultivar (Table 3). Differences in dead heart and white heads in the three regions were related to pest variation in each region. Genotype mean comparison showed that the highest and lowest dead heart and white head rates belonged to Tarom Hashemi and Tarom Molaii. The dead heart rate of Khazar transgenic lines was statistically lower than their non-transgenic parent. The lowest WH belonged to KHT 4 ( Table 2). Mean comparison by slice interaction method showed that in all the three regions, the white heads belonged to Khazar cultivar. Tarom Hashemi and Nemat ranked second with fewer variations. The lowest number of white heads belonged to transgenic Tarom Molaii cultivar and transgenic lines derived from Khazar cultivar, which were diverse in terms of numerical diversity (Table 3).

Stem borer resistance traits
The results of rice genotypes ranking based on dead heart percentage and white heads percentage are presented in Table 4. In terms of dead heart, Tarom Molaii cultivar suffered no damage and was recognized as a resistant cultivar. KHT 4 line got placed in group I (relatively resistant), whose percentage of infection was below 10%. KHT 2 and KHT 3 lines and Nemat cultivar got placed in group III (medium resistance), where dead heart was below 20%. Tarom Hashemi and Khazar cultivars stood in group V (low resistance), where dead heart was equal to 24% and 17% (Table 4).

Correlation analysis
According to findings, there was a negative correlation between FS and white heads in the three sites (Fig. 2a). We observed negative correlations between paddy yield with dead heart and white heads in the three sites (Fig. 2b, c). Negative correlation between these traits showed that with the increase of dead heart and white heads, FS and paddy yield had decreased. This can be a hindrance to aggregating the high value of these variables in a genotype, especially if this negative relationship is genetic forms and due to genetic connectivity.

Path coefficient analysis
Path analysis was used to determine the direct and indirect effects of the variables introduced in the stepwise regression model in the three regions (Table 5). Therefore, genetic correlation was used for the estimation of direct and indirect traits on paddy yield. The most positive direct effect on PY was about FS per panicle. The total direct and indirect effects of this trait on PY (genetic correlation with paddy yield) were equal to 0.93. The indirect effect of FS per panicle by DF, PM, and HI was 1.179, 1.126, and 1.139, respectively (Table 5).
According to the direct effect of FS on PY, its indirect effect by DF and PM can be a proper trait to achieve higher yield which is obtained by reducing the negative indirect effect of DH and WH. The maximum negative indirect effect regarding FS per panicle was achieved by PH and WH (Table 5). However, PM has a negative direct effect on PY; but it has a positive indirect effect on FS, which can be considered as a useful trait for selecting higher yield and used in indirect selection discussion ( Table 5). The most negative indirect effect on paddy yield is about PM and WH. The sum of direct and indirect effects of these traits was 0.91 and − 0.78, respectively. The indirect effect of PM by all investigated traits was obtained to be more than one. Indirect effect of WH was positive by PM, FS, and HI; also, it has negative effects by DF, PH, ITH, and DH (Table 5).

Discussion
In all the three regions, transgenic lines derived from Khazar cultivar were similar to their parent in terms of phenological growing. According to other researchers' findings, phenology is one of the major topics in the field of ecology, and its aim is to study of changes in the vital stages of plants. Determining the vital stages of plants, including the evaluation of phenological traits, to improve yield and help to make decisions to maximize available plant sources, is very important [41]. The phenology and growth period of rice is one of the main factors determining the agronomical and ecological suitability of cultivars in the cultivated area [41]. Therefore, proper prediction of phenological stages of crops is important for optimizing management activities in the field and better adaptation of crop schedule [16,19]. In fact, phenological developmental stages in plants are affected by temperature, photoperiod, and vernalization [11,34]. Adaptation of the growth period of cultivars to achieve acceptable yields, increase productivity, reduce losses, and maintain their cultivation is always one of the important and influential topics. Therefore, the development of medium and late cultivars cultivation in paddy fields in the northwestern part of Iran is highly dependent on the crop rotational calendar, especially seeding date in nursery, seedling age, transplanting date, and growth period of rice cultivars. In fact, flag leaf, with its location close to the panicles, has a major contribution to the transfer of photosynthetic material and grain filling. Hence, in order to achieve maximum yield, the amount of leaf in the shoots of the plant is essential. Davatgar et al. [4] revealed that most of the reduction in leaf area appears to be the consequence of revealed cell expansion, the closing of stomata, and inhibition of photosynthesis.
Genotype mean comparison showed that transgenic lines derived from Khazar cultivar compared to their non-transgenic parent had a higher number and percentage of FS with lower number and percentage of ITH. The most TGW was observed in Nemat cultivar. The TGW weight belonged to Khazar cultivar and its transgenic lines and local cultivars, Tarom Hashemi and Tarom Molaii, all of which were on a same statistical level. The studied trait variation in different regions was due to the effects of environmental conditions (soil and climate) and other factors (planting and harvesting dates). Yield components are affected by management practices, genotypes, and environment and are often used to justify the cause of yield reducing or increasing. Also, an ecosystem affects the ability of a plant to exhibit genetic potential. In addition, incorrect crop management, water usage, nutrition, temperature, light, and other inappropriate environmental factors can reduce one or more components of the yield components. Also, the number of spikelet per panicle depends on genetic ability of cultivars, climatic conditions at panicle formation, pollination, size, and activity of photosynthetic system at panicle formation period and photosynthetic material transferring capacity to panicle and competition between individual plants. One thousand-grain weight is a genetic trait and is less affected by environmental conditions. This trait depends on the size and duration of the activation of carbon dioxide fixation in plant shoots, photosynthetic material capacities loading to grain, duration of grain formation, weather conditions, mineral amount in grain filling period, and pests and diseases' occurrence.
In order to overcome the issue of low productive tillers, some researchers have proposed new plant type rice with low-tillering and large panicle size [20]. Among rice tillers, the panicle development pattern is hierarchical and grain yields decrease for each successive tiller [28]. In such circumstances, the yield contribution of the inferior tillers could not be ignored because they possess high productivity potential in theory, as the totipotency of rice coleoptile tissues [26]. If the limiting factors between the inferior and superior tillers could be reduced, or even removed, the yields of the inferior tillers might be dramatically improved. In addition, hightillering rice cultivars possess a good capacity for functional compensation, and late emerging tillers may suitably compensate for yield losses when nascent tillers are subjected to environmental stresses [38]. Studying the influence of tiller heterogeneity on yield components of rice grown under a different nitrogen treatment, Wang et al. [37] indicated that the quantitative proportions and yield components were decreased in superior tillers and increased in the inferior tillers. The present study suggested that the enhancement of grain filling and weight of the inferior tillers would be an impressive approach to further improve rice yield per acre [37]. Panicle is one of the important traits for rice plant type because it can directly affect grain yield as a factor of sink size [42]. To increase the value of each panicle, an effective strategy to improve grain's function should be adopted [6]. Previous researchers have indicated that grain weight was the most stable yield component, which was a highly heritable trait less affected by environmental factors [13].
The difference between paddy yield and biomass was due to environmental conditions (climate and soil) and pests and diseases variation in each region. In fact, the formation of vegetation that has a high performance depends on controlling components that interfere with plant growth and development. Achieving the highest yield in transgenic cultivars was the result of providing optimal growth conditions. Therefore, yield loss due to damage of stem borer in transgenic cultivars was much lower than non-transgenic cultivars. Paddy yield of Nemat and Khazar cultivars and transgenic lines derived from Khazar was higher than Tarom Hashemi and Tarom Molaii cultivars. Because Nemat and Khazar cultivars are dwarf and there is the small distance between source and sink the majority of photosynthetic materials are transferred to grain. The semi-dwarf characteristic of these cultivars reduces the competition between vegetative and reproductive organs for photosynthetic materials. Their paddy yield was higher in comparison with two tall cultivars Tarom Hashemi and Tarom Molaii. This is because in tall cultivars, distance between source and sink is high, and most energy and nutrition that plays a decisive role in the transfer of photosynthetic materials to grain is used for vegetative growth. Remobilization of photosynthetic materials in rice is such that in the pollination and heading stages, produced photosynthesis material is more than the need for these two processes. The excess of photosynthetic material is transferred to stem and stored as starch, and in grain filling stage, starch materials are converted to sugar and transferred to these grains. Therefore, because of the small distance between source and sink, semi-dwarf cultivars are more likely to transfer nutrition supply to grain and have more success compared to tall cultivars. By field trials of two fertilizer levels and 18 modern cultivars on the paddy yield of rice, Li et al. [22] revealed that variation among cultivars had a great effect on paddy yield. Close correlations were observed between paddy yield and effective panicles and dry matter production. Harvest index, the ratio of grain weight to total shoot weight, is an important trait associated with the dramatic increases in crop yields. Most progress in improving harvest index occurred following the introduction of semi-dwarf traits into rice in the 1960s [30]. However, the scope for continued harvest index increases in modern rice cultivars was limited by the need for maintaining sufficient leaf area and stem biomass for interception of solar radiation, physical support, and storage of assimilates and nitrogen used in grain filling [3]. Jianchang et al. [17] revealed that post-anthesis dry matter production of high-yield rice differs among different yield categories.
Rice baking traits are mainly determined by the properties of starch, which comprises up to 90% of endosperm of white rice and consists of two amylose and amylopectin components. The high amylose rice is completely separated after drying and hardened after cooling. Most rice with lower amylose content has spikelet with a minimum amount of blank grain and number of broken grains [36]. Measurement of gel consistency complements the amylose test, which indicates the amount of baked rice gel movement. Gel consistency shows the exact hardness and adhesion of cooked rice. Determination of endosperm starch gelatinization temperature is an important test for determining the quality of baking in rice. The temperature of the gelatinized starch is an important quality factor that is related to the baking time and the nature of baked rice; in fact, the higher the gelatinization temperature, the higher the ranking is in the qualitative categorization. Studying effects of genotype and environment on bread making quality in wheat, Dencic et al. [5] showed that cultivar and environment interaction had significant effect on all quality traits. Also, variances of quality traits associated with genetic factors (cultivar) were generally larger than those for cultivar by environmental interaction effects.
In terms of white heads, transgenic lines derived from Khazar cultivar got placed in group 7 (relatively sensitive) which had white heads between 16-25%. Other cultivars stood in group 9 (sensitive); it means that the percentage of white heads of these cultivars was above 25%. By increasing the number of larvae in stem, their survival rate in the stem increases and yield decreases in the same way. Also, the reduction of tillers per hill increases the probability of stem damage. However, if there are several tillers, the tillers are not contaminated will produce the crop, and plants can recover from the damage caused by pests. By comparing three transgenic rice lines, Moghaieb [25] reported that on transgenic rice, no surviving larvae were recorded, and a larval mortality of up to 100% was found 4 days after infestation with stem borers. Conversely, on control rice cultivars, all the recovered larvae were observed healthy without any larval mortality. Moreover, the transgenic rice was very toxic to the rice stem borer (Chilo agamemnon) under laboratory condition [25].
The correlation of traits can be due to genes' connection or the existence of a genetic interaction with an environmental component. Ye et al. [40] also did not observe significant correlations between harvest index and grain yield among modern wheat varieties. With relatively little possibility for increases in grain yield by improving harvest index, greater yield potential must come from increases in net primary productivity [3]. Li et al. [22] revealed that it was important to increase net dry matter production over the entire growing season to achieve a high grain yield. It was consistent with the opinion of Donmez et al. [7].
Path coefficient analysis revealed that direct selection for number of fertile spikelet per panicle and harvest index would likely be effective for increasing grain yield [29]. Direct selection of number of effective tillers per plant, grain yield per plant, hundred grain weight, grain breadth, grain length, and grain thickness would increase harvest index. This study also indicated that there is no common causal factor directly influencing both grain yield and harvest index, though one hundred grains' weight, grain length, grain breadth, and grain thickness could be augmented in selection criteria for the simultaneous improvement of both the traits [29]. Path coefficient analysis in Gravois and Helms [12] studies revealed that panicle density had the largest positive direct effect on rice yield. Even under low rates of seeding, the number of filled grain per panicle increased, which compensated for the reduced panicles density. Direct effects for filled grain per panicle and grain weight were of secondary and/or tertiary importance in determining rice yield. Unfilled grain per panicle had negligible effects on rice yield. To achieve optimum rice yields and grain quality in a directseeded cultural system, adequate panicle density per unit area of uniform maturity must be achieved. Path analysis for paddy yield indicated that the number of spikelet per panicle and flag leaf length had positive direct effects and days to complete maturity, and plant height had negative direct effects on paddy yield under optimum irrigation condition, while flag leaf width and number of filled grains per panicle had positive direct effects, and days to 50% flowering had negative direct effect on paddy yield under drought stress condition [1].

Conclusion
According to combined variance analysis, there were significant differences among genotypes for all studied traits. According to results of mean comparison, the existence of minimum and maximum values for different traits in the same applied conditions for genotypes in each region indicates genetic variation. This means that by plant breeding, we can improve that trait. Therefore, this issue can provide a basis for determining the best cultivar according to the conditions of the region for cultivation. To determine the direct and indirect effects of the variables entered into the stepwise regression model, path analysis was used. The most positive direct effect on paddy yield was related to number of filled spikelet per panicle. The total direct and indirect effects of this trait on paddy yield (genetic correlation with paddy yield) were 93%. Regarding the direct effect of number of filled panicle per panicle on paddy yield and its indirect effect by day to 50% flowering and number of panicle per square meter, this trait can be a proper one for achieving higher yields, achieved by reducing the negative indirect effect of dead heart and white heads.