Association mapping of malting quality traits in UK spring and winter barley cultivar collections

Unlocking the genetic diversity and population structure of the newly introduced two-row spring European HerItage Barley collecTion (ExHIBiT)

In the last century, breeding programs have traditionally favoured yield-related traits, grown under high-input conditions, resulting in a loss of genetic diversity and an increased susceptibility to stresses in crops. Thus, exploiting understudied genetic resources, that potentially harbour tolerance genes, is vital for sustainable agriculture. Northern European barley germplasm has been relatively understudied despite its key role within the malting industry. The European Heritage Barley collection (ExHIBiT) was assembled to explore the genetic diversity in European barley focusing on Northern European accessions and further address environmental pressures. ExHIBiT consists of 363 spring-barley accessions, focusing on two-row type. The collection consists of landraces (~14%), old cultivars (~18%), elite cultivars (~67%) and accessions with unknown breeding history (~1%), with 70% of the collection from Northern Europe. The population structure of the ExHIBiT collection was subdivided into three main clusters primarily based on the accession's year of release using 26,585 informative SNPs based on 50k iSelect single nucleotide polymorphism (SNP) array data. Power analysis established a representative core collection of 230 genotypically and phenotypically diverse accessions. The effectiveness of this core collection for conducting statistical and association analysis was explored by undertaking genome-wide association studies (GWAS) using 24,876 SNPs for nine phenotypic traits, four of which were associated with SNPs. Genomic regions overlapping with previously characterised flowering genes (HvZTLb) were identified, demonstrating the utility of the ExHIBiT core collection for locating genetic regions that determine important traits. Overall, the ExHIBiT core collection represents the high level of untapped diversity within Northern European barley, providing a powerful resource for researchers and breeders to address future climate scenarios.

Genomic resources for a historical collection of cultivated two-row European spring barley genotypes

Barley genomic resources are increasing rapidly, with the publication of a barley pangenome as one of the latest developments. Two-row spring barley cultivars are intensely studied as they are the source of high-quality grain for malting and distilling. Here we provide data from a European two-row spring barley population containing 209 different genotypes registered for the UK market between 1830 to 2014. The dataset encompasses RNA-sequencing data from six different tissues across a range of barley developmental stages, phenotypic datasets from two consecutive years of field-grown trials in the United Kingdom, Germany and the USA; and whole genome shotgun sequencing from all cultivars, which was used to complement the RNA-sequencing data for variant calling. The outcomes are a filtered SNP marker file, a phenotypic database and a large gene expression dataset providing a comprehensive resource which allows for downstream analyses like genome wide association studies or expression associations.

Detection of consensus genomic regions and candidate genes for quality traits in barley using QTL meta-analysis

Improving barley grain quality is a major goal in barley breeding. In this study, a total of 35 papers focusing on quantitative trait loci (QTLs) mapping for barley quality traits published since 2000 were collected. Among the 454 QTLs identified in these studies, 349 of them were mapped onto high-density consensus maps, which were used for QTL meta-analysis. Through QTL meta-analysis, the initial QTLs were integrated into 41 meta-QTLs (MQTLs) with an average confidence interval (CI) of 1. 66 cM, which is 88.9% narrower than that of the initial QTLs. Among the 41 identified MQTLs, 25 were subsequently validated in publications using genome-wide association study (GWAS). From these 25 validated MQTLs, ten breeder's MQTLs were selected. Synteny analysis comparing barley and wheat MQTLs revealed orthologous relationships between eight breeder's MQTLs and 45 wheat MQTLs. Additionally, 17 barley homologs associated with rice quality traits were identified within the regions of the breeder's MQTLs through comparative analysis. The findings of this study provide valuable insights for molecular marker-assisted breeding and the identification of candidate genes related to quality traits in barley.

Genomic traces of Japanese malting barley breeding in two modern high-quality cultivars, 'Sukai Golden' and 'Sachiho Golden'

Two modern high-quality Japanese malting barley cultivars, 'Sukai Golden' and 'Sachiho Golden', were subjected to RNA-sequencing of transcripts extracted from 20-day-old immature seeds. Despite their close relation, 2,419 Sukai Golden-specific and 3,058 Sachiho Golden-specific SNPs were detected in comparison to the genome sequences of two reference cultivars: 'Morex' and 'Haruna Nijo'. Two single nucleotide polymorphism (SNP) clusters respectively showing the incorporation of (1) the barley yellow mosaic virus (BaYMV) resistance gene rym5 from six-row non-malting Chinese landrace Mokusekko 3 on the long arm of 3H, and (2) the anthocyanin-less ant2 gene from a two-row Dutch cultivar on the long arm of 2H were detected specifically in 'Sukai Golden'. Using 221 recombinant inbred lines of a cross between 'Ishukushirazu' and 'Nishinochikara', another BaYMV resistance rym3 gene derived from six-row non-malting Japanese cultivar 'Haganemugi' was mapped to a 0.4-cM interval on the proximal region of 5H. Haplotype analysis of progenitor accessions of the two modern malting cultivars revealed that rym3 of 'Haganemugi' was independently introduced into 'Sukai Golden' and 'Sachiho Golden'. Residual chromosome 5H segments of 'Haganemugi' surrounding rym3 were larger in 'Sukai Golden'. Available results suggest possibilities for malting quality improvement by minimizing residual segments surrounding rym3.

Identification of the Key Transcription Factors Regulating the Expression of the Genes Associated with Barley Malt Quality during Malting

Barley malt is produced through a malting process; it begins with steeping followed by germination and kilning, in which dramatic changes happen for a large number of physiological and biochemical traits in barley seeds. The objectives of this study were to comprehensively investigate the phenotypic changes during malting, and identify the key regulators that modulate the expression of genes associated with malt quality traits. The results showed that there was a significant positive correlation between gibberellic acid (GA) content and the activities of some hydrolytic enzymes, including α-amylases, β-amylases, and limit dextrinase (LD), and a significant negative correlation between GA and β-glucan content. Starch content had little change, but starch granules were pitted severely during malting. Weighted gene coexpression analysis (WGCNA) identified the genes associated with the greatest changes of the examined malt traits during malting. The correlation analysis and protein-protein interaction (PPI) analysis detected several key transcriptional factor (TF) regulating genes associated with malt quality. These genes and TFs regulating malting traits are potentially useful in barley breeding for malt quality improvement.

Natural variation in HvAT10 underlies grain cell wall-esterified phenolic acid content in cultivated barley

The phenolic acids, ferulic acid and p-coumaric acid, are components of plant cell walls in grasses, including many of our major food crops. They have important health-promoting properties in grain, and influence the digestibility of biomass for industrial processing and livestock feed. Both phenolic acids are assumed to be critical to cell wall integrity and ferulic acid, at least, is important for cross-linking cell wall components, but the role of p-coumaric acid is unclear. Here we identify alleles of a BAHD p-coumaroyl arabinoxylan transferase, HvAT10, as responsible for the natural variation in cell wall-esterified phenolic acids in whole grain within a cultivated two-row spring barley panel. We show that HvAT10 is rendered non-functional by a premature stop codon mutation in half of the genotypes in our mapping panel. This results in a dramatic reduction in grain cell wall-esterifed p-coumaric acid, a moderate rise in ferulic acid, and a clear increase in the ferulic acid to p-coumaric acid ratio. The mutation is virtually absent in wild and landrace germplasm suggesting an important function for grain arabinoxylan p-coumaroylation pre-domestication that is dispensable in modern agriculture. Intriguingly, we detected detrimental impacts of the mutated locus on grain quality traits where it was associated with smaller grain and poorer malting properties. HvAT10 could be a focus for improving grain quality for malting or phenolic acid content in wholegrain foods.

Malting quality and preharvest sprouting traits are genetically correlated in spring malting barley

Malt for craft "all-malt" brewing can have high quality, PHS resistance, and malted in normal timeframes. Canadian style adjunct malt is associated with PHS susceptibility. Expansion of malting barley production into non-traditional growing regions and erratic weather has increased the demand for preharvest sprouting (PHS) resistant, high quality malting barley cultivars. This is hindered by the relatively unknown relationships between PHS resistance and malting quality. Here we present a three-year study of malting quality and germination at different after-ripening durations post physiological maturity. Malting quality traits alpha amylase (AA) and free amino nitrogen (FAN) and germination rate at six days post PM shared a common association with a SNP in HvMKK3 on chromosome 5H in the Seed Dormancy 2 (SD2) region responsible for PHS susceptibility. Soluble protein (SP) and soluble over total protein (S/T) both shared a common association with a marker in the SD2 region. Significant genetic correlations between PHS resistance and the malting quality traits AA, FAN, SP, S/T were detected across and within HvMKK3 allele groups. High adjunct malt quality was related to PHS susceptibility. Selection for PHS resistance led to a correlated response in malting quality traits. Results strongly suggest pleiotropy of HvMKK3 on malting quality traits and that the classic "Canadian-style" malt is caused by a PHS susceptible allele of HvMKK3. PHS susceptibility appears to benefit the production of malt intended for adjunct brewing, while PHS resistance is compatible with all-malt brewing specifications. Here we present our analysis on the effect of combining complexly inherited and correlated traits with contrasting goals to inform breeding practice in malting barley, the general principles of which can be extended to other breeding programs.

CGIAR Barley Breeding Toolbox: A diversity panel to facilitate breeding and genomic research in the developing world

Breeding programs in developing countries still cannot afford the new genotyping technologies, hindering their research. We aimed to assemble an Association Mapping panel to serve as CGIAR Barley Breeding Toolbox (CBBT), especially for the Developing World. The germplasm had to be representative of the one grown in the Developing World; with high genetic variability and be of public domain. For it, we genotyped with the Infinium iSelect 50K chip, a Global Barley Panel (GBP) of 530 genotypes representing a wide range of row-types, end-uses, growth habits, geographical origins and environments. 40,342 markers were polymorphic with an average polymorphism information content of 0.35 and 66% of them exceeding 0.25. The analysis of the population structure identified 8 subpopulations mostly linked to geographical origin, four of them with significant ICARDA origin. The 16 allele combinations at 4 major flowering genes (HvVRN-H3, HvPPD-H1, HvVRN-H1 and HvCEN) explained 11.07% genetic variation and were linked to the geographic origins of the lines. ICARDA material showed the widest diversity as revealed by the highest number of polymorphic loci (99.76% of all polymorphic SNPs in GBP), number of private alleles and the fact that ICARDA lines were present in all 8 subpopulations and carried all 16 allelic combinations. Due to their genetic diversity and their representativity of the germplasm adapted to the Developing World, ICARDA-derived lines and cultivated landraces were pre-selected to form the CBBT. Using the Mean of Transformed Kinships method, we assembled a panel capturing most of the allelic diversity in the GBP. The CBBT (N=250) preserves good balance between row-types and good representation of both phenology allelic combinations and subpopulations of the GBP. The CBBT and its genotypic data is available to researchers worldwide as a collaborative tool to underpin the genetic mechanisms of traits of interest for barley cultivation.

Analysis of historical selection in winter wheat

KEY MESSAGE

Modeling of the distribution of allele frequency over year of variety release identifies major loci involved in historical breeding of winter wheat. Winter wheat is a major crop with a rich selection history in the modern era of crop breeding. Genetic gains across economically important traits like yield have been well characterized and are the major force driving its production. Winter wheat is also an excellent model for analyzing historical genetic selection. As a proof of concept, we analyze two major collections of winter wheat varieties that were bred in Western Europe from 1916 to 2010, namely the Triticeae Genome (TG) and WAGTAIL panels, which include 333 and 403 varieties, respectively. We develop and apply a selection mapping approach, Regression of Alleles on Years (RALLY), in these panels, as well as in simulated populations. RALLY maps loci under sustained historical selection by using a simple logistic model to regress allele counts on years of variety release. To control for drift-induced allele frequency change, we develop a hybrid approach of genomic control and delta control. Within the TG panel, we identify 22 significant RALLY quantitative selection loci (QSLs) and estimate the local heritabilities for 12 traits across these QSLs. By correlating predicted marker effects with RALLY regression estimates, we show that alleles whose frequencies have increased over time are heavily biased toward conferring positive yield effect, but negative effects in flowering time, lodging, plant height and grain protein content. Altogether, our results (1) demonstrate the use of RALLY to identify selected genomic regions while controlling for drift, and (2) reveal key patterns in the historical selection in winter wheat and guide its future breeding.

Trend, population structure, and trait mapping from 15 years of national varietal trials of UK winter wheat

There are now a rich variety of genomic and genotypic resources available to wheat researchers and breeders. However, the generation of high-quality and field-relevant phenotyping data which is required to capture the complexities of gene × environment interactions remains a major bottleneck. Historical datasets from national variety performance trials (NVPT) provide sufficient dimensions, in terms of numbers of years and locations, to examine phenotypic trends and study gene × environment interactions. Using NVPT for winter wheat varieties grown in the United Kingdom between 2002 and 2017, we examined temporal trends for eight traits related to yield, adaptation, and grain quality performance. We show a non-stationary linear trend for yield, grain protein content, Hagberg Falling Number (HFN), and days to ripening. Our data also show high environmental stability for yield, grain protein content, and specific weight in UK winter wheat varieties and high environmental sensitivity for HFN. We also show that UK varieties released within this period cluster into four main population groups. Using the historical NVPT data in a genome-wide association analysis, we uncovered a significant marker-trait association peak on wheat chromosome 6A spanning the NAM-A1 gene that have been previously associated with early senescence. Together, our results show the value of utilizing the data routinely collected during national variety evaluation process for examining breeding progress and the genetic architecture of important traits.

Trends of genetic changes uncovered by Env- and Eigen-GWAS in wheat and barley

Variety age and population structure detect novel QTL for yield and adaptation in wheat and barley without the need to phenotype. The process of crop breeding over the last century has delivered new varieties with increased genetic gains, resulting in higher crop performance and yield. However, in many cases, the alleles and genomic regions underpinning this success remain unknown. This is partly due to the difficulty of generating sufficient phenotypic data on large numbers of historical varieties to enable such analyses. Here we demonstrate the ability to circumvent such bottlenecks by identifying genomic regions selected over 100 years of crop breeding using age of a variety as a surrogate for yield. Rather than collecting phenotype data, we deployed 'environmental genome-wide association scans' (EnvGWAS) based on variety age in two of the world's most important crops, wheat and barley, and detected strong signals of selection across both genomes. EnvGWAS identified 16 genomic regions in barley and 10 in wheat with contrasting patterns between spring and winter types of the two crops. To further examine changes in genome structure, we used the genomic relationship matrix of the genotypic data to derive eigenvectors for analysis in EigenGWAS. This detected seven major chromosomal introgressions that contributed to adaptation in wheat. EigenGWAS and EnvGWAS based on variety age avoid costly phenotyping and facilitate the identification of genomic tracts that have been under selection during breeding. Our results demonstrate the potential of using historical cultivar collections coupled with genomic data to identify chromosomal regions under selection and may help guide future plant breeding strategies to maximise the rate of genetic gain and adaptation.

Experimental comparison of two methods to study barley responses to partial submergence

BACKGROUND

Crop yield is dependent on climate conditions, which are becoming both more variable and extreme in some areas of the world as a consequence of global climate change. Increased precipitation and flooding events are the cause of important yield losses due to waterlogging or (partial) submergence of crops in the field. Our ability to screen efficiently and quickly for varieties that have increased tolerance to waterlogging or (partial) submergence is important. Barley, a staple crop worldwide, is particularly sensitive to waterlogging. Screening for waterlogging tolerant barley varieties has been ongoing for many years, but methods used to screen vary greatly, from the type of soil used to the time at which the treatment is applied. This variation makes it difficult to cross-compare results.

RESULTS

Here, we have devised a scoring system to assess barley tolerance to waterlogging and compare two different methods when partial submergence is applied with either water or a starch solution at an early developmental stage, which is particularly sensitive to waterlogging or partial submergence. The use of a starch solution has been previously shown to result in more reducing soil conditions and has been used to screen for waterlogging tolerance.

CONCLUSIONS

Our results show that the two methods provide similar results to qualitatively rank varieties as tolerant or sensitive, while also affecting plants differently, in that application of a starch solution results in stronger and earlier symptoms than applying partial submergence with water.