Barley yield formation under abiotic stress depends on the interplay between flowering time genes and environmental cues

Since the dawn of agriculture, crop yield has always been impaired through abiotic stresses. In a field trial across five locations worldwide, we tested three abiotic stresses, nitrogen deficiency, drought and salinity, using HEB-YIELD, a selected subset of the wild barley nested association mapping population HEB-25. We show that barley flowering time genes Ppd-H1, Sdw1, Vrn-H1 and Vrn-H3 exert pleiotropic effects on plant development and grain yield. Under field conditions, these effects are strongly influenced by environmental cues like day length and temperature. For example, in Al-Karak, Jordan, the day length-sensitive wild barley allele of Ppd-H1 was associated with an increase of grain yield by up to 30% compared to the insensitive elite barley allele. The observed yield increase is accompanied by pleiotropic effects of Ppd-H1 resulting in shorter life cycle, extended grain filling period and increased grain size. Our study indicates that the adequate timing of plant development is crucial to maximize yield formation under harsh environmental conditions. We provide evidence that wild barley alleles, introgressed into elite barley cultivars, can be utilized to support grain yield formation. The presented knowledge may be transferred to related crop species like wheat and rice securing the rising global food demand for cereals.


Results and Discussion
HEB-YIELD exhibits strong phenotypic variation as well as environmental and treatment variation. The wild barley introgression population HEB-YIELD comprises a diverse subset of lines selected from the NAM population HEB-25 18 (Supplementary Table S1a). We studied eleven agronomically traits in a HEB-YIELD trial conducted in Dundee, Halle, Al-Karak, Dubai and Adelaide ( Fig. 1; Supplementary Table S2a), where climate data for day length, temperature and precipitation varied considerably between locations Supplementary Figs S1 and S2). The parameters studied included developmental and yield-related traits, used to capture growth variation among HEB-YIELD lines (Supplementary Table S3). At each location, the traits were measured under site-specific abiotic stress conditions, i.e. nitrogen deficiency in Dundee and Halle, drought stress in Al-Karak and Adelaide and salt stress in Dubai. In total, 3,207 field plots were evaluated over all sites, seasons, treatments, and replicates (Supplementary Tables S2c and S4a). Considerable phenotypic variation within locations and treatments was observed for all investigated traits ( Fig. 2; Supplementary Table S5).
The ANOVA revealed that all investigated factors (genotype, year and location) were significant for all traits except plant height (HEI) where the year effect was not significant (Supplementary Table S6a). Interestingly, only the traits shoot elongation phase (SEL) and HEI showed comparable values across locations, whereas for the majority of traits pronounced location effects were observed (Supplementary Table S6c). For instance, flowering time varied from 57 to 144 days and grain yield from 0.14 dt/ha to 74 dt/ha (Supplementary Table S5), reflecting a strong diversity in yield potential among the trial sites (Fig. 2).
Irrespective of the diverging agricultural practices at the trial sites, developmental trait heritabilities were high with an average of 0.87, ranging from 0.10 (ripening phase (RIP) under control treatment in Al-Karak) to 0.99 (shooting (SHO) under control treatment in Adelaide as well as flowering (HEA) and SEL under both treatments in Dubai, Supplementary Table S5). In general, yield-related traits revealed lower heritabilities with an average of 0.65. The most complex trait, grain yield (YLD), revealed average heritabilities of 0.73, ranging from 0.05 (YLD under stress treatment in Dubai) to 0.93 (YLD under control treatment in Dundee).

Trait performance in HEB-YIELD is usually a linear transformation from control to stress treatments indicating low genotype by treatment interaction.
To gain insights into how abiotic stresses may affect plant development and grain yield, we cultivated HEB-YIELD under contrasting stress conditions, which are relevant for the respective test locations (Supplementary Table S2c). The applied stresses exhibited only minor effects on plant development traits except for HEI. In contrast, strong effects on all measured yield-related traits were observed at all test locations, for instance, reducing yield under stress between 16% in Halle and 65% in Adelaide ( Fig. 2; Supplementary Table S5). We observed a weak trend, that HEB-YIELD lines under drought and salt stress exhibited an accelerated plant development, presumably to escape the stress condition, which is in agreement with other studies in cereals 32,41 . Based on our findings we suggest that plant development in the wild barley population HEB-YIELD is mainly determined by genetic factors and to a lesser extend modified by abiotic stresses. This is further supported by the observation that plant developmental traits showed a nearly linear shift between control and stress conditions, as indicated by high correlation coefficients (0.99 > r > 0.59) between stress and control treatments of developmental traits, except for SEL in Adelaide (r = 0.12; Table 1; Supplementary  Table S7a).
Grain yield correlations indicate that yield formation depends on a location-specific interplay between developmental traits and yield components. We observed Pearson correlations coefficients between plant developmental stages shooting, flowering and maturity ranging from r = 0.67 to r = 0.96 (apart www.nature.com/scientificreports www.nature.com/scientificreports/ from shooting correlations in Dubai Supplementary Tables S6b and S7b), indicating a nearly colinear regulation of plant developmental phases. Thus, HEB-YIELD lines early or late in shooting have the tendency to stay early or late respectively until maturity. This observation is in agreement with previous findings in the wild barley NAM population HEB-25, studied in Halle 20 and Dundee 25 . Consequently, early developmental stages may be used as an indirect criterion to select HEB-YIELD lines for early or late maturity.
Following these findings, we explored the relationship between plant development and yield formation in HEB-YIELD (Table 2). We observed a trend that late plant development is beneficial for increased grain yield under Dundee, Halle and Adelaide growth conditions, indicated by positive correlation coefficients of r(HE-AxYLD) = 0.59/0.66, 0.32/0.20 and 0.57/0.51, respectively, under control/stress treatments. This trend fits the general observation that late lines have the potential to exploit a prolonged growing season if the environmental conditions including temperature and precipitation are beneficial 29,42,43 . In contrast, under the harsh environmental conditions at Al-Karak and Dubai, HEB-YIELD lines with accelerated plant development were favored. Consequently, we observed negative correlations between flowering and grain yield at Al-Karak and Dubai with r(HEAxYLD) = −0.30/−0.72 and −0.51/−0.44, respectively, under control/stress treatments (Table 2). Here, elevated temperatures and low rainfall restricted plant growth to a few months and thus earliness is a major breeding goal 44,45 . In future, this situation may intensify, since climate change is expected to further shorten the growing period in drought and heat prone locations like in Jordan 26,32,46 .
We also observed strong location-specific correlations between flowering time and yield components. For example, in Halle and Dundee, flowering time was positively correlated with grain number per ear (GNE), with r(HEAxGNE) = 0.67/0.67 and 0.71/0.64, respectively (Table 2). Here the extended vegetative growth phase allowed more spikelet primordia to be maintained. In contrast, in Al-Karak flowering time negatively affected thousand grain weight (TGW) with r(HEAxTGW) = −0.56/−0.68, reflecting a grain filling penalty for later www.nature.com/scientificreports www.nature.com/scientificreports/ flowering genotypes. These findings suggest that flowering time controls final grain yield to a certain degree. By comparing correlations between grain yield and yield components, we observed that apparently GNE is the key determinant of grain yield in HEB-YIELD -irrespective of location and treatment. This result is in agreement with earlier studies 47,48 , indicating that any increase in number of grains may also improve grain yield 49,50 . We thus reason that improving GNE may offer the best route to increase grain yield in HEB-YIELD independent of the environmental conditions.
The highest positive correlations of yield were found with harvest index (HI; scored only in Dubai and Adelaide with r(YLDxHI) = 0.87 and 0.83, respectively, Table 2). A previous study noted the importance of increasing harvest index to improve yield during the past century 51 . However, a further improvement of grain yield through raising harvest index may be a dead end, since barley is supposed to have reached an optimum with a harvest index of approximately 0.62 51,52 . Therefore, future grain yield improvements may be achieved through increasing plant biomass 51,53 . This suggestion is in accordance with our finding that grain yield exhibited a slightly positive correlation with shoot elongation phase in those environments where lateness was beneficial to increase yield (Table 2). During shoot elongation, which captures the growth period between establishing awn primordia and ear emergence, the leaf growth rate and the potential grain number per area are defined 47,49,54 . An extended shoot elongation phase may thus improve grain yield by increasing leaf size, i.e. biomass, and grain number per area. On the other hand, ripening phase under drought stress exhibited positive and negative correlations with grain yield in Al-Karak and in Adelaide, respectively. Whereas the Adelaide finding fits the assumption that early Adelaide (ADE) are indicated with blue, grey, green, red and yellow box-whiskers, respectively, and, in addition, at the bottom of the plot. Empty and filled boxes refer to control and stress treatments, respectively. Significant differences between treatments are indicated with red asterisks above boxes with *P < 0.05, **P < 0.01 and ***P < 0.001. The relative increase/decrease (in %) of the stress treatment compared to the control treatment is given below the asterisks.
www.nature.com/scientificreports www.nature.com/scientificreports/ maturity and thus a short ripening phase may improve grain yield under terminal drought, the Al-Karak finding is unexpected. Under drought stress conditions in Al-Karak, an extended ripening phase was associated with an increase in grain weight, ultimately resulting in elevated grain yields. We conclude that fine-tuning of plant development, especially their sub-phases, may contribute to a better adaptation of improved varieties to their target environment. The latter notion is supported by the finding that in the first instance climate change is expected to impair flowering time 32 , which is crucial for plant adaptation and yield formation 49,55 . In addition, our stress treatments confirmed the known association between inflorescence development and stress tolerance/avoidance 32,56 . This offers the possibility to use the genetically relatively well-understood trait flowering time as a proxy to select for improved grain yield under abiotic stresses 57 . Flowering time genes exhibit pleiotropic effects on yield formation in HEB-YIELD. In order to explore the interplay between flowering time regulation and yield formation, we investigated the effects of four major flowering time genes, Ppd-H1, Sdw1, Vrn-H1 and Vrn-H3, in HEB-YIELD (Supplementary Table S8a). The relevance of these candidate genes has been reported in various studies 34 www.nature.com/scientificreports www.nature.com/scientificreports/ cultivated alleles at these four flowering time loci (Supplementary Table S1a). In the following, we report on the pleiotropic effects associated with the four flowering genes studied.

Ppd-H1.
Flowering under long days is promoted by the photoperiod responsive Ppd-H1 (PHOTOPERIOD-H1) allele, an orthologue of the Arabidopsis pseudoresponse regulator gene PRR7, which is present in wild barley and winter barley cultivars 37 . In contrast, spring barley cultivars like Barke possess the recessive non-responsive ppd-H1 allele, resulting in late flowering. Allelic variation in Ppd-H1 follows a latitudinal cline, where a large proportion of the dominant, long day responsive allele is present in cultivars from Southern Europe and the recessive, non-responsive allele, predominates in Northern Europe 58 . Selection of the non-responsive alleles found in elite spring barley allowed barley cultivation in northern latitudes, under more humid and colder conditions 59 . During early plant development, the photoperiod signal is transmitted from the circadian clock oscillator Ppd-H1 through mediation of the CONSTANS (CO) protein to the floral inducer Vrn-H3, an orthologue of the Arabidopsis FLOWERING LOCUS T (FT) gene.
In order to locate sequence variants in Ppd-H1 discriminating between wild and cultivated barley in the set of HEB-YIELD parents, we used exome capture sequence data (Supplementary Table S1b 37 . SNP 22 separates winter and wild barley from spring barley and is located in the CCT domain of Ppd-H1, where the G to T substitution causes a glycine to tryptophan change resulting in a photoperiod non-responsive mutant phenotype 37 . In addition, 13 Ppd-H1 variants (15, 27, 30, 37, 39, 42, 51, 62, 64, 88, 99, 112, 116 46) discriminated between the tested wild barleys and all tested spring barleys except HEB-25 donor HID138, originating from Iran but not present in HEB-YIELD 18 Supplementary Table S3. The color of the bars represents the location, blue for Dundee, grey for Halle, green for Al-Karak, red for Dubai and yellow for Adelaide. Ppd-H1 wild allele effects under control and stress treatments are depicted with a bright blue (top) and a bright red background (bottom), respectively. Statistically significant wild allele effects are indicated by red asterisks above or below the bars with *P < 0.05, **P < 0.01, ***P < 0.001. The height of the bars indicates the size of the Ppd-H1 wild allele effect, obtained by calculating the difference between the mean performances of HEB-YIELD lines carrying two wild alleles versus two elite alleles.
www.nature.com/scientificreports www.nature.com/scientificreports/ significant effect of the Ppd-H1 wild allele was observed in Dubai and Adelaide. Most wild barley accessions carry the dominant allele, which is responsive to a long day photoperiod, accelerating plant development through upregulating of Vrn-H3/HvFT1 62,63 . One possible explanation for contrasting effects between locations is the different day lengths at these sites. Dundee and Halle are more than 5,700 km distant from the equator and are clearly exposed to long day conditions indicated by average day lengths of more than 15 hours during shooting phase, which is necessary to trigger the effect of Ppd-H1 37 . In Al-Karak, with a day length of appr. 12 hours during shooting phase, we still observed strong Ppd-H1 effects, although this location is more than 2,000 km closer to the equator than Halle. In Dubai, where day length is shorter with less than 11 hours during shooting phase, only one minor Ppd-H1 effect on plant development and yield formation was observed (Supplementary Table S8b). Apparently, this is because Ppd-H1 is only active under long-day condition 37 . Presumably, the short-day signal to initiate flower development in Dubai is transmitted through Ppd-H2, the FLOWERING LOCUS T3 (FT3) gene 64 . The dominant functional allele of Ppd-H2 promotes spikelet formation under short-day in spring barley and winter barley originating from Southern Europe whereas the recessive non-functional allele carries a large deletion in the transcribed coding region, which is typically present in Northern European winter barley cultivars [65][66][67] . Based on exome capture sequence data, we found nine sequence variants in Ppd-H2 (Supplementary Table S1b). No sequence variant could be identified, which perfectly discriminated between cultivated and wild barley. However, Ppd-H2 variant 3 at chr1H_514.098.364 bp and variant 7 at chr1H_514.098.702 bp discriminated between the tested winter barleys and the tested spring barleys except Morex. The wild barley HEB-YIELD donors showed SNPs in common but not restricted to spring or winter alleles.
In future, follow-up field studies with double HIFs (heterogeneous inbred families) may assist to further characterize the interplay between the Ppd-H1 and Ppd-H2 photoperiod receptors under short-day and long-day conditions 68 . For this, four HEB-25 lines, HEB_05_044, HEB_08_149, HEB_16_063 or HEB_22_039, which were simultaneously heterozygous at both loci in generation BC1S3, can be chosen based on Maurer et al. 18 . Subsequently, four homozygous allele combinations at the two loci can be selected in the available HEB-25 selfing generation BC1S3:11. After seed multiplication, epistatic effects on flowering time and yield formation can be tested in replicated field trials using the resulting four nearly isogenic double HIFs.
In HEB-YIELD, Ppd-H1 acted in a location-specific manner on yield-related traits. The most pronounced Ppd-H1 effect was present in Al-Karak where the day length-sensitive wild barley allele was associated with an increase of grain yield by 4.1 dt/ha (+15%) and 5.3 dt/ha (+30%) under control and drought stress conditions, respectively ( Fig. 3; Supplementary Table S8a). The yield effect may be explained through pleiotropic effects of the wild barley Ppd-H1 allele, which shortened the overall growing season, increased the period of grain filling (RIP) and increased grain size (TGW). A tendency of the Ppd-H1 wild barley allele towards enhanced grain yields was also observed in Dubai and Adelaide, however, only significant in Adelaide under drought stress (+4.8 dt/ha = +29%). Usually, the location-specific effects of Ppd-H1 on yield-related traits are in agreement with the preferred length of the growing period. At those locations where earliness is beneficial, the responsive wild allele of Ppd-H1 exerted increasing effects on yield-related traits, for example in Al-Karak, where early plants escaped higher temperatures and terminal drought at the end of the growing season. On the other hand, where lateness is preferable to achieve higher yields, the elite barley ppd-H1 allele increased yield-related traits, for example, in Dundee and Halle. At those locations late HEB-YIELD lines benefited from the extended growing period since the environmental conditions supported plant growth under suitable conditions. Sdw1. Sdw1 belongs to the group of so-called semi-dwarfing genes 40 , which are responsible for yield elevations during the 'Green Revolution' 69 . Wild barley accessions possess the functional and dominant Sdw1 allele, a gibberellic acid 20 oxidase (GA20ox) gene, which promotes plant growth. In contrast, the recessive, GA-deficient sdw1 allele 70,71 is present in barley cultivars like Barke, causing a semi-dwarf phenotype. Several studies have shown that semi-dwarfs exhibit reduced plant height, late maturity, increased tiller numbers and an improved harvest index, ultimately resulting in elevated grain yields 40,72,73 . Based on exome capture sequence data, we found 46 sequence variants in Sdw1 (Supplementary Table S1b). No Sdw1 sequence variant could be identified to completely discriminate cultivated and wild barley. However, variant 130 at chr3H_634.078.282 bp discriminated all tested cultivated spring and winter barleys from the tested wild barleys except of the three HEB-25 donors HID003, HID099 and HID114, originating from Iraq, Syria and Lebanon, respectively 18 . It is, thus, possible that new GA20ox mutations are present in wild barley.
The reported pleiotropic effects of Sdw1 are also supported by HEB-YIELD field data ( Fig. 4; Supplementary Table S8a). Throughout plant development, we detected an accelerating effect of the wild barley Sdw1 allele in HEB-YIELD, accelerating grain maturity by 4.0 to 8.9 days in Dundee, Al-Karak, Dubai and Adelaide, compared to the semi-dwarfing allele of Barke. Most striking was the pronounced delay of development in Adelaide under the control condition (precipitation = 484 mm), with up to 13 days for SHO. Whereas under stress (precipitation = 159 mm) the effects were on a similar level as in the other locations. The Adelaide effect might be explained by different environmental cues between the two years, resulting from the earlier sowing date and the prolonged growing period of 50 days in 2016. So far, there is no evidence that day length or precipitation affects the function of Sdw1 40,71 . However, a wheat survey under controlled conditions already reported that temperature can modify GA dependent responses, where elevated temperatures increase the abundance of GA 74 .
The most prominent effect of semi-dwarfing genes is their control of plant architecture, in particular, plant height 40,75 . We confirmed this effect in HEB-YIELD since the wild barley allele increased plant height at all locations and under both treatments with a maximum increase of 33.2 cm in Dundee under control condition. The dominance of semi-dwarf genes in modern crop cultivars indicates their global importance for agriculture 69,76 .
This notion is also confirmed in HEB-YIELD where the Barke semi-dwarf allele was associated with an increase in grain yield ( Fig. 4; Supplementary Table S8a). In turn, the wild barley allele significantly reduced grain yield, for instance under control conditions in Dundee, Halle and Adelaide by up to 15.8 dt/ha. Under drought www.nature.com/scientificreports www.nature.com/scientificreports/ stress conditions in Adelaide, the Barke semi-dwarf allele revealed the strongest impact, accounting for 60% of the final yield level. The observed yield increase may be attributed to an accumulation of several positive effects, including an extended growing period, more tillers, a higher harvest index and less lodging and head loss ( Fig. 4; Supplementary Table S8a).

Vrn-H1.
In addition to the aforementioned photoperiod and GA dependent pathways, flowering time is also regulated through the vernalization pathway, where exposure to cold temperatures accelerates flowering 33,77,78 . In barley, the response to cold temperatures is mainly controlled by interaction of the two vernalization genes, Vrn-H1 79  Based on exome capture sequence data, we found a huge number of 377 sequence variants in Vrn-H1 (Supplementary Table S1b). A number of Vrn-H1 sequence variants discriminated between cultivated and wild barley. For example, variants 475 and 499 at chr5H_599.131.041 bp and chr5H_599.131.479 bp, respectively, discriminated the tested spring barleys from the tested wild and winter barleys. Interestingly, we found one variant, 443 at chr5H_599.130.360 bp, which discriminated the tested spring and winter barleys from the tested wild barleys. In addition, ten Vrn-H1 variants (187, 194, 289, 296, 335, 340, 388, 452, 457, 498), discriminated the tested winter barleys from the tested spring and wild barleys. The named variants are exclusively located in introns between exon 1 and exon 4. It was already known that deletions in the first intron of Vrn-H1 result in spring type cultivars, lacking the vernalization need to initiate flowering 82 . Our findings indicate that additional intron regulatory elements may be present in Vrn-H1 to differentiate winter barley, spring barley and wild barley.
In future, follow-up field studies using nearly isogenic HIFs may assist to characterize developmental and yield formation effects of individual Vrn-H1 variants, which are present in HEB-YIELD lines 68 . Likewise, the epistatic interaction between selected alleles of Vrn-H1 and Vrn-H2 may be characterized in double HIFs, which can be developed from any of five double heterozygous HEB-25 lines, HEB_09_101, HEB_16_095, HEB_16_099, HEB_23_061 or HEB_24_066, as mentioned before.  Supplementary Table S3. The color of the bars represents the location, blue for Dundee, grey for Halle, green for Al-Karak, red for Dubai and yellow for Adelaide. Sdw1 wild allele effects under control and stress treatments are depicted with a bright blue (top) and a bright red background (bottom), respectively. Statistically significant wild allele effects are indicated by red asterisks above or below the bars with *P < 0.05, **P < 0.01 or ***P < 0.001. The height of the bars indicates the size of the Sdw1 wild allele effect, obtained by calculating the difference between the mean performances of HEB-YIELD lines carrying two wild alleles versus two elite alleles.
www.nature.com/scientificreports www.nature.com/scientificreports/ In HEB-YIELD, Vrn-H1 exhibited considerable effects on nearly every trait in Dubai and, to a lesser extent, in Adelaide ( Supplementary Fig. S3; Supplementary Table S8a). HEB-YIELD lines carrying the wild barley allele at this locus delayed flowering time and maturity by more than 10 days in Dubai. In Adelaide, pronounced effects on plant development were restricted to the control condition (i.e. the Adelaide growing period 2016). Most likely, this effect is caused by warmer temperatures and therefore less vernalization stimuli at the beginning of the growing season (Supplementary Figs S1 and S2; Supplementary Table S9). However, the late development effect of the wild barley Vrn-H1 allele in Adelaide diminished during cultivation from +12 days at shooting, +6 days at flowering to, finally, +4 days at maturity. This tendency was also present in Halle and Al-Karak, although on a much lower level. In contrast, the late development effect of the wild barley Vrn-H1 allele remained stable throughout plant cultivation in Dubai. This may be because the temperature in Dubai never reached a vernalization-triggering level. In Dubai, the HEB-YIELD lines possessing a wild barley winter allele at Vrn-H1 thus responded to the lack of vernalization with a late plant development.
In addition to its developmental effects the wild barley allele of Vrn-H1 exerted significant reducing effects on all yield components of around 25% in Dubai. Consequently, the final grain yield in Dubai was reduced by 3.2 dt/ha under control conditions, which corresponds to 37% of the total yield. At locations where earliness is the preferred breeding goal and vernalizing conditions are rare, the use of the dominant elite barley allele of Vrn-H1 is highly recommended.

Vrn-H3.
As mentioned before, the expression level of Vrn-H1 increases with exposure to cold temperatures, resulting in flower induction through repression of Vrn-H2 and activation of Vrn-H3 29,82 . Vrn-H3 corresponds to the HvFT1 gene, which is an orthologue of the Arabidopsis FT gene, the so called 'florigen' 33,[83][84][85] Vrn-H3 plays a central role in flower induction integrating photoperiod and vernalization signals 84 . Barley alleles of Vrn-H3 vary regarding the first intron sequence, promoter sequence and copy number. They are widely distributed over winter and spring growth habits [84][85][86] . Unfortunately, we could not identify Vrn-H3 variants in our exome capture sequence data. Presumably, Vrn-H3 produced no variants since all reads were identified as multi-mappers, located at two or more genomic regions simultaneously, and, hence, were ignored by the variant caller. (M. Bayer, personal comm.).
HEB-YIELD field data validated the role of Vrn-H3 on plant development throughout the growing period in all locations except from Dubai ( Supplementary Fig. S4; Supplementary Table S8a). The wild barley allele of Vrn-H3 slowed down plant development between 2.2 and 6.6 days. Generally, winter genotypes are characterized by carrying a recessive Vrn-H3 allele, which displays a reduced expression 85 . Most wild barleys possess a winter type 62 and probably harbor a recessive vrn-H3 allele, which explains the decelerating developmental effects.
Although Vrn-H3 plays an important role for plant development, we identified only weak, mostly non-significant, impacts on yield-related traits. Only in Al-Karak, the wild allele showed significant reducing effects on grain number per ears (under both treatments) and on grain yield under drought stress (−3.3 dt/ha).

The best wild barley HEB-YIELD lines match the yield performance of high-yielding local check
cultivars. The usefulness of wild accessions, related to crop species has been proposed and demonstrated frequently 14,15,87 . Wild barley accessions, in particular H. v. ssp. spontaneum, the progenitor of cultivated barley have been used to improve disease resistance 24,88 and abiotic stress tolerance 22,73,89,90 , as well as plant developmental traits 18,20,25 and quality traits 88,91,92 . The successful use of wild relatives to increase grain yield of barley has not been reported frequently, some exceptions are available 73,93,94 . This may be because of the negative impacts of linked deleterious wild alleles, a phenomenon generally referred to as 'linkage drag' 95 . The HEB-YIELD lines offer the possibility to estimate potentially positive wild allele effects in an adapted genetic background, since they are embedded through backcrossing into the modern elite barley cultivar Barke. In addition, the elite genetic background enables the direct use of HEB-YIELD lines in barley breeding programs.
Based on our two-year field trials, we identified five high yielding HEB-YIELD lines, which showed acceptable grain yield performance, comparable to the recurrent elite parent Barke, across the tested locations. These HEB-YIELD lines are 01_132, 01_104, 10_184, 10_173 and 05_043 (Supplementary Fig. S5; Supplementary Tables S10a-c). They possessed higher grain yields than Barke in Al-Karak (except HEB_10_184). In addition, HEB-YIELD lines 01_132 and 10_184 surpassed the Barke grain yield in Dundee under both stress and control treatments. Furthermore, we identified HEB-YIELD lines, which reached or surpassed the yield level of locally adapted check cultivars (Supplementary Fig. S5; Supplementary Tables S10a-c). These HEB lines are 10_184 and 01_132 in Dundee ( Supplementary Fig. S6), 01_132 and 01_104 in Halle ( Supplementary Fig. S7), 05_043 and 10_173 in Al-Karak (Fig. 5), 15_082 and 06_116 in Dubai (Supplementary Fig. S8) and 10_184 & 01_132 in Adelaide ( Supplementary Fig. S9). For instance, HEB_01_132 surpassed the grain yield of the established local check cultivar 'Navigator' under stress treatment in Adelaide. In addition, under both treatments it was comparable to 'Compass' and 'La Trobe' , which have become the dominant commercial cultivars in South Australia. HEB_01_132 also surpassed the grain yield of the local check '58/1 A' under control treatment in Dubai, indicating that this line may be directly suited for cultivation in the respective environments. Likewise, HEB_05_043 and HEB_10_173 outperformed the check cultivar 'Rum' in Al-Karak under drought stress.
Interestingly, HEB-YIELD lines adapted different yield formation strategies at each test location. Compared to the local check cultivars, HEB-YIELD lines had increased numbers of ears (EAR) at Al-Karak, increased thousand grain weights at Dubai in almost all cases, and increased grain numbers per ear under stress at Dundee and Adelaide in many cases (Supplementary Figs S10-S12). This offers the possibility of achieving future yield improvements following a location-specific adaptation route.
The challenges of climate change demand that cultivars need to re-adapt to changing environmental conditions, for instance shorter growing seasons, higher average temperatures during cultivation and more frequently occurring drought periods 26 Table S5), which offers the potential to use this variation to adapt new cultivars to changing environmental conditions by backcrossing favorable HEB-YIELD donor lines with locally adapted elite cultivars. In those areas where drought and heat affect plant development and grain maturation, early maturing lines like HEB_05_043, HEB_15_082 or HEB_10_173 may be beneficial because of their fast development ( Fig. 5; Supplementary Fig. S9; Supplementary  Table S10a). Moreover, the increased tiller capacity of HEB-YIELD lines may be promising to achieve an improved canopy cover, reducing moisture losses 26,45 and to increase biomass yield. The latter trait may be high value in the eastern part of the Mediterranean basin where straw and grains of barley are mainly used for animal feeding 97,98 .
In the past, Australian varieties followed a strong focus on earliness but changes in agricultural practices have resulted in earlier sowing dates and thus an extended growing season. The earlier sowing allowed later genotypes to benefit from a longer growing period, enabling HEB-YIELD lines HEB_03_085, HEB_10_184, 20_064 and HEB_04_135 to surpass the grain yield of the local check cultivar 'Navigator' in Adelaide under control condition (Supplementary Fig. S9; Supplementary Table S10a). We also identified HEB-YIELD lines that performed quite well in the high yielding environments of Dundee and Halle. Here, lines HEB_01_104, HEB_01_132 and HEB_10_184 accomplished reasonable yields. HEB_10_184, for instance, achieved a maximum grain yield of 74.0 dt/ha under control condition in Dundee, which was almost on par with the local check cultivar 'Odyssey' (−0.3%) and 5.1% higher than the recipient cultivar Barke. Our findings indicate that wild barley HEB-YIELD lines can be used as pre-breeding material to further improve plant development and yield formation of elite barley. (Supplementary Figs S1 and S2; Supplementary Tables S9a and S10a).

Conclusion
It is expected that the impact of climate change necessitates the adaptation of our established crop cultivation systems to harsher environmental conditions 26,96 . Stress avoidance is one promising approach to increase stress tolerance. We explored this relationship by studying the wild barley-derived model population HEB-YIELD in a field experiment, ranging from Dundee in Scotland to Adelaide in South Australia, where the effects of nitrogen deficiency, drought and salinity on plant development and yield-related traits were investigated.
Our findings confirm the crucial relationship between flowering time, plant development and grain yield 99 . The exact timing of the switch from vegetative to reproductive growth under favorable conditions 32 , the length of the growing period and the duration of the sub-phases of plant development are crucial to secure yield under abiotic stress conditions. We suggest that adjusting plant development may be a promising breeding strategy to cope with abiotic stresses. To optimize breeding programs, it is thus advisable to first predict the environment-dependent impact of flowering time genes on yield formation and then to select locally advantageous alleles for sustainable crop improvement.
Our HEB-YIELD data indicate that wild germplasm may serve as a resource to increase genetic diversity 14,20,22 and to enable the above mentioned adaptation to abiotic stresses, through selection of early or late development alleles of known major flowering time genes, e.g. Ppd-H1, Sdw1, Vrn-H1 and Vrn-H3. We showed that allelic variants of these flowering time genes strongly react to environmental cues. This information can be used to www.nature.com/scientificreports www.nature.com/scientificreports/ design novel breeding strategies such as precise backcrossing of suitable developmental genes into regionally adapted cultivars. Our data also provide evidence that wild barley germplasm may be useful to improve yield in low-yielding environments, for instance, in the Middle East, as well as in high-yielding environments, for instance, in Northern and Central Europe. This knowledge may be transferred to related crop species like wheat and rice to secure the rising global food demand for cereals. The HEB-YIELD subset consists of 48 HEB-25 lines that were selected from HEB-25 to ensure the absence of brittleness and a good threshability enabling accurate yield estimation in field trials. In addition, the final HEB-YIELD lines were selected to independently segregate at four major flowering time loci, which exhibited major plant developmental effects in HEB-25: Ppd-H1, Sdw1, Vrn-H1 and Vrn-H3 18,20,25 . Genotypic data. The complete HEB-25 population was genotyped in generation BC 1 S 3 using the barley Infinium iSelect 9k SNP chip (see) 18 . The diagnostic markers i_BK_16, i_12_30924, i_11_10705 and i_12_10218, co-segregating with the four flowering time genes Ppd-H1, Sdw1, Vrn-H1 and Vrn-H3, respectively, were used for selection of HEB-YIELD lines carrying homozygous elite versus homozygous wild barley alleles. In total, HEB-YIELD includes wild barley chromosomal segments derived from 20 wild barley donors. (Supplementary  Table S1a). These donors contributed wild barley alleles at each of the studied four flowering time genes. Supplementary Table S1b provides exome capture-based sequence data of four flowering time genes, Ppd-H1, Ppd-H2, sdw1 and Vrn-H1, collected through the WHEALBI consortium (https://www.whealbi.eu/) and kindly provided by Drs. Micha Bayer and Joanne Russell, The James Hutton Institute, Dundee, UK. At each location, a control treatment and a site-specific stress treatment was applied. Stress treatments were nitrogen deficiency in Dundee and Halle, drought stress in Al-Karak, salt stress in Dubai and drought stress in Adelaide (see Supplementary Table S2c). Therefore, lines of the stress treatment received no nitrogen fertilizer in Dundee and Halle, no drip irrigation in Al-Karak and a saline water drip irrigation in Dubai. In Adelaide, only one treatment was applied per season due to lack of seeds. In this case, the two contrasting seasons represented the treatments where 2015 was regarded as the drought stress treatment with only 159 mm precipitation during the growing period and 2016 as the control treatment with 484 mm precipitation.

Materials and
On average, each HEB-YIELD line was replicated three to four times per treatment. A randomized complete block design was chosen as test design for the trials, with the exception of Dubai and Adelaide where a completely randomized design within each treatment was applied. The trials were conducted in accordance to local practices regarding tillage, fertilization and pest management. Additional information on plant cultivation is provided in Supplementary Table S2b. phenotypic data. Eleven developmental and yield related traits were investigated. A description of where and how each trait was measured is given in Supplementary Table S3. statistical analyses. All statistical analyses were carried out with SAS 9.4 (SAS Institute Inc., Cary, NC, USA) 100 . Variance components (defined as random) were estimated with PROC VARCOMP and broad sense heritabilities (h 2 ) for each trait within locations and treatments were calculated across years following the formula: The analysis of variance (ANOVA) across locations was calculated with PROC MIXED to test for the presence of genotype, location and year effects. For this purpose, the main effects (genotype, location and year), as well as their corresponding interaction effects were treated as fixed effects in the following model: Fixed effects are written in bold Best linear unbiased estimators (BLUEs) were estimated using the PROC MIXED procedure. The BLUEs for each HEB-YIELD line were computed across years and for each treatment level and location separately. Genotype and treatment were modelled as fixed effects and year as a random effect: Fixed effects are written in bold For location Adelaide BLUEs were calculated within years and the model was restricted to a fixed effect of genotype.
Pearson correlation coefficients (r) between trait BLUEs were calculated via PROC CORR. Furthermore, to test for significant treatment effects a simple t-test (PROC TTEST) and an ANOVA within locations were performed (PROC MIXED). The ANOVA model included the main effects (genotype, treatment and year) and their corresponding interaction effects as fixed effects (comparable to model III). In addition, a one-factorial ANOVA was computed to test for significant location effects within treatments where only the main effect (location) was included, followed by a Tukey test (PROC GLM).
Performance of the HEB-YIELD lines was compared to an adapted check cultivar from the corresponding location (see field trials above) by conducting a Dunnett test 101 (PROC MIXED). To enable an easier comparison between the lines the relative performance (RP) was calculated as: To check for significance and estimate effects of the four flowering candidate genes Ppd-H1, Sdw1, Vrn-H1 and Vrn-H3, a simple linear regression model (PROC GLM) was fitted for each candidate gene applying BLUEs across years. Each model included a single locus-specific SNP mentioned above, modeled as a quantitative variable representing the wild allele dosage 20 .