Skip to main content

Protein content and HvNAM alleles in Nordic barley (Hordeum vulgare) during a century of breeding

This article has been updated

Abstract

Background

Barley has been bred for more than a century in the Nordic countries, with dramatic improvements of yield traits. In this study we investigate if this has come at the cost of lower grain protein and micronutrient (iron, zinc) content, by analysing 80 accessions representing four different improvement stages. We further re-sequenced the two grain protein content associated genes HvNAM-1 and HvNAM-2 in full and performed expression analyses of the same genes to search for genetic associations with nutrient content.

Results

We found higher thousand grain weight in barley landraces and in accessions from the late improvement group compared to accessions from the mid of the twentieth century. Straw length was much reduced in late stage accessions. No significant temporal decrease in grain protein, iron or zinc content during twentieth century Nordic crop improvement could be detected. Out of the 80 accessions only two deviant HvNAM-1 sequences were found, represented by one accession each. These do not appear to be correlated to grain protein content. The sequence of HvNAM-2 was invariable in all accessions and no correlations between expression levels of HvNAM-1 and HvNAM-2 and with grain protein content was found.

Conclusions

In contrast to studies in wheat, where a strong negative correlation between straw length and grain protein and micronutrient content has been found, we do not see this relationship in Nordic barley. The last 60 years of breeding has reduced straw length but, contrary to expectations, not protein and micronutrient content. Variation in grain protein and micronutrient content was found among the Nordic barley accessions, but it is not explained by variation of HvNAM genes. This means that HvNAM is an unexploited source of genetic variation for nutrient content in Nordic barley.

Introduction

Barley, Hordeum vulgare, is a major cereal crop with a world-wide production of 141 million tonnes in 2018 (FAOSTAT, http://www.fao.org/faostat/en/). Although some barley grain is used for human consumption its main uses are as animal feed and for malting, with the use for malting and brewing being the most economically important [1]. Barley intended for different uses require different levels of protein in the grain, and, although other protein sources normally complement barley when used as the main carbohydrate source, a high protein content would result in high food and feed quality. In contrast, in malting and brewing too high protein content results in a slow and uneven water uptake during the germination process which lowers the amount of malt extract [2]. In addition, a high protein content in malting barley results in cloudiness in the beer and a reduced shelf life, although sufficient grain protein is needed for the yeast to grow during fermentation and for the production of preferred foam qualities. A suitably low or moderate protein content is thus required from quality malting barley [3].

Grain protein content (GPC) in barley is influenced both by environmental and genetic factors [4, 5]. Mapping studies of GPC points to polygenic control of the trait [3, 6]. In particular, two homologs of the well-studied wheat NAM-B1, HvNAM-1 and -2, have been shown to be associated with GPC [6,7,8,9,10]. In wheat, the NAM-B1 (Gpc-B1) gene has been shown to influence GPC. The gene codes for a NAC transcription factor which has been shown to have pleiotropic effects [11]. The wild-type allele facilitates the remobilisation of nutrients from flag leaves into the maturing grain and increases protein and mineral content while null mutations delay senescence and often result in higher grain yield [12,13,14,15,16,17,18].

In barley, the genetic diversity and physiological effects of the two NAM-genes are less known. Within HvNAM-1, Distelfeld et al. [10] reported two single nucleotide polymorphisms (SNPs) associated with amino acid differences between the high- and low-protein varieties ‘Karl’ and ‘Lewis’, respectively. The associated deposited sequences (EU368851 and EU368852), however, differ at three positions. Fritsch et al. [19] compared 'Karl' with 'Clipper' and reported three SNPs in HvNAM-1. Although no significant association between polymorphism in HvNAM-1 and GPC could be detected, Cai, et al. (2013) showed that a single SNP within HvNAM-2 was associated with GPC. In the same study significant correlations between HvNAM-1, HvNAM-2 and GPC were found among a world-wide set of accessions [6]. NAC genes constitute a large, plant-specific family of transcription factors, many of which have been identified as regulators of senescence in monocotyledons [20]. Similarly to wheat, pleiotropic effects of the HvNAM genes on micronutrient content, senescence and yield have been suggested in barley [21, 22]. In gene expression studies, Christiansen and Gregersen [23] found both HvNAM-1 and HvNAM-2 to be upregulated during senescence.

In the Nordic countries (Denmark, Sweden, Norway, Finland and Iceland) barley has been cultivated since Neolithic times [24]. The region contains some of the northernmost areas in the world where barley can be cultivated and in the far north of the Nordic region, barley has been the only crop that can be grown successfully. The short growing season in the Nordic countries imposes particular selection pressures on cultivated crops, and natural selection for early maturation has been suggested as a force conserving the NAM-B1 wildtype allele in many Nordic spring wheats [25, 26]. It is not known whether natural selection has affected HvNAM genes similarly.

Both animal feed and malt have been important uses of barley in the Nordic countries during the twentieth century [27]. Consequently, barley varieties have primarily been developed based on agronomic and malting quality traits with much less attention paid to food or feed quality [28]. In wheat, selection for yield was probably responsible for the gradual decrease in the frequency of the NAM-B1 wild-type allele in Nordic varieties released during the twentieth century [26]. Should HvNAM genes have a pleiotropic effect similar to that of NAM-B1 a corresponding shift in genetic diversity can be expected.

In this study we have investigated the genetic diversity of HvNAM-1 and HvNAM-2 in Nordic barley, a part of the global germplasm where the HvNAM genes have not been studied previously. The aim was to correlate it to differences in protein and micronutrient content as well as grain size. We also investigated whether the expression of the two HvNAM genes showed any correlation with protein and micronutrient content or grain size. In addition, we were interested in how nutritional and yield traits changed in Nordic barley varieties released during the course of the twentieth century.

Results

Nutritional content and its association with Thousand Grain Weight

Nutritional content was measured in 80 Nordic accessions from four different plant improvement periods and in a high and a low GPC control accession (Table 1). Grain protein content, measured as percentage dry matter, ranged from 6.81% for the accession NGB11311, released in 1990, to 15.25% for the accession NGB9424, released in 1897 (Table 1). Iron (Fe), measured as mg/kg dry matter, ranged from 27.23 for the accession NGB303, released in 1966 to 53.83 for the accession NGB9424. Zinc (Zn), also measured as mg/kg dry matter, ranged from 32.72 for the accession NGB2075, released in 1918, to 49.39 for the accession NGB13022, released in 1994.

Table 1 Information about the accessions used in the study. Measurements of nutrient content in grains, thousand grain weight (TGW) and plant height are averages of three replicates

The low GPC control accession, CIho15487, had the second lowest protein content among the accessions. In contrast, the high control accession, CIho15856, had a protein content that was lower than the mean protein content of all accessions. The high control accession had a Fe content that was in the highest quartile but the low control accession had a Fe content that was higher than the mean among all accessions. For Zn both control accessions had values that laid between the first and the third quartile.

Thousand Grain Weight (TGW) was negatively correlated with protein and Zn content in two-row accessions (protein c = -0.652, p < 0.001; Zn c = -0.570, p < 0.001) and it was positively correlated with Fe (c = 0.421, p < 0.01) and Zn (c = 0.681, p < 0.001) in six-row accessions. Average plant height was in six-row accessions negatively correlated with both Zn (c = -0.404, p < 0.05) and protein (c = -0.366, p < 0.05). Remaining nutritional measurements were not significantly correlated with either TGW or average plant height (all p > 0.05).

Agronomic trait and nutritional differences among improvement stages and row-types

We next compared how agronomic traits and nutrient content varied among improvement groups and row-types in Generalized Linear Model (GLM) analyses (Fig. 1). TGW was significantly higher in both landraces (mean 44.92 g, p < 0.05) and in late improvement accessions (mean 49.96 g, p < 0.001) than in the early improvement accessions (mean 40.38 g), and higher in two-row (mean 46.67 g) than in six-row barley (43.30 g, p < 0.05) in a GLM including both improvement group and row type (glm(TGW ~ improvement_group + row_type)) (Fig. 1a).

Fig. 1
figure 1

Boxplot of improvement groups and row types for a thousand grain weight; b plant height; c protein; d Fe e Zn. Improvement groups consist of landraces and cultivars belonging to an old (released 1890–1940), mid (released 1941–1970) and late (released 1971- 2000) group respectively. The boxes display median value and the first and third quartile, the whiskers the maximum and minimum value

Plant height was lowest in the accession NGB4011 (46.00 cm) and highest in the accession NGB6929 (99.00 cm). In the corresponding GLM of average plant height, height differed between two-row (mean 66.27) and six-row barley (mean 75.46, p < 0.01) and was significantly lower in mid (mean 66.15; p < 0,01) and late (mean 55.51; p < 0.001) improvement accessions than in early improvement accessions (mean 79,00) (Fig. 1b). Neither protein, Fe nor Zn content differed between improvement stages or row-types in GLM analyses (all p > 0.05) (Fig. 1c, d, e).

Haplotypes and genetic diversity

Sequencing of the complete HvNAM-2 gene revealed that all the Nordic accessions were genetically identical to the high GPC control accession CIho15856. HvNAM-1 also had low genetic diversity, but three different haplotypes were detected among the Nordic HvNAM-1 sequences (Table 2). The majority of the sequenced accessions were identical to the high GPC control accession, CIho15856. In the two additional haplotypes, NGB9343 (an accession from the early improvement group) carried a C at position 1063 and NGB6929 (a landrace accession) carried an A at position 1167 using the coordinate system of GenBank accession DQ869678.1. The former of the two substitutions was synonymous (Pro—> Pro) while the second resulted in a tolerated amino acid substitutions (Gly—> Asp), with a SIFT score of 0.15. While NGB6929 had average mineral and protein content and kernel size, values for NGB9343 fell below the lowest quartile for N, Fe, Zn and kernel size. The three differing nucleotides in the low-GPC control accession CIho15487 were not found in any other accession.

Table 2 Consensus sequence and variable positions for three detected haplotypes in HvNAM-1

Gene expression

Gene expression of HvNAM-1 and HvNAM-2 was determined in a subset of eight accessions, representing low and high grain protein content, using qRT-PCR. In the flag leaf, probed 22 days after spike emergence, HvNAM-1 was expressed significantly more (had lower Cq values) than HvNAM-2 (paired t-test, p < 0.001). There was no significant correlation between the expression levels of the two genes (p = 0.16, c = 0.23). Among the analysed accessions, HvNAM-1 had the highest expression (lowest Cq values) in the low protein content accession NGB9942 and the lowest expression in the high protein content accession NGB4613 (Fig. 2). HvNAM-2 had the highest expression in the low protein content accession NGB2070 and the lowest in another low protein content accession, NGB2105. The expression levels of neither HvNAM-1 nor HvNAM-2 were correlated with either nutrient content (protein, Fe, Zn), or TGW (all p > 0.05).

Fig. 2
figure 2

Relative expression of a HvNAM-1 and b HvNAM-2 respectively in flag leaves probed 22 days after spike emergence. Dark grey boxes denote high protein content accessions and light grey boxes denote low protein content accessions. For details see Supplementary Table 1

Discussion

Nordic barley improvement changes during the twentieth century

Barley was one of the first target species when modern plant improvement was initiated in the Nordic region in the 1890s [27, 29]. Improvement of two- and six-row barley have since been carried out independently of each other, with the former being preferred by the malting industry and for cultivation in the south, and the latter primarily as feed barley and for more northern cultivation areas [27]. In contrast with the six-row barley grown in the north, breeding of malting barley also needed to meet quality demands from the malting industry such as low protein, high extraction capacity and even germination ability [28]. Among the most recent cultivars, GPC is also lower in two-row than in six-row varieties. Of the traits measured in this study, TGW and plant height were the only traits for which the row-types differed among the improvement groups, with larger grain and shorter straw-length for two-row barleys. Neither of the values for nutritional contents showed any significant differences in GLM analysis.

A reduction in plant height has been an improvement target not only in barley but also for other cereals during the twentieth century [30]. Plant height in our set of Nordic barley, although following this general trend, was not consistently reduced when comparing more recent improvement groups with older ones. Bertholdsson and Brantestam [31], however, showed that much of the detected effects of plant improvement occurred late in the history of plant improvement, which is supported here by shorter plant heights in accessions from the late improvement group, in both row type groups. In our set of accessions, really tall varieties (> 90 cm) were primarily released before 1950. Accessions released during the second half of the twentieth century were more uniform in height with ever shorter straw lengths (Fig. 1b).

In wheat, worldwide breeding from the 1960s and onwards drastically reduced straw length and increased harvest index and grain yield [32]. As grain yield is often negatively correlated with GPC in wheat [33,34,35,36,37,38], the raised harvests came at the cost of lower GPC, as well as lower and micronutrient content [39, 40]. Similar trends for GPC have been reported in Swedish spring wheats [17].

In our set of Nordic barley we could not detect any differences in nutrient content among the improvement groups, or between row-types. As in wheat, harvest index increased dramatically during twentieth century Nordic barley breeding [31]. Although we were unable to obtain estimates of yield, we note high TGW and short straw length in the late improvement group. Our results indicate that high TGW has not been obtained at the cost of lower GPC. We found no consistent negative correlation between nutritional value and TGW in the accessions studied, although TGW was negatively correlated with protein and Fe in two-row barley.

Links between HvNAM genetic diversity and nutrient content

We found low diversity of HvNAM in our material, in agreement with Jamar et al. [9], in their screen of European cultivars, but contradicting the Wang et al. [41] study of HvNAM-1 in European landrace barley. All but two of the Nordic accessions carried the haplotype Hap7 of Wang et al. [41], the second most common haplotype of that study. The lack of diversity in our material reflects the homogeneity of Nordic barley in a European context [42]. The low genetic diversity prevented us from investigating associations between genetic diversity in either HvNAM gene and nutrient content. Furthermore, the only accessions where SNPs were observed are unlikely to elucidate the function of HvNAM-1. In NGB9343, the SNP did not cause an amino acid substitution and thus, although this accession had low protein, Fe and Zn, we cannot link this to variation in HvNAM-1. The SNP in NGB6929 caused a tolerated aminio acid substitution, but had intermediate grain protein and micronutrient content. The substitution hence caused no detectable phenotypic change in these traits.

Failing to find genetic diversity in eleven studied barley accessions Jamar et al. [9] suggested that differences in protein content could be caused by differential expression of in particular HvNAM-1. Similar to Distelfeld et al. [21] our analysis of the expression of the HvNAM genes, carried out in a subset of our accessions, showed higher expression of HvNAM-1 than HvNAM-2. However, we could not detect any significant association between gene expression and nutrient content.

Evolution of the HvNAM genes

The low genetic diversity of the HvNAM genes meant that no loss of genetic diversity during twentieth century breeding, similar to NAM-B1 in wheat as described by [17] could be seen. The low diversity, in particular in the intronic regions which were invariant, of both HvNAM genes tentatively suggest a past selective sweep, similar to the one described for NAM-B1 in wheat [43]. What the selected trait(s) could have been can only be speculated upon, but if selection was at least moderately strong, the sweep is expected to be quite wide as a consequence of the selfing habit of barley, possibly containing multiple genes. It is not unlikely that regulatory regions of the two genes also show limited genetic diversity and that the differences in expression are expected to be minor. What little diversity was detected occurred in a landrace and an early accession tentatively suggesting more diversity could be present in unimproved Nordic barley. In their more variable set of landrace barley, Wang et al. [41] were, however, unable to detect any statistically significant signals of selection on HvNAM-1.

The GPC of 49 of the 80 the accessions studied here was higher than that of the high protein control accession, CIho15856, used in this and other studies of GPC in barley [7, 8, 10]. Only a single accession, NGB11311, had a lower protein content than the low protein control. This is surprising given that barley improvement schemes in Nordic primarily have targeted malting rather than feed quality [28] and that high protein content is not necessarily desired for malting. In wheat it has been suggested that selection for fast maturation time acting on NAM-B1 has resulted in preservation of the high protein content allele [17, 26]. Earlier senescence in high GPC accessions has similarly been suggested as a cause for negative correlations between GPC and yield in barley [22]. Many barley NAC genes, including HvNAM-1 are also active in the senescence process [20, 23] and selection for early maturation may have pleiotropically kept protein content high. This, however, seems unlikely to be the full explanation. We note that the Hap7 of Wang et al. [41], the most common one in our Nordic barley, is frequently occurring also in climates where much longer growing seasons are expected than in the Nordic region [41] and where fast maturation is not needed to the same extent.

Although HvNAM-1 and HvNAM-2 were almost completely void of genetic diversity, our panel of accessions was not invariant in terms of GPC or micronutrients. It is likely that at least some nutrient variation is caused by genetic effects, although it is clear that the differences detected cannot be attributed to genetic variation in the HvNAM genes. Instead, it is likely that evolution in other loci has given rise to these differences. In their GWAS of protein content Cai et al. [6] not only found associations with HvNAM-1 on chromosome 6H and HvNAM-2 on chromosome 2H. In addition, multiple markers, in many instances within 10 cM of each other, on chromosomes 1H, 3H, 5H and 7H were also found to be associated with protein content. QTL mapping of protein content in Australian and North American two-row barley has similarly showed the trait to be under polygenic control with QTLs identified mainly on 2H, 4H, 5H and 6H [3]. Future studies should evaluate these regions in Nordic barley and explore the material for associations between markers and protein content, to gain a more complete picture of the evolution of differences in nutrient variation.

Conclusions

We find reduced straw length and increased harvest index during twentieth century Nordic barley breeding but no corresponding reduction of grain protein or micronutrient (Zn, Fe) content during the same time period. The genes HvNAM-1 and HvNAM-2 have very little genetic diversity in the Nordic barley gene pool and their expression was not correlated with grain nutrient content. The material nevertheless does show substantial variation in grain nutrient content that could be exploited for further breeding. Genetic variation of HvNAM-1 and HvNAM-2 could also be introduced from other sources, or altered by mutation breeding [44], to contribute additional genetic control of grain protein content.

Methods

Plant material

Eighty accessions of barley from the Nordic region, selected from the set investigated by Kolodinska Brantestam et al. [45], were obtained from NordGen (the Nordic Genetic Resource Center) (Table 1). The accessions consisted of two- and six-row barley from four improvement groups: Landraces and cultivars belonging to an old (released 1890–1940), a mid (released 1941–1970) and a late (released 1971- 2000) group respectively. Accessions from all Nordic countries (with the exception of Iceland) were equally represented within each improvement group.

Two barley varieties, CIho15487 (‘Karl’) and CIho15856 (‘Lewis’), received from the Germplasm Research facility of the Small Grains Collection, US were used as controls. CIho15487 has low GPC while CIho15856 has high GPC. CIho15487 and CIho15856 are known to differ in two SNPs (Single Nucleotide Polymorphism) located in HvNAM-1 [10].

From the set of accessions, four high protein accessions (NGB4613, NGB9343, NGB16752 and NGB9944) and four low protein accessions (NGB2070, NGB2105, NGB2072 and NGB9942) were chosen for expression analyses of HvNAM-1 and HvNAM-2. These eight accessions represented all Nordic countries and three of the four improvement groups. The oldest cultivar group was represented by a single accession.

Plant cultivation

The cultivation experiment was laid out in a randomized block design with three replicates and conducted in a greenhouse with supplementary lighting to create a daylength of 16 h. Day temperature was set to 22 °C and night temperature to 16 °C. The plants were grown in pots 12 × 12 × 12 cm in size, filled with a fertilized peat compost mixture originally containing ~ 200 mg N per pot and essential macro- and micronutrients (g m−3, N 140, P 85, K 150, Fe 4.7, Zn 1.2). Three seeds were sown per pot and reduced to one plant per pot after 17 days. After 17 days the plants were also supplied with additional fertilizer (Algomin Professional Gardening) containing (per pot) 40 mg N, 50 mg P and 110 mg K. The plants were watered with tap water as needed to keep the soil moist throughout the experiment.

Height was measured from the soil surface to the base of the primary spike at full maturation. Grains were harvested manually at full maturity. Three spikes from each plant (less if enough spikes were not available) were harvested. The grains were then dried in a heating cupboard for 48 h at 60 °C and the dry weight was measured. The dried grains were threshed to remove all chaff and counted to obtain 30 grains for estimation of thousand grain weight (TGW). The grain samples were analysed for nitrogen (Leco Corp, Lakeview Avenue, United States), iron and zinc content (Inductively Coupled Plasma atomic emission spectrometry) at Agrilab AB in Uppsala, Sweden. Nitrogen content was converted to protein content through multiplication with 6.25.

The eight accessions used for RNA analysis (NGB4613, NGB2072, NGB9343, NGB2070, NGB2105, NGB16752, NGB9942, NGB9944) were grown as above in single pots, one plant per pot, with six replicates placed in a complete randomized design and surrounded by one row of plants to minimize edge effects. The time points for spike emergence were noted and 22 days later the flag leaves were collected from all replicates, flash frozen in liquid nitrogen and stored at -80 °C until used.

Genetic analysis

DNA was extracted from the different accessions using the Qiagen DNAeasy Plant Mini Kit (Qiagen, Germany) according to the manufacturer's recommendation. HvNAM-1 and HvNAM-2 were amplified in all accessions using PCR. The PCR reaction contained 1 µl of template DNA, 16 µl MilliQ water, 2 µl 10 × DreamTaq Buffer (Thermo Scientific), 0.4 µl 10 mM dNTP (Thermo Scientific), 0.2 µl of 10 µM of each forward and reverse primer and 0.2 µl DreamTaq polymerase (Thermo Scientific). The primers used for amplification of HvNAM-1 were: HvNAM1F 5’-ATGGGCAGCCCGGACTCATCCTCC-3’ and HvNAM1R 5’-TACAGGGATTCCAGTTCACGCCGGAT-3' and for HvNAM-2: HvNAM2F 5’-ATGGGCAGCTCGGACTCATCTTCC-3’ and HvNAM2R 5’-TCAGGGATTCCAGTTCACGCCGGA-3’. These primers, together with internal sequencing primers (HvNAM1SeqF 5’-GCATGAGTACCGCCTCAC-3’ and HvNAM1SeqR 5’-GTGAGGCGGTACTCATGC-3’ for HvNAM-1 and HvNAM2SeqF 5’-GCAGTAACCGATCTCCGTATTT-3’ and HvNAM2SeqR 5’-GGAGATCGGTTACTGCTTGAC-3’ for HvNAM-2, respectively) were used for Sanger sequencing. Sequencing was carried out by Macrogen Europe, the Netherlands.

RNA was isolated from frozen (-80 °C) barley flag leaf tissue using RNeasy Mini Kit as described in [46]. The extracted RNA was diluted to 20 ng/µl and samples with a lower concentration were discarded, leaving four to six biological replicates per accession (Supplementary Table 1). The primers used for qPCR were for HvNAM-1: HvNAM-1QF 5'- CGGCAGTATGTCGCTGTCATCC-3' and HvNAM-1QR 5'-ATGGCGTTCACCGCATTGCC-3', for HvNAM-2: HvNAM-2QF 5'-CGTATGCCACAGCGTGCATGAC-3' and HvNAM-2QR 5'-CTGGTGATGGAGCAGTGAAGCG-3' and for SplFact2: SplFact2QF 5'- GAAGGATGAGTAGGCGCTGG-3' and SplFact2QR 5'-CTGGGAGGTTCCCAACGTAA-3'. QuantiFast SYBR Green RT-PCR (Qiagen AB, Germany) kit was used for the one step qPCR with the following cycling conditions: 10 min at 50 °C for reverse transcription followed by 5 min at 95 °C for Taq activation, and thereafter 40 cycles of 10 s at 95 °C followed by 30 s at 60 °C. To verify for absence of amplifying DNA in the RNA samples negative controls with no reverse transcriptase added in the reaction mixture were run for all samples. Negative controls without template were also run for each primer combination. Amplification of HvNAM-1, HvNAM-2 and the reference gene SplFact2 were run separately but with two repeated reactions for cross referencing between plates. Negative control qPCRs all resulted in no amplification. Melt curve analyses were normal except for one replicate from the accession NGB4613, which was left out from downstream analyses.

Data analysis

DNA sequences were aligned and manually corrected using the software Geneious (v 6.0.5). SIFT scores for predicting effects of nonsynonomous sequence variation was calculated using the SIFT web server [47]. R (v 3.5.1) was used to create plots and for performing statistical analyses. Data on thousand grain weight (TGW) and nutrient content could not be obtained for the accessions NGB2658 and NGB13068 and these accessions were hence left out of analyses of TGW and nutrient content. Nutrient measurements carried out on samples weighing less than 30 g were also excluded. Correlations were tested using the cor.test command. Generalized Linear Models (GLM, glm command) was used to test for differences between row types and improvement groups.

For the qPCR data analysis, the LC480Conversion program (https://www.medischebiologie.nl/files/) was used to convert the raw expression data from the Light Cycler 480 (Roche). The converted data was further analysed with LinRegPCR program to obtain PCR efficiencies and Cq values. Cq values were calibrated across plates and then normalized against the SplFact2 gene expression using GenEx (MultID Analyses, Sweden).

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the GenBank repository, https://www.ncbi.nlm.nih.gov under accession numbers OM416539 - OM416620 (HvNAM-1) and OM416621 - OM416702 (HvNAM-2).

Change history

  • 14 May 2022

    Following the original publication of this article, the funding declaration of the article has been corrected with the addition of the following note: Open access funding provided by Stockholm University.

References

  1. Newton AC, Flavell AJ, George TS, Leat P, Mullholland B, Ramsay L, et al. Crops that feed the world 4. Barley: a resilient crop? strengths and weaknesses in the context of food security. Food Secur. 2011;3(2):141–78.

    Article  Google Scholar 

  2. Briggs DE. Malts and malting: Springer Science & Business Media; 1998.

  3. Emebiri LC, Mood DB, Horsley R, Panozzo J, Read BJ. The genetic control of grain protein content variation in a doubled haploid population derived from a cross between Australian and North American two-rowed barley lines. J Cereal Sci. 2005;41(1):107–14.

    CAS  Article  Google Scholar 

  4. Eagles HA, Bedggood AG, Panozzo JF, Martin PJ. Cultivar and environmental-effects on malting quality in Barley. Aust J Agric Res. 1995;46(5):831–44.

    CAS  Article  Google Scholar 

  5. Bertholdsson NO. Characterization of malting barley cultivars with more or less stable grain protein content under varying environmental conditions. Eur J Agron. 1999;10(1):1–8.

    CAS  Article  Google Scholar 

  6. Cai S, Yu G, Chen X, Huang Y, Jiang X, Zhang G, et al. Grain protein content variation and its association analysis in barley. BMC Plant Biol. 2013;13(35):35.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. See D, Kanazin V, Kephart K, Blake T. Mapping genes controlling variation in barley grain protein concentration. Crop Sci. 2002;42(3):680–5.

    CAS  Article  Google Scholar 

  8. Mickelson S, See D, Meyer FD, Garner JP, Foster CR, Blake TK, et al. Mapping of QTL associated with nitrogen storage and remobilization in barley (Hordeum vulgare L.) leaves. J Exp Bot. 2003;54(383):801–12.

    CAS  PubMed  Article  Google Scholar 

  9. Jamar C, Loffet F, Frettinger P, Ramsay L, Fauconnier ML, du Jardin P. NAM-1 gene polymorphism and grain protein content in Hordeum. J Plant Physiol. 2010;167(6):497–501.

    CAS  PubMed  Article  Google Scholar 

  10. Distelfeld A, Korol A, Dubcovksy J, Uauy C, Blake T, Fahima T. Colinearity between the barley grain protein content (GPC) QTL on chromosome arm 6HS and the wheat Gpc-B1 region. Mol Breed. 2008;22:25–38.

    CAS  Article  Google Scholar 

  11. Uauy C, Brevis JC, Dubcovsky J. The high grain protein content gene Gpc-B1 accelerates senescence and has pleiotropic effects on protein content in wheat. J Exp Bot. 2006;57(11):2785–94.

    CAS  PubMed  Article  Google Scholar 

  12. Brevis JC, Dubcovsky J. Effects of the chromosome region including the locus on wheat grain and protein yield. Crop Sci. 2010;50(1):93.

    CAS  Article  Google Scholar 

  13. Brevis JC, Morris CF, Manthey F, Dubcovsky J. Effect of the grain protein content locus Gpc-B1 on bread and pasta quality. J Cereal Sci. 2010;51(3):357–65.

    CAS  Article  Google Scholar 

  14. Kade M, Barneix AJ, Olmos S, Dubcovsky J. Nitrogen uptake and remobilization in tetraploid “Langdon” durum wheat and a recombinant substitution line with the high grain protein gene Gpc-B1. Plant Breed. 2005;124:343–9.

    CAS  Article  Google Scholar 

  15. Distelfeld A, Cakmak I, Peleg Z, Ozturk L, Yazici AM, Budak H, et al. Multiple QTL-effects of wheat Gpc-B1 locus on grain protein and micronutrient concentrations. Physiol Plant. 2007;129(3):635–43.

    CAS  Article  Google Scholar 

  16. Uauy C, Distelfeld A, Fahima T, Blechl A, Dubcovsky J. A NAC Gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science. 2006;314(5803):1298–301.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. Asplund L, Bergkvist G, Leino MW, Westerbergh A, Weih M. Swedish Spring Wheat Varieties with the Rare High Grain Protein Allele of NAM-B1 Differ in Leaf Senescence and Grain Mineral Content. PloS One. 2013;8(3):e59704.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. Velu G, Singh RP, Cardenas ME, Wu B, Guzman C, Ortiz-Monasterio I. Characterization of grain protein content gene (GPC-B1) introgression lines and its potential use in breeding for enhanced grain zinc and iron concentration in spring wheat. Acta Physiol Plant. Acta Physiol Plant. 2017;39(9):1–9.

  19. Fritsch L, Soeur R, Hansen C, Fischer R, Schillberg S, Schröper F. Next-generation sequencing of amplicons is a rapid and reliable method for the detection of polymorphisms relevant for barley breeding. Mol Breed. 2016;36(7):1–9.

    CAS  Article  Google Scholar 

  20. Christiansen MW, Holm PB, Gregersen PL. Characterization of barley (Hordeum vulgare L.) NAC transcription factors suggests conserved functions compared to both monocots and dicots. BMC Res Notes. 2011;4(1):302.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. Distelfeld A, Avni R, Fischer AM. Senescence, nutrient remobilization, and yield in wheat and barley. J Exp Bot. 2014;65(14):3783–98.

    PubMed  Article  Google Scholar 

  22. Jukanti AK, Fischer AM. A high-grain protein content locus on barley (Hordeum vulgare) chromosome 6 is associated with increased flag leaf proteolysis and nitrogen remobilization. Physiol Plant. 2008;132(4):426–39.

    CAS  PubMed  Article  Google Scholar 

  23. Christiansen MW, Gregersen PL. Members of the barley NAC transcription factor gene family show differential co-regulation with senescence-associated genes during senescence of flag leaves. JExp Bot. 2014;65(14):4009–22.

    CAS  Google Scholar 

  24. Sørensen L, Karg S. The expansion of agrarian societies towards the north–new evidence for agriculture during the Mesolithic/Neolithic transition in Southern Scandinavia. J Archaeol Sci. 2014;51:98–114.

    Article  Google Scholar 

  25. Asplund L, Hagenblad J, Leino MW. Re-evaluating the history of the wheat domestication gene NAM-B1 using historical plant material. J Archaeol Sci. 2010;37(9):2303–7.

    Article  Google Scholar 

  26. Hagenblad J, Asplund L, Balfourier F, Ravel C, Leino MW. Strong presence of the high grain protein content allele of NAM-B1 in Fennoscandian wheat. Theor Appl Genet. 2012;125(8):1677–86.

    CAS  PubMed  Article  Google Scholar 

  27. Persson G. Förädling av vårkorn i Sverige från förr till nu. In: Olsson G, editor. Den svenska växtförädlingens historia. Skogs- och lantbrukshistoriska meddelanden 20. Stockholm: KSLA; 1997. p. 187-204.

  28. Kolodinska Brantestam A. A century of breeding-is genetic erosion a reality? Diss. Alnarp: SLU; 2005.

    Google Scholar 

  29. Fischbeck G. Diversification through breeding. In: Von Bothmer R, van Hintum T, Knüpffer H, Sato K, editors. Diversity in barley (Hordeum vulgare). Elsevier; 2003. p. 29–52.

    Chapter  Google Scholar 

  30. Feil B. Breeding progress in small grain cereals - a comparison of old and modern cultivars. Plant Breed. 1992;108(1):1–11.

    Article  Google Scholar 

  31. Bertholdsson NO, Brantestam AK. A century of Nordic barley breeding-effects on early vigour root and shoot growth, straw length, harvest index and grain weight. Eur J Agron. 2009;30(4):266–74.

    Article  Google Scholar 

  32. Hedden P. The genes of the Green Revolution. Trends Genet. 2003;19(1):5–9.

    CAS  PubMed  Article  Google Scholar 

  33. Kibite S, Evans LE. Causes of Negative Correlations between Grain-Yield and Grain Protein-Concentration in Common Wheat. Euphytica. 1984;33(3):801–10.

    Article  Google Scholar 

  34. Peltonen-Sainio P, Jauhiainen L, Nissila E. Improving cereal protein yields for high latitude conditions. Eur J Agron. 2012;39:1–8.

    Article  Google Scholar 

  35. Martre P, He J, Le Gouis J, Semenov MA. In silico system analysis of physiological traits determining grain yield and protein concentration for wheat as influenced by climate and crop management. J Exp Bot. 2015;66(12):3581–98.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. Bogard M, Jourdan M, Allard V, Martre P, Perretant MR, Ravel C, et al. Anthesis date mainly explained correlations between post-anthesis leaf senescence, grain yield, and grain protein concentration in a winter wheat population segregating for flowering time QTLs. J Exp Bot. 2011;62(10):3621–36.

    CAS  PubMed  Article  Google Scholar 

  37. Oury F-X, Godin C. Yield and grain protein concentration in bread wheat: how to use the negative relationship between the two characters to identify favourable genotypes? Euphytica. 2007;157(1–2):45–57.

    CAS  Article  Google Scholar 

  38. Blanco A, Mangini G, Giancaspro A, Giove S, Colasuonno P, Simeone R, et al. Relationships between grain protein content and grain yield components through quantitative trait locus analyses in a recombinant inbred line population derived from two elite durum wheat cultivars. Mol Breed. 2011;30(1):79–92.

    Article  CAS  Google Scholar 

  39. Garvin DF, Welch RM, Finley JW. Historical shifts in the seed mineral micronutrient concentration of US hard red winter wheat germplasm. J Sci Food Agric. 2006;86(13):2213–20.

    CAS  Article  Google Scholar 

  40. Fan MS, Zhao FJ, Fairweather-Tait SJ, Poulton PR, Dunham SJ, McGrath SP. Evidence of decreasing mineral density in wheat grain over the last 160 years. J Trace Elem Med Biol. 2008;22(4):315–24.

    CAS  PubMed  Article  Google Scholar 

  41. Wang Y, Ren X, Sun D, Sun G. Origin of worldwide cultivated barley revealed by NAM-1 gene and grain protein content. Frontiers Plant Sci. 2015;6:803.

    Google Scholar 

  42. Aslan S, Forsberg N, Hagenblad J, Leino MW. Molecular genotyping of historical barley landraces reveals novel candidate regions for local adaption. Crop Sci. 2015;55(6):2766–76.

    CAS  Article  Google Scholar 

  43. Lundström M, Leino MW, Hagenblad J. Evolutionary history of the NAM-B1 gene in wild and domesticated tetraploid wheat. BMC Genet. 2017;18(1):118.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  44. Chapman EA, Orford S, Lage J, Griffiths S. Delaying or delivering: identification of novel NAM-1 alleles that delay senescence to extend wheat grain fill duration. J Exp Bot. 2021. https://doi.org/10.1093/jxb/erab368.

    Article  PubMed  PubMed Central  Google Scholar 

  45. KolodinskaBrantestam A, Bothmer R, Dayteg C, Rashal I, Tuvesson S, Weibull J. Genetic diversity changes and relationships in spring barley (Hordeum vulgare L.) germplasm of Nordic and Baltic areas as shown by SSR markers. Genet Res Crop Evol. 2007;54(4):749–58.

    Article  Google Scholar 

  46. Vanhala T, Normann KR, Lundström M, Weller JL, Leino MW, Hagenblad J. Flowering time adaption in Swedish landrace pea (Pisum sativum L.). BMC genetics. 2016;17:117.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  47. Sim N-L, Kumar P, Hu J, Henikoff S, Schneider G, Ng PC. SIFT web server: predicting effects of amino acid substitutions on proteins. Nucleic Acids Res. 2012;40(W1):W452–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

Sandra Lilja, Sara Ödling, Amanda Norberg and an anonymous student are all acknowledged for their contributions towards the sequencing of HvNAM-1 and HvNAM-2.

Funding

The authors gratefully acknowledge the funding contributed by the Lagersberg Foundation, the Carl Tryggers Foundation (grants CTS 14:169 and CTS 15:189) and the C.F. Lundström Foundation (grant CF2017-0009). The funding bodies played no part in the design of the study, analysis, interpretation of data or in writing the manuscript. Open access funding provided by Stockholm University.

Author information

Affiliations

Authors

Contributions

MWL designed the study together with JH. SM, JH and TV carried out the lab work. JH analysed the majority of the data with additional analyses done by TV. JH wrote the paper together with MWL. All authors read and approved the final manuscript.

Authors’ information

Not applicable

Corresponding author

Correspondence to Matti W. Leino.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Hagenblad, J., Vanhala, T., Madhavan, S. et al. Protein content and HvNAM alleles in Nordic barley (Hordeum vulgare) during a century of breeding. Hereditas 159, 12 (2022). https://doi.org/10.1186/s41065-022-00227-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s41065-022-00227-y

Keywords

  • HvNAM
  • Resequencing
  • Grain protein content
  • Hordeum vulgare
  • Crop improvement