Growing Degree Days during the Late Reproductive Phase Determine Spike Density and Cognate Yield Traits

Abstract

Drought has become more frequent in central Asia causing large losses in cereal yield. To surmount the existing problem, it is imperative to emphasize early maturing varietal development. However, the impact of heat units on spike morphology and its relationship with yield potential is still unclear. Thus, the current investigation was carried out to test wheat lines and varieties for variation in total heat unit’s accretion for anthesis and maturity and to understand the manipulating impact of sunlight on spike morphology, grain yield and its cognate traits. Furthermore, the gene action controlling major traits inheritance, combining ability effects, heritability, and association studies were also estimated. Following the Half Sib/Full Sib approach 27 hybrids along with 12 parents were tested. Results depicted broad variation in genetic stock. Correlation study demonstrated that earliness negatively affects the yield, while positively influencing spike density. Genetic variances were greater than variances due to environment, pointing to higher heritability (>50%) for all the characters except for grain’s weight spike⁻¹. The degree of dominance revealed that the partial and over-dominant type of gene action conditioned inheritance of investigated traits. Thus, earliness can be used as an indirect selection criterion for yield advance.

1. Introduction

The primary concern of many researchers is about yield’s enhancement. During the past few decades, the improvement in production is not as obvious as it was half a century ago, while the increase in population is at its highest [1]. To keep rising wheat yield at a similar pace of population growth, we have to rely on genetic improvements more than ever [2]. No significant new growing areas will be incorporated in the future, nor will the contribution of management practices be noteworthy enough to improve yield, due to environmental and economic concerns. The apparent trend in current production showed a negligible increase in average yields. There is a need to employ alternative approaches for yield improvement. It could be possible either by nutrient supplementation in crop species (biofortification) [3] or by direct or indirect selection criterion for yield improvement. A majority of the crops are supplemented for health benefits with various elements like selenium, zinc, amylose and organic micronutrients [4], while the process of supplementation by means of exogenous fertilizers lead towards positive as well as negative impacts on crops physiology [5]. Environmental and abiotic factors have influential effects on yield and can serve as indirect selection criteria for its improvement. In this context, physiological basis provides help to identify opportunities for future breeding (using either conventional or molecular tools) and may help to break the apparent barriers to keep yield gains as healthy as they used to be during the 1960–1990s [6].

Harsh climatic conditions, especially the damaging impacts of abiotic factors (specifically drought) on the physiological response of plants, manipulate the normal functioning of crop species. In addition, genetic and environmental interactions greatly influence yield. The efforts to attain yield improvement of major field crops are always demanding. For a successful breeding program, the selection of parents is of immense importance. Based on the magnitude and nature of variation among the population and the association of connate characters, better parents can be selected [7]. Production can be enhanced by developing early maturing varieties, which must be productive and can be grown in different stresses and agro-climatic conditions. Selection for the improvement in grain yield can be most effective only when the genetic material displays variability [8]. A technique developed by Kempthorne is a powerful tool to use in pedigree selection, in order to assess combining ability estimates among parents and progenies [9]. The performances do not obligatorily predict the combining abilities of parents as a good or poor combiner. To overcome this arduousness, it is essential to amass knowledge about gene actions [10]. The information about heritability adds another dimension in assessing the natural response. Heritability is determined by the type of gene action and gives information about genetic variability. Hence, it is valuable to predict the selection response in the prospering generations [11]. Selection in early generations with desirable characteristics can be fruitful with higher heritability along with high genetic advance [12]. However, the association studies for earliness and heat units consumed by the crop to attain its physiological maturity, and how it affects the spike density and yield of the plant is still unclear.

Wheat served as an important everyday diet worldwide and belongs to family Poaceae, being the leading cereal worldwide and staple food in Pakistan. It has more nutritional value than any other food source, provides 55% of the carbohydrates and 20% of the calories of the world’s need annually [13]. Diversified climatic and abiotic factors certainly influence its yield. Hence, all the minor factors should also be taken into consideration for future use and are of similar value as major productive traits. However, our understanding of the influence of climatic factors on spike morphology and yield potential is limited, especially in wheat. The present investigation was directed to explore the nature and magnitude of gene action, heritability and association studies for growing degree days and yield cognate traits in Triticum aestivum L. and their possible interactive roles for improved production.

2. Materials and Methods

2.1. Experimental Conditions

The experiment was carried out at coordinates 31.4310° N, 73.0695° E during the 2013/2014 and 2014/2015 cropping season. To test the correlations and nature of gene action, material from our previous report was used [14]. Three male testers viz; E-108, E-113, E-114 (of International Maize and Wheat Improvement Center (CIMMYT) origin) represented broad genetic base and nine locale female parental lines (three varieties and six inbred lines) viz; Punjab-11 abbreviated as PB-11, AAS-11, AARI-11, 9859, 9860, 9861 and 9730, 9731, 9733 were crossed using line × tester mating design and 27 F₁ hybrids developed during the 2013/2014 cropping year. The seeds were sown in three replications using randomized complete block design (RCBD) during the 4th week of November of the 2014/2015 cropping season. The distance maintained for plant × plant and row × row was 15 cm and 30 cm respectively. The data presented for parents were the pooled data of two years (when parental material was used to develop crosses 2013/2014, and next year sown with crosses 2014/2015) and presented as mean average data of 2 years while the data presented for F₁ was 1-year triplicate mean average data. Moreover, triplicate experimental repeats of the material at different places made the results statistically more reliable.

2.2. Agronomic Practices and Data Collection

Data were collected for metric traits, namely: Days to 50% maturity, productive tillers plant⁻¹, ×1000-grain weight (g), and grain weight spike⁻¹ (g) measured in grams [15]. Growing degree days (GDD) (measured in terms of heat units from sowing to complete maturity of the crop) was calculated according to [16] with a little modification i.e.,

GDD = ((Max Temp − Min Temp of a day)/2) − base temperature of the crop,

(1)

Temperature was taken in degrees Celsius (°C) while base temperature of the crop was set as 5.

Spike density (SD) was measured by the formula:

SD = Fertile spikelets per main spike/Spike length (cm),

(2)

Genetic variances were computed, using general combining ability (GCA) and specific combining ability (SCA) values, as [17]

Additive genetic variance (б²_D) = 2 × б²_GCA,

(3)

Dominance genetic variance (б²_H) = б²_SCA,

(4)

б² = Variance, б²_GCA = Variance of GCA, б²_SCA = Variance of SCA.

Expected genetic advance (GA) was evaluated with one selection cycle at 10% selection intensity as [18]

$$GA=K\times \sqrt{{{\sigma }^2}_P}\times h^2,$$

(5)

K = selection differential, being 2.06 and 1.75 at 5% and 10% selection intensity, respectively.

$\sqrt{{{\sigma }^2}_P}$ = standard deviation of the phenotypic variance of the population under selection.

h² = heritability estimates in fraction of the trait under study.

2.3. Statistical Analysis

Data were expressed as means with least significant difference (LSD) in order to separate and compare the means, then subjected to analysis of variance (ANOVA) [19]. General combining ability (GCA) and specific combining ability (SCA) effects were determined as in an earlier report [20]. The t-test (2 and 1 tailed) at p ≤ 0.05 or 0.01 was applied to test the significance for correlation and combining ability estimates. Correlations and heritability analysis were performed using Agricolae package “R” version 3.4.2 [21]. Other statistical analyses were performed using Microsoft Excel 2016 and GenStat (10th statistical package) [22].

3. Results and Discussion

The analysis of variance for combining ability exposed significant differences among all traits for parents and hybrids. The mean values, combining ability effects, nature of genetic control for the inheritance of traits and the proportion of heritable change were accessed. However, correlation studies of earliness influencing spike density and yield were the major concern under study.

3.1. Estimation of Mean Square Values

Concomitant paramount differences were observed between mean square values. Treatment effects were highly significant for all traits under study at probability p ≤ 0.01, while non-significant differences were observed between replications. Parents and crosses depicted highly significant differences for all traits. While parents vs. crosses (interaction) revealed differences that were highly significant for ×1000-kernel weight and spike density. Adequate genetic variability was present in the material to assess combining ability effects (

Table 1. Mean square (MS) values from analysis of variance (ANOVA) for metric traits.

Source of Variation	d.f	DTM	GDD	Tp−1	SD	GWTsp−1	×1000 GWT
Replications	2	0.577 NS	1810.200 *	0.855 NS	0.0097776 NS	0.106 NS	2.943 NS
Treatments	38	32.710 **	8531.742 **	4.086 **	0.027416 **	0.196 **	46.924 **
Parents	11	23.174 **	12281.764 **	3.684 **	0.00302 **	0.202 **	58.984 **
Parents vs. crosses	1	0.145 NS	372.657 NS	0.201 NS	0.313389 *	0.398 *	337.790 **
Crosses	26	37.997 **	7259.005 **	4.405 **	0.017475 **	0.185 **	30.634 **
Lines	8	40.255 **	13944.306 **	4.566 **	0.007112 **	0.278 **	46.898 **
Testers	2	80.658 **	6959.387 **	0.044 NS	0.02436 **	0.093 NS	73.961 **
Lines × Testers	16	31.535 **	3953.807 **	4.870 **	0.036279 *	0.150 *	17.086 **
Error	76	0.349	406.072	0.712	0.006678	0.080	1.477

d.f, degree of freedom; DTM, Days to 50% maturity; GDD, Growing degree days (heat units); Tp⁻¹, Productive tillers per plant; SD, Spike density; GWTsp⁻¹, Grain weight per spike (g); ×1000 GWT, 1000-grain weight (g). Different values derived from ANOVA indicate significant differences at probability; ** = p ≤ 0.01; * = p ≤ 0.05; ^NS = Non-significant.

).

3.2. Study of Mean Values among Parents

Genetic variation and mean performance could be exploited for genotypic evaluation of the parents and hybrids. The average mean differences were 130.9 days, 1584.6 heat units, 8.69 tillers, 1.68, 2.83 g and 44.16 g for parameters like days to 50% maturity, growing degree days, tillers plant⁻¹, spike density, grain weight spike⁻¹ and ×1000-kernel weight respectively (Figure 1). Moreover, overall mean differences revealed that lines and hybrids were 2.15 days earlier and consumed 57.2 fewer heat units to meet their physiological maturity, while hybrids gained ×1000-kernel weight advantage of 3.47 g over parents. Parents along with hybrids depicted similar results for tillers plant⁻¹, spike density and weight of grain spike⁻¹. Among parents, line 9730 performed better with the least mean value (126.03 days) for days to 50% maturity and (1512.12) growing degree days. Hence, early maturing parents could be preferred to overcome the drought problem. Maximum tillers contributed as (10.6) by line 9731. Line 9859 showed a promising advantage for ×1000-grain weight (51.667 g) with the weight of grain spike⁻¹ (3.02 g). While AARI-11 was superior for spike density (1.88).

Unlike parents, the average mean differences for hybrids showed broad variation for almost all the traits: 13.66 days, 189.25 heat units, 5.44 tillers, grain weight spike⁻¹ (0.913 g) and ×1000-grain weight (11.93 g) (Figure 2). The minimum mean values were observed in crosses PB-11 × E-113 and 9731 × E-114 for days to 50% maturity and growing degree days respectively. However, mean performances are not a valid measure to assess variation between parents and hybrids. Screening of genetic stocks should be based on GCA/SCA effects, not just mean values.

3.3. Estimation of GCA and SCA Effects

Among lines, 9733 and PB-11 proved best for days to 50% maturity, growing degree days, tillers per plant and spike density. While testers E-108 and E-114 were best for growing degree days, spike density, grain weight and productive tillers plant⁻¹ [23]. The lines 9731, 9860 and 9861 showed high GCA effects for spike density and grain weight per spike (

Table 2. General combining ability (GCA) effects of parents.

Lines	DTM	GDD	Tp−1	SD	GWTsp−1	×1000 GWT
9730	−0.977 NS	−28.812 NS	−0.918 NS	1.245 *	0.324 NS	4.685 NS
9731	−2.177 NS	−57.585 NS	−0.696 NS	−0.283 NS	0.071 *	3.130 NS
9733	0.004 *	2.132 *	0.504 *	−0.070 NS	0.000 NS	−1.037 NS
9859	1.463 NS	−25.590 NS	1.340 NS	0.306 *	−0.194 NS	−0.704 NS
9860	2.352 NS	23.504 NS	−0.624 NS	0.447 *	0.180 *	−1.004 NS
9861	1.887 NS	45.438 NS	0.449 *	2.191 **	0.026 *	−1.093 NS
AARI-11	1.142 NS	−27.835 NS	−0.252 NS	0.508 *	−0.096 NS	−0.626 NS
AAS-11	−4.163 NS	64.660 NS	0.155 *	−1.640 *	−0.225 NS	−1.393 NS
PB-11	0.468 **	4.088 *	0.042 *	2.063 **	−0.087 NS	−1.959 NS
Testers
E-108	1.767 NS	13.956 **	0.043 *	1.597 *	−0.046 NS	−0.963 NS
E-113	−1.687 NS	−17.546 NS	−0.037 NS	1.448 *	−0.020 NS	−0.948 NS
E-114	−0.079 NS	3.590 *	−0.006 NS	2.222 *	0.066 *	1.911 NS
S.E for lines	0.197	6.717	0.281	0.794	0.094	0.405
S.E for testers	0.114	3.878	0.162	0.797	0.054	0.234

DTM, Days to 50% maturity; GDD, Growing degree days (heat units); Tp⁻¹, Productive tillers per plant; SD, Spike density; GWTsp⁻¹, Grain weight per spike (g); ×1000 GWT, 1000-grain weight (g); S.E = Standard Error; * and ** = significance at the 0.05 and 0.01 levels of probability, respectively (1-tailed); ^NS = Non-significant.

). Both poor and good combiners can contribute to the elevated performance of a specific cross combination, by crossing recessive and dominant alleles from them respectively [24].

Crosses which showed significant positive effects for kernel yield are listed in

Table 3. Specific combining ability (SCA) effects of F₁.

Hybrids	DTM	GDD	Tp−1	SD	GWTsp−1	×1000 GWT
9730 × E-108	−1.433 NS	−1.390 NS	0.625 *	−18.077 NS	0.171 *	0.789 *
9730 × E-113	1.021 NS	0.229 *	−0.295 NS	2.811 NS	0.091 *	0.215 *
9730 × E-114	0.413 *	1.160 *	−0.329 NS	1.153 *	−0.262 NS	−1.004 NS
9731 × E-108	−0.233 *	34.466 NS	−0.488 NS	1.522 **	0.211 *	−1.956 NS
9731 × E-113	1.221 NS	−5.265 NS	−1.408 NS	0.527 *	−0.089 NS	0.904 *
9731 × E-114	−0.987 NS	−29.201 NS	1.895 NS	−13.108 NS	−0.122 NS	1.052 *
9733 × E-108	−2.791 NS	−25.251 NS	−0.358 NS	−8.589 NS	−0.118 NS	−3.256 NS
9733 × E-113	2.700 NS	47.118 NS	0.882 **	0.297 *	−0.024 NS	−0.130 NS
9733 × E-114	0.091 *	−21.867 NS	−0.525 NS	4.840 NS	0.142 *	3.385 NS
9859 × E-108	−0.840 NS	−18.678 NS	−1.194 NS	1.425 **	0.088 *	1.911 NS
9859 × E-113	−2.313 NS	5.473 *	2.002 NS	1.939 NS	0.082 *	−0.896 NS
9859 × E-114	3.153 NS	13.205 *	−0.808 NS	3.440 NS	−0.171 NS	−1.015 NS
9860 × E-108	1.604 NS	2.194 *	1.660 NS	0.767 *	−0.232 NS	0.878 *
9860 × E-113	−2.721 NS	−14.754 NS	−0.480 NS	−3.878 NS	−0.004 NS	1.470 NS
9860 × E-114	1.117 NS	12.560 *	−1.180 NS	3.625 NS	0.236 *	−2.348 NS
9861 × E-108	−5.180 NS	−52.673 NS	−0.080 NS	84.000 NS	0.202 *	1.700 NS
9861 × E-113	0.487 *	−4.354 NS	−0.333 NS	−9.800 NS	−0.024 NS	0.026 *
9861 × E-114	4.693 NS	57.027 NS	0.413 *	−29.697	−0.178 NS	−1.726 NS
AARI-11 × E-108	2.558 NS	20.599 **	0.515 *	0.735 *	0.011 *	−1.467 NS
AARI-11 × E-113	1.155 NS	7.718 *	−1.409 NS	2.047 NS	−0.315 NS	−2.974 NS
AARI-11 × E-114	−3.713 NS	−28.317 **	0.894 **	1.454 **	0.305 **	4.441 *
PB-11 × E-108	0.860 NS	2.522 *	−1.452 NS	−13.591 NS	−0.034 NS	1.867 NS
PB-11 × E-113	−0.126 NS	26.073 NS	0.965 **	−0.149 NS	0.160 *	−0.441 NS
PB-11 × E-114	−0.734 NS	−28.595 NS	0.487 *	0.850 *	−0.126 NS	−1.426 NS
AAS-11 × E-108	5.456 NS	38.210 NS	0.771 *	−1.545 NS	−0.298 NS	−0.467 NS
AAS-11 × E-113	−1.424 NS	−62.238 NS	0.075 *	−2.097 NS	0.122 *	1.826 NS
AAS-11 × E-114	−4.032 NS	24.027 NS	−0.846 NS	−1.108 NS	0.176 *	−1.359 NS
S.E for crosses	0.341	11.634	0.487	0.796	0.163	0.702

DTM, Days to 50% maturity; GDD, Growing degree days (heat units); Tp⁻¹, Productive tillers per plant; SD, Spike density; GWTsp⁻¹, Grain weight per spike (g); ×1000 GWT, 1000-grain weight (g); * and ** = significance at the 0.05 and 0.01 levels of probability, respectively (1-tailed); ^NS = Non-significant.

. The combination AARI-11 × E-114 holds potential to be used as the best hybrid for the traits: Weight of grain, ×1000-kernel weight, tiller plant⁻¹, growing degree days and spike density. Similar significant SCA effects were observed in cross combination AARI-11 × E-108 in the desired direction for productive tillers, GDD, spike density and tillers plant⁻¹. 9731 × E-108 and AARI-11 × E-114 exhibited significant negative SCA effects for days to 50% maturity and growing degree days. Although positive SCA effects were also observed in them for yield traits. A similar trend in cross combinations 9730 × E-113, 9860 × E-108, 9860 × E-114 and 9861 × E-113 was observed for earliness and ×1000-kernel weight. The trend fluctuated among different cross combinations for the traits under study (Table 3). The parents with good GCA can also be used to develop a pure line with the yield improvement due to the additive type of gene even though these cross combinations depicted non-significant SCA effects [25].

3.4. Proportional Contribution of Line and Testers

Maternal influence was predominant for traits like growing degree days (59.10%) and ×1000 kernel weight (47.10%) [26,27], while maternal × paternal interaction was predominant for days to 50% maturity (51.073%), tillers/plant (68.031%), spike density (45.353%) and grain weight/spike (49.935%). The paternal influence was not so obvious for most of the traits (Figure 3). The results depicted that maternal and maternal × paternal interaction contribute more towards genetic variation of cognate traits [27].

3.5. Genetic Components and Degree of Dominance

GCA variance was lower than the SCA variance for all the mentioned traits (

Table 4. Genetic component variations (additive, dominance genetic effects).

Genetic Variation	Days to 50% Maturity	Growing Degree Days	Productive Tillers per Plant	Spike Density	Grain Weight per Spike	×1000-Grain Weight
Variance of GCA	3.336	933.060	0.213	1.002	0.017	3.880
Variance of SCA	10.395	1182.578	1.386	7.122	0.024	5.203
Additive variance	6.671	1866.121	0.427	1.002	0.034	7.761
Dominance variance	10.395	1182.578	1.386	7.122	0.024	5.203
Variance ratio of GCA to SCA	0.321	0.789	0.154	0.141	0.714	0.746
Degree of dominance	1.248	0.796	1.802	2.666	0.837	0.819

). These findings are favored by the ratio of GCA/SCA variance which was smaller than unity [17]. Therefore, it is perceivable that the dominant gene action conditioned parameters inheritance [28]. The non-additive genetic variance for grain yield plant⁻¹ and other cognate characters was also pointed out previously [29,30]. The differences in the results were mainly due to breeding material and the diversified genotype and environment interactions. Among all the studied traits, over-dominance was observed for tillers/plant, spike density and days to 50% maturity. The dominance genetic effects were observed for grain weight, ×1000-kernel weight and duration of the vegetative growth period (growing degree days for earliness), revealing that selection of superior genotypes in an F₁ generation could be useful for producing productive hybrids.

3.6. Heritability, Genetic Gain and Correlation

Heritability and expected genetic advance are mentioned in

Table 5. Heritability, genetic gain and coefficient of variability.

	DTM	GDD	Tp−1	SD	GWTsp−1	×1000 GWT
Ve	0.3	406.1	0.7	0.0067	0.1	1.5
Vg	10.787	2708.557	1.124667	0.006907	0.038667	15.149
Vp	11.136	3114.629	1.836667	0.013585	0.118667	16.626
H2	0.96866	0.869624	0.612341	0.508431	0.325843	0.911163
GA	5.689168	85.41762	1.460566	0.104299	0.197554	6.538871
CV%	0.5	1.3	9.7	5.0	10	2.7

V_e, Environmental variance; V_g, Genotypic variance; V_p, Phenotypic variance; H², Heritability; GA, Genetic advance; CV, Coefficient of variability; DTM, Days to 50% maturity; GDD, Growing degree days (heat units); Tp⁻¹, Productive tillers per plant, SD, Spike density; GWTsp⁻¹, Grain weight per spike (g); ×1000 GWT, 1000-grain weight (g).

. It was found that a major proportion of variability in phenotype was due to genotypic variation and not environmental variation. High heritability was found (>50%) for all traits while grain weight spike⁻¹ (32.5%) was moderately heritable [31]. The heritability values fluctuate between moderate to highly heritable i.e., 32.5% to 96.8% for different parameters [32] shadowed by genetic advance ranged from 0.19 g (grain weight spike) to 85.4 heat units (growing degree days). Hence, beneficial highly heritable characters could be recovered in the very next generations.

It is possible to test whether the direct effects of different photoperiods on the length of spike also translate into grain number differences. At days to 50% and 100% maturity, the number of fertile florets and days taken by genetic stock to attain its maturity were counted, which reflects closely the number of grains in most conditions. The dynamics of floret development was significantly affected by the photoperiod during the spike length elongation period. The longer the length of spike associated with shorter photoperiod, the higher the spikelet fertility will be because distal and less developed floret primordia were able to progress to the fertile floret stage. These findings favoured by correlation studies were the same as previously reported [33]. Early maturing crops can best fit in double pattern cropping season, with benefits over moisture use and avoiding delay seasonal effects, insects and pest damage. That may allow wheat to flourish and best fit with ever-changing demand. Moreover, chemical and pesticides used as fertilizer and irrigation applications could be minimized. Hence, reduction in the maturity time of crop can bring ultimate benefits to cope with ever-increasing challenges [34]. A significant positive relationship among DTM, spike density and growing degree days was found (

Table 6. Genotypic and Phenotypic correlations.

Phenotypic Correlation	Genotypic Correlation
	DTM	GDD	Tp−1	SD	GWTsp−1	×1000 GWT
DTM	1 **	0.45125 **	0.08029 NS	0.2902 **	−0.19064 *	−0.0225 NS
GDD	0.40744 **	1 **	0.13473 NS	0.25388 **	−0.18587 *	−0.28906 **
Tp−1	0.05706 NS	0.13118 NS	1 **	−0.1084 NS	−0.31496 **	−0.02665 NS
SD	0.20942 *	0.19091 *	−0.17791 *	1 **	−0.46559 **	−0.01083 NS
GWTsp−1	−0.10889 NS	−0.08238 NS	−0.16573 NS	−0.24462 **	1 **	0.2847 **
×1000 GWT	−0.01657 NS	−0.27514 **	−0.01835 NS	−0.01083 NS	0.20208 **	1 **

DTM, Days to 50% maturity; GDD, Growing degree days (heat units); Tp⁻¹, Productive tillers per plant, SD, Spike density; GWTsp⁻¹, Grain weight per spike (g); ×1000 GWT, 1000-grain weight (g). * and ** = Correlation is significant at the 0.05 and 0.01 level of probability, respectively (2-tailed); ^NS = Non-significant.

, r_g = 0.25, r_p = 0.19, p < 0.01). The assessments of correlation clearly demonstrated strong and negative association of earliness (DTM, GDD) with ×1000 grain weight (r_g = −0.289, r_p = −0.275) and grain weight/spike (r_g = −0.1906). In this presentation, we attempted to envisage, from published and recent unpublished evidence, using studies carried out under both controlled and field conditions. The rate of crop development was manipulated during the late reproductive phase. The later the crop is harvested, the more negative its yield will be. Thus, harvesting at the right time results in the most fruitful output.

In the case of the field study, changes in heat unit duration also related to changes in the weight of grains set in both parents and hybrids. Averaging across years, grain weight per spike was reduced by exposing the plants at days to 50% and 100% maturity. This could be one of the reasons for low yield in genetic stock under consideration. Moreover, earliness greatly influences the spike density positively. The days required by the crop to meet 50% maturity and the total amount of heat units consumed by crop while reaching its maturity greatly affect spike density. The association was positive and significant at phenotypic and genotypic levels (Table 6). The growing degree days positively regulate the spikelets density, however the environmental changes at the grains filling stage may lead to an earlier onset of physiological maturity resulting in lower transfer of nutrients and essential elements in grains, which ultimately affects yield.

The studies in which the plants were subjected to different photoperiods throughout development also tended to show that the longer the duration of the heat units consumed by crop, the higher the number of grains per spike. However, the increase in the weight of grain was not as obvious as the number of grains, which ultimately affect yield. Although these cases confirm the results from experiments with manipulations of photoperiods focussed on spike density alone, they also may reflect ‘memorised effects’ of the photoperiod [35] from previous phases.

4. Conclusions

From the aforementioned discussion, it can be seen that the late reproductive phase manipulated the wheat yield as a strong and negative association seen. Majorly, parameters were conditioned by the dominant and over-dominant type of gene interaction. Combination AARI-11 × E-114 proved to be the best hybrid, consuming lesser amounts of heat units while reaching its maturity and at the same time had maximum positive SCA estimates for spike density, grain weight and ×1000-kernel weight favored by association studies. Among parents, E-114 and PB-11 revealed significant GCA estimates for yield traits. High heritability was observed for all traits except for the weight of grain spike⁻¹, which was moderately heritable escorted by genetic gain. Hybrid breeding is recommended for improving the quality traits as the dominance variance was predominant for traits under study. These evidences provide experimental support, with plants grown in the field, that sensitivity to photoperiod may actually be used as an indirect tool to further rise wheat yield. This opens the possibility to manipulate the sensitivity to photoperiod during spike elongation and fertile spikelet’s formation (a mirror action of modifying the photoperiods) as an alternative avenue for wheat breeders to increase yield potential.