Introduction
Wheat, being the second most important staple food of the human population, is adapted to a wide range of eco-climatic conditions on six continents of the world. The demand for wheat is rising in tandem with the ongoing population growth and is predicted to increase by 40% by 2030. Thus, in order to maintain long-term food security, wheat production must be drastically increased [1]. Low yields can be caused by a variety of factors, such as low-quality seed, improper weed removal techniques, delayed planting, poor soil management, uneven fertilizer application, diseases, a lack of water, heat and drought stress brought on by climate change, and more [2]. Wheat’s status is crucial among cereal crops due to its high consumption and nutritional value. Wheat scientists now face more problems in breeding wheat varieties with increased yield, quality, and resilience to both biotic and abiotic stressors due to the world’s population growth and improved lifestyle. A diet high in wheat contains more fiber than a diet high in meat [3]. However, research has demonstrated that wheat proteins lack important amino acids like lysine and threonine [4]. Developing new and improved varieties of wheat that can yield more and perform better under various agro-climatic conditions might lead to an increase in wheat output and quality [5]. Drought is one of the many obstacles compromising wheat production. Low rainfall and unpredictable fluctuations in precipitation are the main causes of drought [6]. Water shortages reportedly cause 17% to 70% of yield losses. 50% to 90% less wheat was produced than could have been under irrigation in emerging nations as a result of a water shortage. During the stages of tillering, jointing, booting, anthesis, and filling, wheat plants exhibit a significant response to water deficiency stress. The crucial stage of tillering is when the wheat plant produces its florets, spikelets, primodia, and tillers [7]. At this point, water deficit stress can result in a 46% reduction in wheat production overall. The process of drought resistance is complicated and depends on a number of parameters, including crop type, drought intensity, duration, and plant development stages. Under drought stress, a plant must rely on many systems at once to survive [8]. Recurrent drought is the leading abiotic stress causing reduced production and productivity of wheat, especially in low altitude regions [9].
Vigorous seedlings are a key indicator of a plant’s ability to produce large amounts of food quickly [8, 10]. The amount of chlorophyll in leaves serves as a proxy for the photosynthetic capacity of plant tissues. In situations of drought, the concentration of chlorophyll pigments changes. Plant resilience to drought stress is dynamically influenced by carotenoid levels. Conditions lacking in water prevent the synthesis of chlorophyll a and b, lower the quantity of binding-related proteins, and decrease the amount of pigment-related proteins that are necessary for light harvesting and other proteins related to photosystem II. Various techniques were developed to screen wheat genotypes for drought tolerance during the seedling stage. Relative water content and the root-to-shoot ratio were proposed by [11] as selection criteria for drought resistance in wheat. Relative water content (RWC) is a useful factor for choosing drought-tolerant wheat types when the seedling stage is approaching. The relationship between RWC and cell volume can accurately depict the equilibrium between water absorbed by the plant and water released by transpiration [12].
There are still relatively few cultivars with stable drought tolerance expressed in a variety of environmental conditions, despite the fact that breeding for drought tolerance is acknowledged widely as an important strategy. The identification of superior genotypes and environmental conditions that impact its expression complicate the polygenic characteristic of drought tolerance. Additionally, the identification of superior genotypes under fluctuating moisture conditions is complicated by the unexpected nature of drought manifestation and the various methods that plants have chosen to cope with drought stress. For drought tolerance, field-based empirical selection is still frequently employed [13, 14].
Therefore, considering the aforementioned facts, the current investigation aimed to predict the performance of wheat lines/varieties for grain yield and yield-related parameters under normal and drought-related stress conditions.
Materials and Methods
Experimental Location and Site
The experiment was carried out in the field of the Department of Plant Breeding and Genetics (PBG), Islamia University of Bahawalpur. In this work, 50 genotypes were grown and assessed for physio-morphological characteristics at maturity under both normal and drought-stressed conditions using a Randomized Complete Block Design (RCBD) in a triplicate fashion. The analysis revealed that the soil was sandy clay loam with a pH of 7.6 and 0.85% organic matter. The seeds of 50 wheat genotypes were collected and arranged for sowing purposes. In the beginning, seeds of each genotype were sown in a 2 m long row for each replication under normal and drought conditions. Plant-to-plant and row-to-row distances were maintained at 6 inches and 12 inches, respectively. The experiment included standard agronomic and cultural procedures, with two seeds per hill being sown with the aid of a dibbler. To maintain the plant population, seedlings were thinned to one seedling per hill after they were established. In the normal experiment, recommended irrigation was applied at three critical stages, i.e., at (1) tillering (35 days after sowing, DAS), (2) booting (85 DAS), and (3) milking (112 DAS [15]). In the stressed experiment, drought stress was applied at the tillering stage by upholding (missing) the irrigation treatment. One set of genotypes was irrigated at all three critical stages, while the other set of the same wheat genotype was kept under drought stress, missing the irrigation at the first (tillering) critical stage at 35 DAS. At the time of maturity, from each replication, 5 healthy plants randomly selected for each genotype were grown under normal and drought experimental conditions. The data was collected from tagged plants for yield and yield related traits, namely: No of tillers, Spikelets per spike, Chlorophyll Content, Flag Leaf Area (cm), Plant Height (cm), Peduncle Length (cm), Spike Length (cm), Biomass (g), Main Spike Weight (g), Grain Yield per Spike (g), Yield per plant (g), and 1000 grain weight (g). The analysis of variance (ANOVA) was performed on the scoring data of variance (ANOVA) using statistics 8.1 [16]. Genotypic and phenotypic correlation analyses were conducted using R Studio. For correlation analysis significance levels, = 0.01 was used for highly significant effects and =0.05 was used for significant effects. Principle component analysis (PCA) was performed using the Minitab software to ascertain the relationship between traits and genotypes under drought-stricken conditions [17]. The following formula was used to calculate the reduction percentage:
Reduction percentage = (normal traits mean-drought trait mean)/normal trait mean) *100
Results and Discussion
Analysis of Variance
Analysis of variance revealed the existence of highly significant differences among different wheat genotypes for grain yield and quality-related traits under both normal and drought conditions (Table 1). All the studied characteristics, including NT, SPS, CC, PH, PL, FLA, BM, MSW, MSGW, YPP, and TGW, exhibit highly significant variations among wheat genotypes, according to the conducted analysis. The accessions that showed variance in performance in comparison to the genotypes investigated for the researched qualities under drought conditions were indicated by the mean variability of different seedling traits; those with the best performance were thought to be drought-tolerant (Table 1). In this study, the relationship between genotype and phenotype varies significantly among different genotypes. There were identical results published by Fouad et al. [18].
The mean values of spikelets per spike changed significantly, ranging from 3.14 to 4.66 under drought. The genotypes G45 (3.14), G7 (4.16), G24 (4.27), G38 (4.43), and G16 (4.66) performed best under drought conditions, and the genotypes G49 (16.07), G11 (16.95), G41 (16.96), G23 (17.43), and G32 (19.35) performed poorly under drought conditions, as shown in Table 2. Under drought duress, the characteristics of the shoot and root significantly declined. The decline in shoot characteristics was significantly greater than that of root characteristics. In addition, under drought stress, shoot length, fresh weight, and dry weight all decreased. Jager et. Al. also [19] also reported identical outcomes. One essential characteristic of wheat plants that is affected by water scarcity is their shoot length. Any trait’s phenotypic response results from the interaction of environment and genotype [20]. Under normal circumstances, the mean values for shoot length ranged from 3.08 to 3.98 cm, while under conditions of stress, their dimensions were 11.00 to 25.19 cm.
Under drought circumstances, the genotypes that showed the highest quality of biomass were G45 (-8.04), G24 (-1.27), G16 (3.85), G38 (4.30), and G7 (4.34), while the low quality of biomass was shown by the genotypes G11 (20.29), G32 (21.25), G49 (23.52), G41 (25.03), and G23 (33.03), as shown in Table 2. According to [21], all genotypes exhibited highly significant differences in wheat yield, plant height, and thousand grain weight. Wheat breeders [22] published identical results for agromorphological characteristics and cereal yield under stress and non-stress conditions. The genotypes that showed superlative performance for YPP were G7 (4.46), G16 (9.51), G24 (9.67), G45 (10.19), and G38 (10.22). On the other hand, the genotypes that showed poor performance were G11 (23.11), G49 (23.56), G41 (24.83), G23 (28.50), and G32 (30.89), as shown in Table 2. G24 (13.81), G7 (14.49), G16 (14.93), G45 (14.94), and G38 (14.94) are those genotypes that showed the best under drought conditions for main spike weight, and G24 (1.69), G7 (2.91), G38 (3.25), G16 (3.29), and G45 (3.39) showed the best for 1000 grain weight.
Table 1.
Analysis of variance (ANOVA) for grain yield and related traits under normal and drought stress conditions.
Table 2.
Reduction percentage under drought stress condition.
Table 3.
Best and worst genotypes under drought stress conditions.
Photosynthesis is the primary process in plant cells that controls the low concentration of water culture media [23]. More chlorophyll will result in more effective photosynthesis mechanisms. The genotypes G7 (-53.15), G16 (-28.65), G38 (3.56), G24 (4.65), and G45 (5.31) showed the best performance under drought conditions, and the genotypes G23 (20.55), G11 (20.75), G41 (23, 82), G49 (25.63), and G32 (26.00) performed poorly under drought conditions, as shown in Table 3. Similar results were observed by wheat scientists [24, 25], who reported that drought stress substantially affects chlorophyll and causes leaves to turn yellow rather than green.
Correlation Analysis
The correlation analysis indicated strong positive correlations for most of the traits within the same experiment. In this experiment, evidence of both genotypic and phenotypic correlations of different traits in normal and stress conditions may help advance strategies for the assortment of required varieties with preferred traits (Table 4). According to genotypic correlation, positive and highly significant correlations were found between the number of spikelets per spike (r = 0.56**) with plant height, main spike weight, chlorophyll content, spike length, yield per plant, and thousand grain weights [18]. In addition to a positive and highly significant genotypic correlation with thousand-grain weight (r = 0.70**), there was a positive and statistically significant genotypic correlation with peduncle length, flag leaf area, biomass, main spike grain weight, and yield per plant. Yield per plant (r = 0.70**) showed a positive and significant phenotypic association with thousand-grain weight, as shown in Table 4. However, a positive and statistically significant genotypic connection was found between the area of the flag leaf and the chlorophyll content. Plant biomass (r=0.88**) was positively and statistically significantly associated with main spike weight, main spike grain weight, and plant yields. Javed et al. [26] also found a strong positive link between biomass output, plant height, 1000-grain weight, and grain yield in their investigation of advanced wheat genotypes. The thousand-grain weight (r = 0.06NS), flag leaf area, and peduncle length were shown to have a positive but not statistically significant phenotypic correlation under drought (Table 5). According to Ali et al. [27], there is a strong positive correlation between plant height and grain weight. In the current study, the relationship between genotype and phenotype varies significantly for each genotype. Main spike grain weight (-0.15ns) was found to have a negative and minor phenotypic connection to thousand-grain weight, but a positive and substantial phenotypic connection to yield per plant. Both Javed et al. [28] and Baye et al. [26] showed similar results under drought conditions, which are negative but statistically significant.
Principal Component Analysis
Principal component analysis (PCA) was used to examine the wheat genotypes to characterize the variety of the germplasm and the relationship between the wheat seedling indices under both normal and drought stress conditions. PCA based on the correlation matrix. The statistically significant principal components (PCs) were selected based on the eigenvalues’ importance, as determined by Kaiser [29]. The PCs classified as significant were limited to those with eigenvalues greater than one. Multivariate statistical analysis, known as principal component (PC) analysis, is used to examine and deconstruct big and complicated datasets. The results of Kamel [12, 30] showed that this analysis reduces the number of associated variables from a bigger one to a smaller one.
Table 4.
Genotypic correlation and phenotypic correlations under normal conditions.
Of the twelve main components (Table 6), four showed substantial fluctuation during drought circumstances, and the first five showed Eigenvalues of more than one (significant) under normal conditions. Under normal and stress conditions, the first three PCs showed 0.75% and 0.72% total variation, respectively, in the studied germplasm. Under normal circumstances, the first PC explained 0.24% of the variation, the second, 0.20%, the third, 0.11%, and the fifth, 0.9%. Under drought conditions, as shown in Table 6, the first PC contributed 0.21% of the total variance, the second contributed 0.20%, the third contributed 0.11%, the fourth contributed 0.10%, and the fifth contributed 0.10%.
The PC1 was highly related to spike length (0.48) and main spike weight (0.46) under normal conditions, while under stress conditions It was highly related to peduncle length and chlorophyll content (Table 6). There is a negative association between PC4 and SPS, which suggests that SPS has a major impact on this factor in typical settings. On the other hand, PC2, PC4, and PC5 were found to have unfavorable associations. The fact that PL is positively correlated with both PC3 and PC4 demonstrates that it has an effect on those variables. The fact that YPP is positively correlated with both PC2 and PC4 demonstrates its effect on those subsystems. TGW had a positive connection with PC5 that was statistically significant. Table 7 showed that TGW has a small effect on PC1, PC2, PC3, and PC4 under typical settings, as indicated by the small correlations between these variables.
The PCA extracted twelve principal components (PCs) to assess drought tolerance in the wheat genotypes based on their key morphological and yield related traits (Table 6). Among these 12 PCs, only five had an eigenvalue greater than 1, having variations up to 0.72% of the overall data under drought stress conditions (Table 6 and Fig. 1).
This (Table 7) is a principal component analysis loading plot (PCA). Principle component analysis (PCA) is a statistical method that creates new variables known as principle components (PCs) from the original variables in order to minimize the dimensionality of a data collection. The majority of the variance in the data is captured by the PCs, which are simply linear combinations of the original variables. The loading plot displays the relative importance of each original variable on each PC. The graph has two plots, one for the data under drought conditions and one for the data under normal conditions. Each plot has two axes, PC1 and PC2, which are the first and second principal components. The lines in the plots represent the original variables, such as P1, P2, P3, etc. The length and direction of each line indicate how much each variable contributes to t he PCs. We can see that P1 has a higher loading on PC1 under drought conditions than under normal conditions, which means that P1 explains more variation in the data when there is drought, as shown in Fig. 2.
Table 5.
Genotypic correlation and phenotypic correlation under drought stress conditions.
Table 6.
Variability under normal and drought conditions.
Table 7.
Loading factor of significant variables under normal drought conditions.
When choosing parents from a variety of backgrounds for hybridization and other plant breeding techniques, principal component analysis can be useful [31]. The selection of the various parent groups benefited from the genotype projection on PC1 and PC2. The expected pattern of genotypes on the two PCs indicated the population structure under normal and drought conditions. The genotypes appearing in the same square box in Figs. 3a) and b), had the same performance, while those appearing in different squared boxes showed differing performance. Under normal conditions, G5, G11, G26, and 38 were opposite to G49 and G50. Genotypes G3 and G44 were opposite to genotype G10, while G33 and G17 were opposite to each other and showed clear diversity among all genotypes (Fig. 3a). There was a clear difference between drought-tolerant and drought-susceptible genotypes. Under drought conditions, the genotypes G10 and G37 were opposite to G16, G19, and G32, while G50 was opposite to G11 and G24. Genotype G44 and genotype G10 showed clear diversity from G16 and G20 (Fig. 3b).
Fig. 2.
(A) Principal component (PC) loading plot (projection of variables) of Fig. 2a Principal component (PC) loading plot of seedling traits under normal conditions. (B) Principal component loading plot of seedling traits under drought conditions.

Fig. 3.
(A) Two-dimensional score plot of wheat genotypes on PC1 and PC2 under normal conditions. (B) Two-dimensional score plot of wheat genotypes on PC1 and PC2 under drought conditions.

None of the interventions resulted in statistically significant alterations in the PCA. It is a powerful statistical method for reducing the number of dimensions in a dataset and extracting actionable, fact-based feedback from a complex one [32]. The angle between the cosines of the trait vectors is used in the PCA Biplot to approximately depict the correlation between traits. An angle less than 90 degrees indicates independence between the features, whereas an angle greater than 90 degrees indicates a positive link. Our findings, however, demonstrated unequivocally that the contributions of a trait pair to the Biplot for PCA and the correlations between a trait pair were well-coordinated with the approximate vector angles (Figs. 3, 4). Other researchers have utilized PCA Biplot analysis to successfully screen drought-tolerant wheat cultivars [33]. So it is likely that it is useful here as well.
Conclusion
In this study, a total of 50 wheat genotypes were used to investigate the response of each genotype to drought stress. There are two different sets of plots (one with normal conditions and the other with drought stress). According to the results of this study, five genotypes (G7, G16, G24, G38, and G45) performed best in both normal and drought conditions. The genotypes G11, G23, G32, G41, and G49 did not perform better under both stressful and non-stress conditions due to their high decline in reduction percentage. Among these spring wheat genotypes, those that did better under both conditions were deemed drought-tolerant, and those that performed the worst were deemed drought-susceptible. Under normal conditions, PC1 exhibited 0.24%, whereas PC1 exhibited 0.21% under drought conditions. In order to develop wheat germplasm that is tolerant to drought stress, this study found favorable features under drought stress that might be used in wheat breeding programs. As a result, by utilizing contemporary breeding techniques, these genotypes will be essential in the future development of drought tolerance.

