The Inﬂuence of Parental Heat-Stress Priming on Drought-Tolerant Maize Progenies’ Field Performance

: Maize ( Zea mays L.) is a staple crop and an industrial crop. Improving its thermotolerance will be a crucial step in ensuring food security. The objective of this research was to assess the inﬂuence of the parental growth environment on their progenies in a non-heat stress environment. The progenies evaluated in this ﬁeld study were obtained from a 2 × 3 × 3 factorial in a completely randomised greenhouse experiment. Two contrasting growth environments, three maize varieties, and three soil amendments were used in the greenhouse study. A randomised complete block design experiment with three replications was used to grow the progenies. The progenies were examined for nineteen morphological attributes. In this study, 69.51% of the yield variation was explained by the ﬁrst and second principal component axes. Among the studied attributes, grain weight and cob weight explained more variations in the progenies than the other attributes. The interaction of the parental heat-stress and soil amendment conditions elicited different responses from the drought-tolerant maize progenies. Based on the differences in their yield attributes, the progenies were grouped as poor yielders (Cluster IV), good yielders (Cluster I) and high yielders (Clusters II and III). The parental growth environment inﬂuenced the progenies’ ﬁeld performance in a non-heat-stress environment. Further evaluation of the progenies under a heat-stress environment and molecular analyses are required to establish that a transgenerational effect has occurred. S.M.


Introduction
In the field, plants are exposed to a wide range of stressors at various stages of development. Plants, as sedentary organisms, respond to stress by modifying their physiological and genomic architecture [1]. These adjustments frequently result in decreased resource utilisation efficiency, which can result in yield losses ranging from minor to severe losses, or, in the worst-case scenario, plant death [2][3][4].
Plants use phenotypic alteration as a stress-resistance mechanism. Depending on the type of stress, the severity of the stress, and the developmental phase of the plant during which the stress was encountered, the plant's developmental process can be altered in anticipation of future occurrences of the stress [5][6][7][8]. Stress memory refers to the ability to remember stressful events. Stress memory can last anywhere from hours to months, and it can be passed down to a stressed plant's progeny in the hopes that the progeny will be exposed to stress in the future [6,7]. The parental effect, interorganismal effect, or transgenerational effect refers to the inheritance and expression of stress memory in plants that have not been exposed to the stress (primed) [9,10]. Stress memory inheritance In the 2019/2020 and 2020/2021 summer planting seasons, field evaluations were conducted at the North-West University Experimental Farm (−25.7902166, 25.6187922 3.5 mi) Mafikeng campus, North West Province, South Africa. The experimental soil used for the greenhouse study and field evaluation was composed of 82% sand, 4% silt, and 14% clay (loamy sand), had a pH (KCl) of 4.98, 11 mg P kg −1 (Bray1), 290 mg K kg −1 , 390 mg Ca kg −1 , 163 mg Mg kg −1 , 5 mg Na kg −1 , 376 mg kg −1 total N, and 0.37% organic C [21]. The South African Weather Services provided weather data for the experimental site during the two planting seasons (Figure 1).

Experimental Design, Treatments and Cultural Practices 2.2.1. Parental Evaluation
The treatments evaluated in this field study were developed from a 2 × 3 × 3 factorial in a completely randomised greenhouse experiment [3,21]. The greenhouse experiment is described briefly. The greenhouse study included two contrasting growth environments (heat stress [HS] environment and non-heat stress [NHS] environment), three mediummaturing drought-tolerant maize varieties (WE3128, WE5323, and ZM1523), and three soil amendments (poultry manure application [PM], mineral fertiliser application [MF], and combination of poultry manure/mineral fertiliser (50:50) application [MPM]). The weekly Agriculture 2021, 11, 1229 3 of 15 temperature range for the NHS was from 15 to 33.9 • C and from 16.2 to 32.5 • C for the first and second planting seasons, respectively. In the HS, the weekly temperature ranges for the first and second planting seasons were from 16.38 to 40.88 • C and from 19.83 to 44 • C. Recommended greenhouse maize production protocols were followed. The cobs were harvested and shelled when they reached physiological maturity. The seeds were packaged and stored as per the eighteen treatment interactions.

Parental Evaluation
The treatments evaluated in this field study were developed from a 2 × 3 × 3 factorial in a completely randomised greenhouse experiment [3,21]. The greenhouse experiment is described briefly. The greenhouse study included two contrasting growth environments (heat stress [HS] environment and non-heat stress [NHS] environment), three mediummaturing drought-tolerant maize varieties (WE3128, WE5323, and ZM1523), and three soil amendments (poultry manure application [PM], mineral fertiliser application [MF], and combination of poultry manure/mineral fertiliser (50:50) application [MPM]). The weekly temperature range for the NHS was from 15 to 33.9 °C and from 16.2 to 32.5 °C for the first and second planting seasons, respectively. In the HS, the weekly temperature ranges for the first and second planting seasons were from 16.38 to 40.88 °C and from 19.83 to 44 °C. Recommended greenhouse maize production protocols were followed. The cobs were harvested and shelled when they reached physiological maturity. The seeds were packaged and stored as per the eighteen treatment interactions.

Progeny Field Evaluation
The seeds were evaluated in the field using a randomised complete block design with three replications to assess the impact of the different parental growth conditions on their progenies. The average minimum and maximum temperatures during the field evaluation for 2019/20 planting season were 13.89 °C and 29.78 °C, respectively, while for 2020/21 planting season, the average minimum and maximum temperatures were 17.27 °C and 31.81 °C, respectively.
The land was mechanically ploughed and harrowed, and plots (235 cm by 100 cm) were manually marked. Both 100 cm and 50 cm pathways were maintained between and within blocks, respectively. Three replications (blocks) were created, each with eighteen plots for each treatment interaction. On each plot, three maize seeds were planted per hole

Progeny Field Evaluation
The seeds were evaluated in the field using a randomised complete block design with three replications to assess the impact of the different parental growth conditions on their progenies. The average minimum and maximum temperatures during the field evaluation for 2019/20 planting season were 13.89 • C and 29.78 • C, respectively, while for 2020/21 planting season, the average minimum and maximum temperatures were 17.27 • C and 31.81 • C, respectively.
The land was mechanically ploughed and harrowed, and plots (235 cm by 100 cm) were manually marked. Both 100 cm and 50 cm pathways were maintained between and within blocks, respectively. Three replications (blocks) were created, each with eighteen plots for each treatment interaction. On each plot, three maize seeds were planted per hole at 75 cm by 30 cm inter-and intra-row spacing, according to treatment interactions. Four weeks after planting (WAP), the seedlings were thinned to one plant per hole. A total of sixteen plants were established per plot, resulting in 68,085 plant population per hectare. Mineral fertiliser was applied at a rate of 350 kg N ha −1 . At planting, inorganic NPK grade 13:7:10 (30) + 0.5% Zn + 5% S + 3% Ca supplied 60% of the N requirement, while the remaining 40% was top-dressed with lime ammonium nitrate (28% N) at seven weeks after planting. The plots were manually kept weed free, and irrigation was provided by a sprinkler.

Data Collection
The emergence percentage (EP) was calculated by counting the number of seedlings that emerged in each plot at two WAP over the number of planted seeds. Growth and yield data were collected from the inner four plants in each plot, while phenological data was collected from all the plants in a plot. Plant height (PH) was measured from the soil surface to the tip of the last leaf using a 5 m steel measuring tape. Leaf length and width were measured with a measuring tape from the base of the leaf on the maize stalk to the tip and at the widest part of the leaf. The leaf area (LA) was calculated by multiplying the product of leaf length and width by a constant of 0.75 [22]. The number of leaves (NoL) per plant was counted and the leaf chlorophyll content (LCC) was determined using a portable chlorophyll content metre (CCM-200 plus, Opti-Sciences, Hudson, NY, USA). The stem diameter (SD) was measured using a digital Vernier calliper (153-006-11, ACCUD, Vienna, Austria). The growth traits were measured at 4, 7, and 10 WAP.
Plants were monitored at the start of the reproductive phase to determine the number of days until tassel appearance (DT 50 ) and silk appearance (DS 50 ) in half of the plants in a plot. The tassel silk interval (TSI) was calculated by subtracting the DT 50 from the DS 50 .
The plants were harvested with sharp secateurs when they reached physiological maturity. The number of cobs per plant (CN) and columns per cob (CC) were counted. The cobs were shade-dried to a moisture content of 12%, while the shoots were oven-dried to a constant weight at 65 • C to determine the stover dry weight in grams (SDWt). A measuring tape and a digital Vernier calliper were used to determine the length of the cob (CL) and width of the cob (CW). Individual cob was weighed to determine the gram weight of the cob (CWt). For the grain weight determination, each cob was shelled, and the grains weighed. The obtained CWt and GWt were used in calculating the shelling percentage (SP) as: Numigral seed counter (CHOPIN Technologies, Villeneuve-la-Garenne, France) was used to count the grains to obtain the grain number (GN). Thereafter, the 100-seed weight (SWt) was obtained by weighing 100 seeds from each treatment interaction.

Data Analysis
For the purposes of analysis, the data from the two field evaluations were combined. Analysis of variance (ANOVA) was performed on the data, following the procedure for randomized complete block design using GenStat Software (VSN Int. Ltd., Hemel Hempstead, UK). Protected Fisher' s least significant difference test was used as a post hoc test for mean separation at p ≤ 0.05. Bar charts were plotted using a Microsoft Excel spreadsheet (Microsoft Office 365, Microsoft, Redmond, WA, USA) to show the trend of measured attributes. A principal component biplot was created using GenStat Software (VSN Int. Ltd., Hemel Hempstead, UK) based on the yield attributes of the treatment interactions.

Analysis of Variance
The SD, DT 50 , CW, and CC were all significantly affected by the growth environment, while LA, PH, SD, DT 50 , CW, and SWt were all significantly influenced by the maize variety ( Table 1) Table 1 shows the growth characteristics measured at 10 WAP.

Effect of Parental Heat Stress on the Progenies' Growth and Yield Attributes
Despite being statistically similar, NHS provided higher quantitative EP and LA at 10 WAP, NoL at 4, 7, and 10 WAP, PH at 10 WAP, and SD at 4 WAP than HS (Table 2). HS was higher than NHS in LA at 4 and 7 WAP, LCC at 4, 7 and 10 WAP, PH at 4 and 7 WAP, and SD at 7 and 10 WAP. Although HS took a longer number of days to reach DT 50 and DS 50 , its TSI was lower than NHS ( Figure 2). SDWt, CL, CW, CC, CWt, GN, GWt, SP, and SWt were all higher in the NHS.

Varietal Influence on the Progenies' Growth and Yield Attributes
Variety WE5323 produced the largest LA, which was significantly higher than ZM1523 and WE3128 at 7 and 10 WAP (Table 3). At 4 and 7 WAP, its LCC was higher than ZM1523 and WE3128. It also produced the tallest PH at 4 WAP as well as the widest SD at 7 and 10 WAP. ZM1523 had the greatest LCC at 10 WAP, the highest NoL at 4 and 7 WAP, the tallest PH at 7 WAP, and the widest SD at 4 WAP. WE3128 produced the highest NoL and PH at 10 WAP, and the highest EP. Its LCC at 7 WAP, SD at 4 and 10 WAP were significantly lower than WE5323 and ZM1523.

Varietal Influence on the Progenies' Growth and Yield Attributes
Variety WE5323 produced the largest LA, which was significantly higher than ZM1523 and WE3128 at 7 and 10 WAP (Table 3). At 4 and 7 WAP, its LCC was higher than ZM1523 and WE3128. It also produced the tallest PH at 4 WAP as well as the widest SD at 7 and 10 WAP. ZM1523 had the greatest LCC at 10 WAP, the highest NoL at 4 and 7 WAP, the tallest PH at 7 WAP, and the widest SD at 4 WAP. WE3128 produced the highest NoL and PH at 10 WAP, and the highest EP. Its LCC at 7 WAP, SD at 4 and 10 WAP were significantly lower than WE5323 and ZM1523. WE3128 required the most days to achieve DT 50 and DS 50 ( Figure 3). It had higher CL, CW, CWt, GWt, and SWt than WE5323 and ZM1523. WE5323 produced the highest CC and GN, whereas ZM1523 had the highest TSI, SDWt, CN, and SP.
WE3128 required the most days to achieve DT50 and DS50 ( Figure 3). It had higher CL, CW, CWt, GWt, and SWt than WE5323 and ZM1523. WE5323 produced the highest CC and GN, whereas ZM1523 had the highest TSI, SDWt, CN, and SP.

Effect of Parental Soil Amendment on the Progenies' Growth and Yield Attributes
The MF amendment produced the highest LCC at 7 WAP, while the PM amendment gave the highest LCC at 4 WAP, the largest NoL at 4 WAP, and the widest SD at 4 and 10 WAP ( Table 4). The MPM amendment had the greatest EP, LA, NoL, PH, and LCC at 10 WAP and the widest SD at 7 WAP. It was significantly higher than MF and PM in LA at 7 and 10 WAP, and PH at 10 WAP.

Effect of Parental Soil Amendment on the Progenies' Growth and Yield Attributes
The MF amendment produced the highest LCC at 7 WAP, while the PM amendment gave the highest LCC at 4 WAP, the largest NoL at 4 WAP, and the widest SD at 4 and 10 WAP (Table 4). The MPM amendment had the greatest EP, LA, NoL, PH, and LCC at 10 WAP and the widest SD at 7 WAP. It was significantly higher than MF and PM in LA at 7 and 10 WAP, and PH at 10 WAP. The longest number of days to complete DT 50 and DS 50 was observed in PM amendment ( Figure 4). The PM amendment was also significantly higher than the MF in CN. The MF amendment produced the longest TSI, widest CW, highest SDWt, and SWt. It was significantly lower in GN compared with MPM and PM amendments. The MPM amendment produced the highest CL, CC, CWt, CN, GWt, and SP. The longest number of days to complete DT50 and DS50 was observed in PM amendment ( Figure 4). The PM amendment was also significantly higher than the MF in CN. The MF amendment produced the longest TSI, widest CW, highest SDWt, and SWt. It was significantly lower in GN compared with MPM and PM amendments. The MPM amendment produced the highest CL, CC, CWt, CN, GWt, and SP.

Interaction of Parental Heat Stress, Maize Variety and Soil Amendment on the Progenies' Growth and Yield Attributes
EP ranged from 54.2% in HS-WE3-PM-P 1 to 97.9% in HS-ZM1-MPM-P 1 ( Table 5). At 4 and 7 WAP, the LA ranged from 229.1 cm 2 to 354.8 cm 2 , and from 566.7 cm 2 to 743 cm 2 , respectively. The LA ranged from 597.9 cm 2 in HS-ZM1-PM-P 1 to 743 cm 2 in HS-WE3-MPM-P 1  Among the treatment interactions at 4 WAP, the shortest and tallest plants were produced by HS-WE3-PM-P 1 and HS-WE5-MPM-P 1 , respectively (Table 6). PH ranged from 102.2 cm to 158.8 cm at 7 WAP, while at 10 WAP, it ranged from 223.7 cm to 279.3 cm. SD ranged from 11.2 mm to 20.3 mm at 4 WAP, from 27.4 mm to 38 mm at 7 WAP, and from 25.3 mm to 36.7 mm at 10 WAP.
NHS-ZM1-MPM-P 1 produced the fewest number of days to attain DT 50 and DS 50 , whereas the longest number of days to reach DT 50 and DS 50 came from NHS-WE3-PM-P 1 and HS-ZM1-PM-P 1 . The TSI ranged from 1.3 d to 3.7 d. SDWt ranged from 182 g plant −1 in HS-WE3-PM-P 1 to 424 g plant −1 in NHS-WE5-PM-P 1 .
NHS-ZM1-PM-P 1 had the highest CN, GN and SP, while HS-ZM1-MF-P 1 gave the highest SWt (Table 7). NHS-WE3-PM-P 1 had the highest CL, CW, CC, GWt and CWt. The following ranges were obtained for the yield attributes: The first and second principal component axes explained 69.51% (PC1 = 46.45%, PC2 = 23.06%) of the yield attributes variations in the treatment interactions. The treatment interactions were divided into three groups and an outlier using the principal component biplot ( Figure 5). Group I consist of the following: H1, H2, H8, N2, N5, N6, and N7. The outlier N9 represented Group II, while H5 and N3 comprised Group III. Group IV contained the remaining treatment interactions, which included H3, H4, H6, H7, H9, N1, N4, and N8.      The first and second principal component axes explained 69.51% (PC1 = 46.45%, PC2 = 23.06%) of the yield attributes variations in the treatment interactions. The treatment interactions were divided into three groups and an outlier using the principal component biplot ( Figure 5). Group I consist of the following: H1, H2, H8, N2, N5, N6, and N7. The outlier N9 represented Group II, while H5 and N3 comprised Group III. Group IV contained the remaining treatment interactions, which included H3, H4, H6, H7, H9, N1, N4, and N8.  Group I had the highest SDWt group average (Table 8). It performed better than the population mean in CN, CL, CW, CC, CWt, GN, GWt, and SP. It did, however, provide the least SWt. The CN, GN, and SP of the outlier N9 were the highest. In SDWt, CL, CWt, GWt, and SWt, it performed better than the population mean, whereas in CW and CC, it performed below the population mean. Among all the groups, Group III had the best CL, CW, CC, CWt, GWt, and SWt. Its CN and GN were higher than the population mean, but it had the lowest SDWt and SP. Group IV performed above-population mean in SWt and below-population means in SDWt and SP. It produced the least group means in CN, CL, CW, CC, CWt, GN and GWt. The grain weight and cob weight were the most discriminatory and representative yield attributes, followed by CC, CL, and GN, whereas SP and SWt were the least discriminatory and representative attributes. This categorization was based on each attribute' s vector length and angle to the biplot origin ( Figure 5).

Effect of Parental Heat Stress on the Progenies' Growth and Yield Attributes
Heat waves are expected to become more severe and frequent, posing a risk to global food security [17]. Heat stress has a negative effect on maize plants [3,21]. As plants are exposed to constant temperature fluctuations in the field, improving their thermotolerance is critical for their survival [7]. Plants exposed to sublethal heat stress (heat priming) during an early growth phase improved their thermotolerance [1,[5][6][7]12,23]. Wang, et al. [23] reported that the progeny of heat-primed wheat plants produced higher grain yield, better leaf photosynthesis, higher antioxidant enzyme activity, and less cell membrane damage when exposed to high temperatures than unprimed (control) wheat plants. Heat priming in the first generation induced transgenerational thermotolerance, according to the study [23]. Kushawaha et al. [13] observed that heat priming improved the grain yield of a heatsensitive rice variety (Pusa Basmati I) under heat-stress conditions. Due to the retention of the initial stress memory, heat primed plants outperform unprimed plants under heat stress conditions. The duration of the stress, the severity of the stress, and the developmental stage in which the plant encountered the stress can all influence the progeny's ability to remember the stress [5][6][7][8]. Stress memory transfer occurs as a result of epigenetic inheritance, which involves heritable patterns of phenotypic variation that are unrelated to DNA sequence variation [11].
According to Sultan, et al. [24], depending on the progeny's growth environment, stress memory can have an adaptive or maladaptive effect on progeny development. Plants' development and yield are hampered by the expression of maladaptive stress memory [6,24]. The progenies of drought-primed wheat outperformed their control counterparts under drought conditions, but there was no significant difference when the two progenies were grown under control conditions [8]. A similar observation was made for Arabidopsis thaliana, with no significant difference in progenies that differ in ancestral temperature exposure under control conditions, but the descendants of the heat-primed ancestor outperformed the descendants of the unprimed ancestor in heat stress environments [1,9].
The non-significant differences in the heat primed and control progenies under control environments observed in the aforementioned studies [1,9] corresponded to the nonsignificant differences in this study between the HS and NHS progenies under control environments. The lower EP in the HS progenies corroborated Kaushal, et al. [25] findings that heat stress slows seedling emergence. The HS progenies' higher LA, LCC, and PH could be early adaptive traits, as the plants experienced the highest maximum and minimum field temperatures during the first ten WAP (Figure 1). Wang, et al. [8] found that primed wheat progenies had higher photosynthetic activity, antioxidant capacity, and osmolyte accumulation than control progenies under stress.
When the temperature in the field dropped, the NHS progenies produced more LA, NoL, and PH, which likely resulted in more SDWt, CL, CW, CC, CWt, GN, GWt, SP, and SWt. According to Crisp, et al. [6], the events that occur during the stress recovery period determine whether the acquired stress memory is adaptive or maladaptive. A proper recovery and memory reset will result in resources being allocated to growth and/or reproduction, whereas an impaired recovery may result in resources being allocated to the acclimation process [6]. This could explain the HS progenies' longer vegetative phase as measured by DT 50 and DS 50 , as well as their lower yield attributes when compared to the NHS progenies in this study. Plant genotypes play a role in the creation and retention of stress memory, in addition to abiotic and biotic stress effects [16,26].

Varietal Influence on the Progenies' Growth and Yield Attributes
Environmental stresses affect maize plants in different ways. Therefore, high yield consistency in both favourable and unfavourable conditions is a desirable attribute in maize. According to Edreira and Otegui [27], the identification of stress-tolerant maize varieties is a critical step toward reducing maize yield losses. The significant differences in the maize population's growth and yield attributes observed in this study are indicators of their diversity. Balla, et al. [5] found significant differences between wheat cultivars, not only in yielding abilities under control conditions, but also in heat sensitivity. When stress occurs on a regular basis, genotypic adaptation to stress is beneficial [28]. Heat stress is an environmental challenge that plants in natural populations face [10], and it is expected to become more severe and common in the future [17]. Therefore, improving thermotolerance in maize is of strategic importance to global food security.
In comparison to WE5323 and ZM1523, the most depressive heat stress effect was seen in variety WE3128 during parental heat stress priming. According to Crisp, et al. [6], maternally stressed plants may allocate more resources to their seeds in order to increase the survivability of their progenies, in preparation for the stress experienced by the mother plant. The higher values recorded in important attributes such as EP, CWt, and GWt in WE3128 may be an indication of adaptive stress memory, even though there were no significant differences in most of the yield attributes measured in the maize varieties' progenies. The formation and transfer of stress-memories is genotype-dependent [14,26].

Effect of Parental Soil Amendment on the Progenies' Growth and Yield Attributes
Biotic and abiotic stress conditions, such as nutrient deficiency, can cause epigenetic effects. Plants are subjected to nutrient stress when nutrients are in short supply or are available in forms that are inaccessible to the plants. Under such a stressful environment, plants change their physiology to adapt to the new environment. These phenotypic changes have been shown to influence progenies even in second generations that were not exposed to the stress [9,14,15,29].
Kou, et al. [15] observed that N deficiency in rice (Orzya sativa L.) caused DNA methylation in stressed plants, which was passed down to the progeny. Despite the fact that the N stress was only applied to the parents, the stress memory was expressed in the second generations [15].
In this study, progenies of the MF amendment had the lowest levels of CN, GN, CWt, and GWt, all of which were also low in their parents. Similarly, the highest CL, SP, GN, CWt, and GWt obtained from the MPM amendment progenies were consistent with the observations from their parents. A uniform mineral fertiliser was applied during the progeny evaluations. Hence, similar patterns in yield attributes may be linked to the parental growth environment.
Fan, et al. [14] investigated the effect of N deficiency in rice and found reduced growth and yield in stressed parents and their progenies, even when the progenies were grown in N-sufficient conditions. They hypothesised that the progenies inherited their stressed parents' stress memory, which hampered their growth and development in stressfree environments.

Interaction of Parental Heat Stress, Maize Variety and Soil Amendment on the Progenies' Growth and Yield Attributes
Plants often respond physiologically and genomically to abiotic stresses such as heat waves, droughts, and nutrient deficiency [6]. These stress responses can be reversed or passed down to the progenies of stressed plants after stress recovery. However, the exact mechanisms underlying stress memory formation, transfer, and reversal are not clearly understood [6,8,16]. Studies have shown that plants respond differently to single stress encounters and repeated occurrences of the same stress [5,28,30]. Similarly, a plant's response to a single source of stress differs from its response to multiple stress factors [31,32]. Sun, et al. [32] proposed that plant evolution included both resistance to a single type of stress and the ability to counter multiple types of environmental stresses at the same time.
Balla, et al. [5] used morphological and grain yield-related attributes to identify seven distinct groups in the wheat population exposed to single and repeated heat stress. Since grain yield determines the economic value of maize varieties, the treatment interactions' progenies were clustered in this study based on yield-related attributes. The findings in this study that GWt and CWt were the most informative attributes agreed with previous reports on maize studies [3,33,34].
The combination of heat stress and soil amendment elicited a wide range of phenotypic responses in the progenies. The progenies were not grouped based on a single factor when they were clustered. Instead, the population's variation may have been influenced by a possible interactive effect of the factors. Herman and Sultan [10] observed that stress memory may not be expressed uniformly among close relatives. This lack of alignment in the treatment interactions' progenies was consistent with previous studies [1,8,9] that found no significant differences between the progenies of primed and control plants in a stress-free environment. Clusters II and III contained the best performing treatment interaction for each maize variety among the four identified clusters in the biplot. WE3128 and ZM1523 (NHS-WE3-PM [N3] and NHS-ZM1-PM [N9]) produced their highest yields from PM amendments in the NHS environment, whereas WE5323 produced its highest yield from the MPM amendment in the HS environment (HS-WE5-MPM [H5]).
Phenotypic variation within and between populations can be provided by transgenerational plasticity, allowing for further adaptive evolution [35]. As a result, persistent transgenerational effects in controlled environments can serve as a crop diversity reservoir that plant breeders can tap into if a new environmental condition arises [29]. Increased expression of transposable elements results in more mutations in plants when they are stressed by the environment [26].

Conclusions
The parental growth environment influenced the progenies' field performance in a non-heat stress environment. Further evaluation of the progenies under a heat stress environment and molecular analyses are required to establish that a transgenerational effect has occurred. Author Contributions: Conceptualization, U.P.C., F.R.K. and S.M.; methodology, U.P.C., F.R.K. and S.M.; formal analysis, U.P.C.; investigation, U.P.C.; data curation, U.P.C.; writing-original draft preparation, U.P.C.; writing-review and editing, U.P.C., F.R.K. and S.M.; and supervision, F.R.K. and S.M. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.
Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.