Response of Indian Dwarf Wheat and Persian Wheat to Sowing Density and Hydrothermal Conditions of the Growing Seasons

: The need for foods with high nutritional value has led to the rediscovery of ancient wheat species Triticum sphaerococcum and T. persicum as raw materials with valuable consumption properties, but their reintroduction requires assessment of their productivity under different agricultural practices. The field experiments were carried out for three years (2018–2020) to test the hypothesis that the sowing density of T. sphaerococcum and T. persicum (400, 500 and 600 no m − 2 ) will affect their agronomic traits, yield and occurrence of diseases, but the response will depend on the hydrothermal conditions of the growing seasons. In this study, a significant correlation of the grain yield with the amount of precipitation in tillering, and from booting to the beginning of fruit development was demonstrated. The sowing density of T. sphaerococcum had an impact on the grain yield only under moderate drought stress during the growing season (2019), when the highest yield was obtained at a sowing density of 600 m − 2 . In 2019, the most favorable sowing density was also the highest for T. persicum. In the year with the lowest amount of rainfall during the growing season (2018), the yield of T. persicum was the highest in the lowest sowing density. At the shooting stage, a greater intensity of powdery mildew was observed on T. persicum , especially with higher sowing densities. Increasing the sowing density also increased the occurrence of root rot symptoms in both wheat species in the year that favored the occurrence of this disease (2018). It can be concluded that in the integrated low-input cultivation of T. sphaerococcum and T. persicum, it is justified to use a sowing density of 600 pcs. m − 2 , in an agroclimatic zone with moderate droughts during the growing season.


Introduction
One of the key objectives of long-term EU policy is to increase the production of foods with high nutritional value using low-input cultivation technologies [1]. Ancient wheat species characterized by an increased content of macro and micronutrients [2] and much higher (as compared to other cultivated wheat species) contents of phenolic acids and alkylresorcinols, with proven pro-health effects [2], may be useful for these purposes. Examples of ancient species of cereals with increased nutritional value are Indian dwarf wheat (Triticum sphaerococcum Perc. = T. aestivum L. ssp. sphaerococcum (Perc.) Mac Key) and Persian wheat (Triticum persicum Vav. = T. carthlicum Nevski) [3,4]. The grains of both genotypes are a rich source of phenolic acids and alkylresorcinols. Moreover, in these wheats, an enhanced share of ferulic acid was found, accompanied by higher quantitative and qualitative variability of homologues within sterols, tocols and carotenoids [3]. Indian The incidence of fungal diseases is also a known factor significantly affecting the size of wheat crops [26][27][28]. Among them, the ones most often observed on the leaves of wheat plants are powdery mildew (caused by Blumeria graminis), Septoria leaf spot (caused by Zymoseptoria tritici, synonym Septoria tritici) and wheat leaf rust (caused by Puccinia recondita). Foot and root rot diseases, in properly cultivated spring forms of wheat, do not pose a great threat. However, the symptoms of root rot (caused by Gaeumannomyces graminis, Fusarium spp., Rhizoctonia spp.), Fusarium foot rot (caused by Fusarium spp.) and eyespot (caused by Oculimacula acuformis, and Oculimacula yallundae) can be observed during the growing season [29,30]. Diseases such as powdery mildew, Fusarium foot rot and root rot are of particular importance in the cultivation of Persian wheat and Indian dwarf wheat [31]. The incidence of disease symptoms in wheat significantly depends on weather conditions. Rainfall, not only the total, but also the frequency, distribution in particular stages of plant development, and moisture in the canopy play a special role here [32]. For pathogens whose spores are easily spread by raindrops splashing from the lower infected parts of the plant to the upper leaves, the closer the leaves are to each other, the easier the spores are spread to the upper parts [33]. The incidence of disease symptoms in wheat crops also depends on various agrotechnical practices, including the date and density of sowing [28,34]. This effect may be direct, but often the sowing method indirectly influences the incidence of diseases on plants [35].
The aim of the research was to evaluate the response of yield, biometric traits and incidence of diseases on the ancient wheat species, Indian dwarf and Persian wheat, using integrated, low-input (limited N rate and pesticide consumption) technology depending on the sowing density under varied hydrothermal conditions of the growing seasons. It was assumed that 1) increasing the sowing density would increase the number of production shoots per area unit and increase the grain yield, but at the same time would contribute to greater disease pressure; 2) the response of yield, plant biometric traits and the incidence of diseases will depend on the hydrothermal conditions during the growing season.

Site Description and Crop Management
Field experiments were located in Mochełek, Kuyavian-Pomeranian voivodeship (53°13′ N; 17°51′ E), Poland. The soils at experimental sites were characterized as Alfisol (USDA). The abundance of available macronutrients, pH and C-org content in the site of the research is presented in Table 1. The precipitation and thermal conditions in the subsequent years of field experiments were very diverse ( Table 2). The year 2018 was warm (except March) with a large amount of rainfall in April and July. The hottest month in 2019 was June, with heavy rainfall only in May. 2020 was cool, with average rainfall in May and heavy rainfall in June and July. At the beginning of August, immediately after harvesting of the previous crop (triticale), the catch crop of a pea was sown. This catch crop was plowed at the end of November, after the above-ground part had frozen over. At the beginning of March, in each year of the research, pre-sowing mineral fertilization was performed at the following rates: 30 P2O5 kg ha −1 (superphosphate), 50 K2O kg ha −1 (potassium salt) and 30 N kg ha −1 (ammonium sulphate) and additionally, 30 N kg ha −1 at the beginning of the shooting stage. In early spring (middle of march), the soil was cultivated with a tillage unit. Sowing of the wheat was performed 9 April 2018, 1 April 2019 and 23 March 2020 using OYORD plot seeder, in a row spacing of 16 cm. In determining the seeding rate, the 1000-grain weight and germination capacity as well as plant establishment (the ratio of seedlings vs. seeds), a value of 90% was set. Persian wheat was harvested at the end of July, and Indian dwarf wheat 10 days later, using a Wintersteiger plot harvester.
In our study, the cultivar 'Trispa' of Indian dwarf wheat was used [36]. This cultivar has long (65-80 cm), stiff shoots. The stems are covered with wax coating. The spike is short (5 cm) and awned. The spikelets are one-or two-grained. The grain of Indian dwarf wheat is rounded and relatively small (average 1000-grain weight is 30 g). For Persian wheat, the cultivar 'Persa' was used [36]. This cultivar develops thin and long stems (65-80 cm), which are susceptible to lodging. There is no wax coating on the leaves. The medium long spike (7-8 cm) is brownish, loose, and awned. The spikelets are two-or threeflowered, and one-or two-grained. The grain of Persian wheat is elongated and small (average 1000-grain weight is 28 g). Both Persian wheat and Indian dwarf wheat grains are naked.

Measurements of Morphological and Agronomic Traits
At the heading stage, the number of generative and vegetative tillers per plant was determined on subsequent plants collected from 1 linear meter from each plot. At full maturity, the number of generative fertile tillers (with grains) and sterile generative tillers (with spikes without grains) was determined on test plots of 1m 2 on each plot. Moreover, 50 spikes were randomly collected from each plot, on which the number of grains per spike were measured. During harvest, the grain and straw yields as well as their moisture content were determined. Then, yields were converted into a constant humidity of 15%. The 1000-grain weight was determined based on 500 pcs, in four replications of each treatment.

Diseases Occurrence Assessment
The plant health assessment included an analysis of the severity of the most dangerous leaf diseases, the stem base and spring wheat roots, i.e., powdery mildew, Fusarium foot rot and root rot. The severity of disease symptoms caused by mildew was determined at stem elongation (Se) and at fruit development (Df). On the first assessment date, carried out prior to the fungicide application, the health of the lower leaves (L3 and L4) was assessed, and on the second date, the health of the flag leaf (L1), which was described as a percentage of the leaf area that was covered with Blumeria graminis mycelium or disease symptoms caused by this pathogen, was assessed. On the second assessment date, the severity of disease symptoms was also determined based on damage to wheat stalks caused by Fusarium fungi (Fusarium foot rot) and the severity of disease symptoms on the roots caused by a complex of pathogens, i.e., Gaeumannomyces graminis, Fusarium spp., Rhizoctonia spp. (root rot). The assessment of the severity of powdery mildew was performed directly on the experimental plots and the severity of Fusarium foot rot and root rot was determined in the laboratory. One hundred generative tillers of each combination were assessed for health status and the occurrence or severity of disease symptoms. In laboratory conditions, after separating of the tillers and removal of the leaf sheaths, the degree of infection of the stem base and roots was assessed on a scale of 0-4, where 0 meant no symptoms of a given disease, and 4 meant severe infection, when the symptoms of Fusarium foot rot covered the stem base around followed by its rot or symptoms of root rot were observed on more than 60% of the surface of the assessed roots. Infection rates were transferred to the disease index (DI) according to the Townsend and Heuberger formula and expressed as a percentage.

Statistical Analyses
The basic statistical descriptors included arithmetic mean ( ̅ ) and standard deviation (SD). Analysis of the results was performed using statistical inference methods [37]. Normality was examined by the Kolmogorov-Smirnov test and homogeneity of variance with the Levene's test. The data were analyzed by two-way analysis of variance (ANOVA), in which the sowing density and study years (temperature and precipitation were included as passive variables) were the main factors, and then to determine significant differences between the means, the Tukey test was used. Based on the same data set, the analysis of the structure of agronomic traits and disease occurrence, depending on the years of study and sowing density, was performed via Canonical Correspondence Analysis (CCA). Pearson's correlation coefficient was used to measure the correlation between grain yield, biometric features of T. sphaerococcum and T. persicum and Sielianinov's hydrothermal coefficient. Quadratic polynomial regression was used to relate the grain yield and sowing density of T. sphaerococcum and T. persicum in the study years 2018-2020. The analyses were performed using PAST 3.2 software, CANOCO 4.5, MS EXCEL 365 and STATISTICA 13.3 (Hammer UiO, 2018, Microcomputer Power, USA, Microsoft 2020, Statsoft 2019). The significance level for all analyses was assumed to be minimal, α = 0.05.

Yield and Biometric Features of Triticum sphaerococcum
The grain yield of Indian dwarf wheat varied over the years of the research and ranged from 1.60 to 4.79 Mg ha −1 ( Table 3). The least favorable year for the growth and development of wheat was 2018, where low rainfall in May and June resulted in a strong reduction of the yield ( Table 2). The most favorable, in turn, was 2020, where heavy rainfall in May and June was beneficial for the yield. In 2019, there was heavy rainfall only in May. Under these conditions, the grain yield was average (2.50 Mg ha −1 ). A significant effect of the sowing density on the grain yield of Indian dwarf wheat was demonstrated only in 2019. Increasing the sowing density from 400 or 500 pcs. m −2 to 600 resulted in a significant increase in the yield (Table 3, Figure 1A). There was no significant effect of sowing density on the number of grains per spike in any of the study years. In contrast, 1000-grain weight in 2019 was higher at the sowing density of 400 pcs. m −2 compared to the size of this parameter at the sowing density of 500 pcs. m −2 . Indian dwarf wheat developed the lowest straw weight in 2018 and the highest in 2020. The straw weight was higher at the sowing density of 600 pcs. m −2 compared to the straw weight at other sowing densities only in 2019. The variation coefficients of the grain and straw yield and the number of grains per spike of Indian dwarf wheat were highest in 2018. In 2018, Indian dwarf wheat developed the least generative tillers and the most vegetative tillers per plant (Table 4). In this unfavorable year for plant growth and development, the number of generative tillers per plant was higher at lower sowing densities (400 and 500 pcs. m −2 ). Only in 2019, at a sowing density of 600 pcs. m −2 , the number of vegetative tillers per plant was lower compared to the number of such tillers at the sowing density of 400 pcs. m −2 . Moreover, in 2019, the number of generative tillers per m 2 and the number of total generative tillers were higher at the sowing density of 600 pcs. m −2 compared to that recorded in the sowing density of 400 pcs. M −2 . In 2019, the number of sterile generative tillers per m −2 and their share in total generative tillers was the highest. Among the analyzed biometric features of Indian dwarf wheat, the number of sterile tillers per m 2 was characterized by the highest variation coefficient.

Yield and Biometric Features of Triticum persicum
The highest grain yield of Persian wheat was obtained in 2020 (3.45 Mg ha −1 ) ( Table  5). It was 2.4 times higher than the grain yield in 2019. In 2018, the most unfavorable for plant growth, the average grain yield of Persian wheat was only 0.41 Mg ha −1 . In 2018, the grain yield was higher at the sowing density of 400 pcs. M −2 compared to that obtained in the density of 600 pcs. M −2 . In contrast, in 2019, the highest grain yield was obtained at a density of 600 pcs. M −2 , the yield was significantly lower at the sowing density of 500 pcs. M −2 and the lowest at a sowing density of 400 pcs. M −2 (Table 5, Figure 1B). The number of grains per spike was the highest in 2020 and the lowest in 2018. In 2019, increasing the sowing density from 500 to 600 pcs. M −2 had a positive effect on the number of grains per spike of Persian wheat. The Persian wheat grain obtained in 2018 had the highest 1000grain weight. There was no significant effect of Persian wheat sowing density on 1000grain weight or straw weight found in any of the study years. In 2018, Persian wheat developed the least generative tillers and the most vegetative tillers per plant ( Table 6). In the same year, the number of sterile generative tillers per m 2 was higher at the sowing density of 600 pcs. m −2 compared with the number of such tillers at the sowing density of 400 pcs. m −2 . In 2019, it was shown that the number of generative tillers per plant, as well as the number of such tillers per m 2 , was the highest at the highest sowing density of 600 pcs. m −2 . The coefficients of variation of the grain yield and most of the analyzed biometric features were higher in 2018 compared to 2019 and 2020. Only the coefficients of variation of the number of vegetative tillers per plant and sterile generative tillers per m were the highest in 2020.

Occurrence of Diseases
On Indian dwarf wheat, the severity of powdery mildew varied over the years of the study and observation dates (Table 7). On the first assessment date, most symptoms of the disease on the lower leaves (L3 and L4) were found in 2019 (45.58-56.67% of the area of assessed leaves with symptoms of the disease), where the high total rainfall in May was conducive to the development of this disease (Table 2). However, the least were found in 2018 (1.0-1.75% of the leaf area), where rainfall in May was the lowest. On the second assessment date, the severity of powdery mildew symptoms on the flag leaf (L1) was small and was observed only in 2019 and 2020. The occurrence of powdery mildew did not depend on the sowing density. The large amount of rainfall in May, during the period of intensive growth and development of plants, limited the occurrence of roots and stem base diseases. Hence, the fewest symptoms of root rot were observed in 2019 and the most in 2018 (Table 7). The severity of these disease symptoms, with the exception of 2019, significantly depended on the sowing density. In the year with the greatest number of symptoms of root rot (2018), the severity of the disease increased with an increase in sowing density. In 2020, the lowest number of disease symptoms was recorded at the sowing density of 500 pcs. m −2 , while at the sowing densities of 400 and 600 pcs. m −2 there was a significant increase in infection. Relatively few symptoms of Fusarium foot rot were recorded on Indian dwarf wheat, and in 2019 they were observed only sporadically. Only in 2020, a significant effect of the sowing density on the severity of the symptoms of Fusarium foot rot was demonstrated. Increasing the sowing density of Indian dwarf wheat from 400 or 500 pcs. m −2 to 600 pcs. m −2 resulted in a significant increase in the infestation of the stem base by fungi of the genus Fusarium. Many symptoms of powdery mildew were observed on Persian wheat, especially before the fungicide was applied. Their intensity depended on the year of the study, the observation date and the sowing density (Table 8). At the first assessment date, on average, the greatest number of disease symptoms caused by B. graminis was found in 2019 (34.92-58.75%), which was characterized by the highest total rainfall in May. There were slightly fewer of them in 2020 (41.26-49.42%), and clearly the fewest in 2018 (2.75-7.50%). In all years of observation, the most symptoms of powdery mildew were found at the highest sowing density (600 pcs. m −2 ). In 2019, many symptoms of the disease were also found at the sowing density of 500 pcs. m −2 . At the second assessment date, the severity of powdery mildew symptoms on Persian wheat plants was low. It was observed only in 2019 and 2020 and was not dependent on the sowing density. The large amount of rainfall in May 2019, during the period of intensive plant growth and development, limited the occurrence of symptoms of foot and root rot diseases on Persian wheat plants. Hence, the fewest symptoms of Fusarium foot rot were found in 2019, when they were observed only occasionally, and the most were observed in 2018. The incidence of this disease significantly depended on the sowing density. In 2019, the most symptoms of Fusarium foot rot were observed at the sowing density of 500 pcs. m −2 , and in 2020, at the sowing density of 400 pcs. m −2 . In general, the highest sowing density was not conducive to the occurrence of this disease. Furthermore, the symptoms of root rot were least observed in 2019 and the most in 2018, when the disease index value was high and ranged from 50.25 to 63.25% (Table 7). Despite such a high intensity of disease symptoms, there was no significant effect of sowing density on the severity of symptoms of root rot.

Variable Dependency Analysis
Among the grain yield structure components of Indian dwarf wheat and Persian wheat, biometric features such as the number of grains per spike and the number of fertile generative tillers per m 2 were significantly positively correlated with the grain yield (Table  9). These yield components and the grain yield were significantly positively correlated with the Sielianinov's hydrothermal coefficient in the first ten days of May, and with the values of this coefficient in June (Tables S1 and S2). The number of sterile generative tillers per m 2 was negatively correlated with the grain yield and Sielianinov's hydrothermal coefficient in June. Both studied wheat species showed a significant positive correlation between grain and straw yield. Canonical correspondence analysis (CCA) of variables from individual research years, including grain and straw yield, number of grains per spike, 1000-grain weight, number of generative tillers per plant, fertile generative tillers per m 2 and the incidence of diseases (powdery mildew, root rot and Fusarium foot rot), clearly indicates a very strong relationship with the research year (weather conditions in the growing season) and, to a lesser extent, with the sowing density ( Figure 2). Biometric features, yield of wheat and incidence of diseases are distributed horizontally, according to vectors: (longer) to the left -Year and (shorter) to the right -Sowing densities (the vectors-green lines -almost coincide with the line of the first axis). Axis 1 explains the system at 100% and is statistically significant, while the second axis is statistically insignificant. According to the conducted analysis, the significance of weather conditions (years of research) on shaping the analyzed features of the studied wheat species was highest in 2019, mainly at the lowest sowing density (400 pcs. m −2 ). A separate group is formed by 2018, where the effect of sowing density on the analyzed features is clearly visible, while the species is also important (Triticum sphaerococcum and T. persicum are grouped separately). In 2020, T. sphaerococcum also forms a separate group.  Canonical correspondence analysis of variables, on the average of three years of research, indicates that powdery mildew in flag leaf development significantly affects the formation of biometric features, more in the case of T. sphaerococcum, especially at the sowing densities of 400 and 500 pcs. m −2 (Figure 3). Powdery mildew in steam elongation, in turn, as well as root rot and Fusarium foot rot affect the biometric features of T. persicum to a greater extent, especially at the densities of 600 and 500 pcs. m −2 . Axis 1 of the CCA analysis statistically significantly explains the system to a degree over 90%, while the second axis is statistically insignificant.

Morphological and Agronomic Traits
Canonical correspondence analysis of grain yield and agronomical traits (number of grains per spike, 1000-grain weight, number of generative tillers per plant and fertile generative tillers per m 2 ) of Indian dwarf wheat and Persian wheat indicates a significant relationship with the year of the study (weather conditions in the growing season) and, to a lesser extent, with the sowing density. Auškalnienė and Auškalnis [10] indicate the greater importance of weather conditions compared to the sowing density of common wheat. Tokatlidis [38] concludes that the optimal wheat population is primarily affected by the water regime, while sowing date, heat, terminal drought, frost and type of soil are also indicators.
The least favorable year for the growth and yield of Indian dwarf wheat and Persian wheat was 2018, with low rainfall at the end of April, in May and June (from the stem elongation stage to the beginning of the grain development stage). Under severe water deficit, both wheat species exhibited reduced fertile generative tillers per m 2 , number of grains per spike, and the number of generative tillers per plant, morphological traits significantly correlated with the grain yield. Moreover, the water deficit and high air temperatures in May limited tillering, which resulted in a reduction in the number of generative tillers per plant.
According to Zhang et al. [23], drought decreases wheat grain yield the most during the complete growth cycle, which is achieved from a combination of decreases during the tillering stage (−27.4%), the stem elongation stage (−21.4%), the heading stage (−16.8%), the flowering stage (−17.7%), and the development of grain stage (−16.3%). In the study of Triticum aestivum it was proven that water stress affects morpho-physiological traits such as plant height, number of tillers per plant, number of spikelets per spike, number of grains per spike, 1000-grain weight, flag leaf area, relative water content, and chlorophyll content (d). Moreover, drought or heat stress enhances enzymatic antioxidant activity (catalase, dismutase, peroxidase), and proline content. These chemicals help to stabilize membranes, sub-cellular structures and cellular redox potential by destroying the free radicals [21]. In a study of the ancient wheat Triticum boeoticum, some proteins (MICP, in the group of lipid metabolism proteins) were found [24]. These proteins may play an important role in the metabolic pathways of wheat and its drought tolerance [24].
In our study, the number of generative tillers per plant of Indian dwarf wheat in 2018 was significantly higher at lower sowing densities (400 and 500 pcs. M −2 ). Similarly, Silva et al. [12], in the study of the response of common wheat cultivars to sowing density (75, 150, 225, 300, 375 pcs. m −2 ), showed an increase in the number of generative tillers per plant along with a decrease in sowing density. In the study by Beavers at al. [16] the number of spikes per plant was significantly larger at a sowing density of 330 pcs. m −2 compared to the number found at densities of 412 and 495, and the smallest at a sowing density of 660 pcs. m −2 . In 2018, the grain yield of Persian wheat was higher at the lowest sowing density (400 pcs. m −2 ) compared to that obtained in the sowing density of 600 pcs. m −2 . At the lowest sowing density, the number of sterile generative tillers per m −2 was significantly lower and the number of fertile generative tillers per m 2 was slightly higher, compared to the number of such tillers at a sowing density of 600 pcs. m −2 . It is worth noting that, at the lowest sowing density, the share of sterile generative tillers in the total of generative tillers was lower (12%) compared to the share of such tillers at sowing densities of 500 and 600 pcs. m −2 (20 and 21%, respectively). The reason for the increase in the number of sterile spikes at higher sowing densities was the increased competition for water under conditions of heavy shortage, which started in May and deepened in June. The study by Pinheiro et al. [39] clearly shows an increase in the number of non-productive tillers as the sowing density of wheat increases. Marinho et al. [40] prove that competition between plants and tillers can increase the incidence of infertile tillers, especially at high sowing densities.
The year 2019 was characterized by a moderate water deficit. Heavy rainfall occurred in the second half of May (at the stem elongation stage), while a rainfall deficit and high air temperatures were recorded in June (from the booting to the beginning of the grain development stage). According to Sattar et al. [21], heat stress along with low water availability at the reproductive stage are major contributing factors towards lower wheat production. In 2019, increasing the sowing density from 400 or 500 pcs. m −2 to 600 pcs. m −2 resulted in a significant increase in the grain yield of Indian dwarf wheat. Similarly, in Persian wheat, the highest grain yield was obtained at the highest sowing density (600 pcs. m −2 ), while it was significantly lower at a density of 500 pcs. m −2 and was lowest at the density of 400 pcs. m −2 . The increase in grain yield at the highest sowing density was mostly due to the increased number of fertile generative tillers per m 2 , the most important yield component. Moreover, there was an increase in generative tillers per plant in Persian wheat and a reduction in number of vegetative tillers per plant in Indian dwarf wheat. In the study by Hussain et al. [13] among the yield components, only the number of generative tillers per m −2 increased in proportion to the increase in the grain yield at higher seed rates. According to Wu et al. [41], high sowing density as a purposeful agricultural practice increases yields, reduces weed infestation and increases sustainability by increasing biomass and soil cover during the growing season of the plants, which reduces nutrient loss and erosion. In 2019, there was a significant increase in the straw weight of Indian dwarf wheat at the sowing density of 600 pcs. m −2 . This weight consists mainly of generative tillers, hence the significant correlation of the straw yield with the grain yield. Furthermore, we have noted that the 1000-grain weight was higher at the lowest sowing density of 400 pcs. m −2 compared to the size of this parameter at a sowing density of 500 pcs. M −2 . Poutala et al. [9] showed a reduction in the 1000-grain weight of common wheat due to an increase in sowing density. In the study by Tahir et al. [8], a higher grain weight in common wheat was found at seed rates of 100 kg ha −1 and lower at 175 kg ha −1 .
In Persian wheat grown in 2019, despite the increase in the number of fertile generative tillers per m 2 at the sowing density of 600 pcs. m −2 , the number of grains per spike in this treatment was higher compared to wheat in a sowing density of 500 pcs. m −2 . Opposite relationships are presented by Tahir et al. [8], who showed a reduction in the number of grains per spike as a result of increasing the density of common wheat sowing from 125 to 175 kg ha −1 . The different response of the compared wheat species could result from a different structure of the spike. Common wheat in the studies by Tahir et al. [8] had a number of spikelets per spike twice as high as Persian wheat in 2019, which was significantly influenced by drought stress during flowering, heading and the beginning of grain development stages.
In the present study, the most favorable year for the growth of Indian dwarf wheat and Persian wheat was 2020, where heavy rainfalls in May and June favored the yield. The distribution of rainfall favored the development of numerous fertile generative tillers per m 2 and grains per spike -yield components significantly correlated with the grain yield. It is noteworthy that under these weather conditions, almost all generative tillers of Indian dwarf wheat and Persian wheat remained fertile (the shares of fertile generative tillers as a percentage of total generative tillers per m 2 were 99% and 98%, respectively).
In that year, no significant effect of the sowing density on any of yield components or the grain yield of both studied wheat species was demonstrated. The obtained results confirm the conclusion of Fisher et al. [11] on the lack of response of common wheat grain yield to the sowing density under favorable weather conditions.
Our study did not show a significant effect of sowing density on the number of grains per spike of Indian dwarf wheat in any of the study years. However, a correlation was observed between the number of generative tillers per m 2 and the number of grains per spike. If, in a given year, the number of generative tillers per m 2 was low, the number of grains per spike was also low. Therefore, there was no compensation of the yield components. The opposite conclusions were formulated by Lloveras et al. [42] who claimed that common wheat can compensate for low populations by modifying the number of grains per spike. However, in our research, the number of grains per spike in 2018 and 2019 was mostly limited by drought from the booting stage to the beginning of the grain development stage, which made it impossible to compensate for the limited density of generative shoots in the canopy.

Occurrence of Diseases
The analyzed Indian dwarf wheat and Persian wheat plants showed a similar intensity of disease symptoms as in previous studies on those wheat species conducted in organic farming crops, in which powdery mildew, Fusarium foot rot and root rot were also mainly observed [31]. On the first observation date, a high intensity of powdery mildew symptoms was observed on the assessed lower leaves of both wheat species, and its incidence was significantly correlated with the total rainfall. According to Te Beest et al. [43], the disease is favored by warm and humid conditions, but even under favorable conditions for the development of the disease, appropriate chemical protection can significantly reduce the disease's severity. Hence, on the first observation date, the fewest symptoms were observed in 2018, with low rainfall at the end of April and in May ( Table 2). The lesser severity of powdery mildew symptoms on the second observation date, found on the analyzed flag leaves of wheat, was influenced by the use of the Tango Star fungicide at the fully developed flag leaf stage. According to Semaškienė et al. [44], a mixture of fenpropimorph and epoxiconazole, active substances contained in this preparation, are highly effective in limiting the development of powdery mildew on spring wheat.
In the present authors' research, the occurrence of powdery mildew on Persian wheat also depended on the sowing density, which was observed only at the stem elongation stage. Generally, with increasing sowing density, the severity of the symptoms of this disease increased. This is consistent with the observations made on the form of winter common wheat, as an increase in plant density strongly contributed to the increased occurrence of powdery mildew symptoms [34,45]. It is known that dense sowing brings the leaves closer to each other, which facilitates the spread of pathogen spores. In addition, the sowing density may also indirectly affect the occurrence of pathogens and severity of diseases by changing the microclimate inside the canopy. In dense crops, the temperature tends to be more uniform, the air humidity is higher, and water remains on the leaves for longer, favoring infection by B. graminis [32,35]. Moreover, the distance between the infected and healthy tissues is also of great importance. The closer the leaves are to each other, the easier it is for the fungus spores to splash, and thus for pathogens to spread through the canopy [33]. On the other hand, Finckh et al. [46] found that, in the case of susceptible wheat cultivars, the severity of powdery mildew symptoms increased with the decrease in sowing density. They did not observe such an effect for resistant wheat genotypes. Such relationships were also not observed in previous studies conducted in organic crops, where in both Persian wheat and Indian dwarf wheat the sowing density had no effect on the severity of powdery mildew [31].
In the present authors' research, the severity of the symptoms of foot and root rot diseases significantly depended on the sowing density. On Indian dwarf wheat, the most symptoms of root rot and Fusarium foot rot were found in the highest sowing density of 600 pcs. m −2 . The sowing density of Persian wheat, however, had an ambiguous effect on the occurrence of Fusarium foot rot, and usually the fewest disease symptoms were found at the highest sowing density. Moreover, in this wheat, there was no significant effect of sowing density on the severity of root rot disease. Other authors, studying the winter forms of common wheat, also observed an increase in the severity of the symptoms of diseases observed on the roots and stem base at a higher sowing density [47,48], which was especially noticeable at the early development stages of the analyzed plants. They found that in the case of foot and root rot diseases, an increased plant density means a shorter distance between the host plant and the inoculum source, which increases the probability of the pathogen reaching the adjacent plant through the spores or mycelium and infecting it. In previous studies, when Indian dwarf wheat or Persian wheat were grown on organic farms, no significant influence of their sowing density on the severity of foot and root rot diseases, including the severity of Fusarium foot rot and root rot, was found [31]. In addition, Eken et al. [49] found that there was no significant effect of the sowing density on the severity of root and crown rot caused by Fusarium fungi in common wheat. In organic crops, due to the lower intensity of cultivation, the growth and development of plants is less intense, the canopy architecture is different, and thus the influence of factors on the development of pathogens is different as well [26,32].

Conclusions
Our research presents, for the first time, the characteristics of morphological and agronomic features and the intensity of disease symptoms in Indian dwarf wheat and Persian ancient wheat species depending on the sowing density under cultivation using lowinput technology. The three-year field experiments demonstrated that the grain yield and occurrence of diseases depended primarily on the weather conditions in the growing season and to a lesser extent on sowing density. Increasing the sowing density to 600 pcs. m −2 increased the number of fertile generative tillers per m −2 and the grain yield in the year in which moderate drought stress occurred during the growing season. However, the increased sowing density also contributed to greater disease pressure. In low-input integrated cultivation, Indian dwarf wheat and Persian wheat sowing should be performed at a sowing density of 600 pcs. m −2 , in an agroclimatic zone with moderate droughts during the growing season. Further research is needed on the physiological and biochemical responses of ancient species of cereals such as Indian dwarf wheat and Persian wheat on sowing density under various agroclimatic conditions. Supplementary Materials: The following are available online at www.mdpi.com/article/10.3390/ag-riculture12020205/s1, Table S1: Significant correlation coefficients (p ≤ 0.05) between Sielianinov's hydrothermal coefficient and grain yield, and biometric features of Triticum sphaerococcum, Table S2