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.

A GDSL-motif Esterase/Lipase Affects Wax and Cutin Deposition and Controls Hull-Caryopsis Attachment in Barley.

The cuticle covering aerial organs of land plants is well known to protect against desiccation. Cuticles also play diverse and specialised functions, including organ separation, depending on plant and tissue. Barley shows a distinctive cuticular wax bloom enriched in beta-diketones on leaf sheaths, stem nodes and internodes, and inflorescences. Barley also develops a sticky surface on the outer pericarp layer of its grain fruit leading to strongly adhered hulls, 'covered grain', important for embryo protection and seed dispersal. While the transcription factor-encoding gene HvNUDUM (HvNUD) appears essential for adherent hulls, little is understood about how the pericarp cuticle changes during adhesion or whether changes in pericarp cuticles contribute to another phenotype where hulls partially shed, called 'skinning'. To that end, we screened barley lines for hull adhesion defects, focussing on the Eceriferum (= waxless, cer) mutants. Here, we show that the cer-xd allele causes defective wax blooms and compromised hull adhesion, and results from a mutation removing the last ten amino acids of the GDSL-motif esterase/lipase HvGDSL1. We used severe and moderate HvGDSL1 alleles to show that complete HvGDSL1 function is essential for leaf blade cuticular integrity, wax bloom deposition over inflorescences and leaf sheaths, and pericarp cuticular ridge formation. Expression data suggests that HvGDSL1 may regulate hull adhesion independently of HvNUD. We found high conservation of HvGDSL1 among barley germplasm, so variation in HvGDSL1 unlikely leads to grain skinning in cultivated barley. Taken together, we reveal a single locus which controls adaptive cuticular properties across different organs in barley.

A guide to barley mutants.

BACKGROUND: Mutants have had a fundamental impact upon scientific and applied genetics. They have paved the way for the molecular and genomic era, and most of today's crop plants are derived from breeding programs involving mutagenic treatments. RESULTS: Barley (Hordeum vulgare L.) is one of the most widely grown cereals in the world and has a long history as a crop plant. Barley breeding started more than 100 years ago and large breeding programs have collected and generated a wide range of natural and induced mutants, which often were deposited in genebanks around the world. In recent years, an increased interest in genetic diversity has brought many historic mutants into focus because the collections are regarded as valuable resources for understanding the genetic control of barley biology and barley breeding. The increased interest has been fueled also by recent advances in genomic research, which provided new tools and possibilities to analyze and reveal the genetic diversity of mutant collections. CONCLUSION: Since detailed knowledge about phenotypic characters of the mutants is the key to success of genetic and genomic studies, we here provide a comprehensive description of mostly morphological barley mutants. The review is closely linked to the International Database for Barley Genes and Barley Genetic Stocks ( bgs.nordgen.org ) where further details and additional images of each mutant described in this review can be found.

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.

Major chromosome 5H haplotype switch structures the European two-rowed spring barley germplasm of the past 190 years.

KEY MESSAGE: Selection over 70 years has led to almost complete fixation of a haplotype spanning ~ 250 Mbp of chomosome 5H in European two-rowed spring barleys, possibly originating from North Africa. Plant breeding and selection have shaped the genetic composition of modern crops over the past decades and centuries and have led to great improvements in agronomic and quality traits. Knowledge of the genetic composition of breeding germplasm is essential to make informed decisions in breeding programs. In this study, we characterized the structure and composition of 209 barley cultivars representative of the European two-rowed spring barley germplasm of the past 190 years. Utilizing high-density SNP marker data, we identified a distinct centromeric haplotype spanning a ~ 250 Mbp large region on chromosome 5H which likely was first introduced into the European breeding germplasm in the early to mid-twentieth century and has been non-recombining and under strong positive selection over the past 70 years. Almost all cultivars in our panel that were released after 2000 carry this new haplotype, suggesting that this region carries one or several genes conferring highly beneficial traits. Using the global barley collection of the German Federal ex situ gene bank at IPK Gatersleben, we found the new haplotype at high frequencies in six-rowed spring-type landraces from Northern Africa, from which it may have been introduced into modern European barley germplasm via southern European landraces. The presence of a 250 Mbp genomic region characterized by lack of recombination and high levels of fixation in modern barley germplasm has substantial implications for the genetic diversity of the modern barley germplasm and for barley breeding.

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.

Conserved signalling components coordinate epidermal patterning and cuticle deposition in barley.

Faced with terrestrial threats, land plants seal their aerial surfaces with a lipid-rich cuticle. To breathe, plants interrupt their cuticles with adjustable epidermal pores, called stomata, that regulate gas exchange, and develop other specialised epidermal cells such as defensive hairs. Mechanisms coordinating epidermal features remain poorly understood. Addressing this, we studied two loci whose allelic variation causes both cuticular wax-deficiency and misarranged stomata in barley, identifying the underlying genes, Cer-g/ HvYDA1, encoding a YODA-like (YDA) MAPKKK, and Cer-s/ HvBRX-Solo, encoding a single BREVIS-RADIX (BRX) domain protein. Both genes control cuticular integrity, the spacing and identity of epidermal cells, and barley's distinctive epicuticular wax blooms, as well as stomatal patterning in elevated CO2 conditions. Genetic analyses revealed epistatic and modifying relationships between HvYDA1 and HvBRX-Solo, intimating that their products participate in interacting pathway(s) linking epidermal patterning with cuticular properties in barley. This may represent a mechanism for coordinating multiple adaptive features of the land plant epidermis in a cultivated cereal.

Multi-environment genome -wide association mapping of culm morphology traits in barley.

In cereals with hollow internodes, lodging resistance is influenced by morphological characteristics such as internode diameter and culm wall thickness. Despite their relevance, knowledge of the genetic control of these traits and their relationship with lodging is lacking in temperate cereals such as barley. To fill this gap, we developed an image analysis-based protocol to accurately phenotype culm diameters and culm wall thickness across 261 barley accessions. Analysis of culm trait data collected from field trials in seven different environments revealed high heritability values (>50%) for most traits except thickness and stiffness, as well as genotype-by-environment interactions. The collection was structured mainly according to row-type, which had a confounding effect on culm traits as evidenced by phenotypic correlations. Within both row-type subsets, outer diameter and section modulus showed significant negative correlations with lodging (<-0.52 and <-0.45, respectively), but no correlation with plant height, indicating the possibility of improving lodging resistance independent of plant height. Using 50k iSelect SNP genotyping data, we conducted multi-environment genome-wide association studies using mixed model approach across the whole panel and row-type subsets: we identified a total of 192 quantitative trait loci (QTLs) for the studied traits, including subpopulation-specific QTLs and 21 main effect loci for culm diameter and/or section modulus showing effects on lodging without impacting plant height. Providing insights into the genetic architecture of culm morphology in barley and the possible role of candidate genes involved in hormone and cell wall-related pathways, this work supports the potential of loci underpinning culm features to improve lodging resistance and increase barley yield stability under changing environments.

Identification of candidate MYB transcription factors that influence CslF6 expression in barley grain.

(1,3;1,4)-ß-Glucan is a non-cellulosic polysaccharide required for correct barley grain fill and plant development, with industrial relevance in the brewing and the functional food sector. Barley grains contain higher levels of (1,3;1,4)-ß-glucan compared to other small grain cereals and this influences their end use, having undesirable effects on brewing and distilling and beneficial effects linked to human health. HvCslF6 is the main gene contributing to (1,3;1,4)-ß-glucan biosynthesis in the grain. Here, the transcriptional regulation of HvCslF6 was investigated using an in-silico analysis of transcription factor binding sites (TFBS) in its putative promoter, and functional characterization in a barley protoplast transient expression system. Based on TFBS predictions, TF classes AP2/ERF, MYB, and basic helix-loop-helix (bHLH) were over-represented within a 1,000 bp proximal HvCslF6 promoter region. Dual luciferase assays based on multiple HvCslF6 deletion constructs revealed the promoter fragment driving HvCslF6 expression. Highest HvCslF6 promoter activity was narrowed down to a 51 bp region located -331 bp to -382 bp upstream of the start codon. We combined this with TFBS predictions to identify two MYB TFs: HvMYB61 and HvMYB46/83 as putative activators of HvCslF6 expression. Gene network analyses assigned HvMYB61 to the same co-expression module as HvCslF6 and other primary cellulose synthases (HvCesA1, HvCesA2, and HvCesA6), whereas HvMYB46/83 was assigned to a different module. Based on RNA-seq expression during grain development, HvMYB61 was cloned and tested in the protoplast system. The transient over-expression of HvMYB61 in barley protoplasts suggested a positive regulatory effect on HvCslF6 expression.

Measuring the frequency and distribution of meiotic crossovers in homozygous barley inbred lines.

We report a novel approach for establishing the number and position of CO events in individual homozygous inbred plants by combining low level EMS mutagenesis, speed breeding, whole genome shotgun sequencing and sliding window analysis of the induced molecular variant data. We demonstrate the approach by exploring CO frequency and distribution in self-fertilised progeny of the inbred barley cultivar Bowman and compare these observations to similar data obtained from a Bowman nearly isogenic line (BW230 Hvmlh3) containing a mutation in the DNA mismatch repair gene HvMLH3. We have previously shown that Hvmlh3 decreases both plant fertility and recombination by ~50%. We compare our results to those from previously published traditional genetic analysis of F3 families derived from multiple F2 lines containing WT or mutant alleles of HvMLH3, revealing a high level of correspondence between analyses. We discuss possible applications of the approach in streamlining the assessment of recombination in plant meiosis research.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis.

BACKGROUND: Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. RESULTS: We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. CONCLUSIONS: AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

The evolutionary patterns of barley pericentromeric chromosome regions, as shaped by linkage disequilibrium and domestication.

The distribution of recombination events along large cereal chromosomes is uneven and is generally restricted to gene-rich telomeric ends. To understand how the lack of recombination affects diversity in the large pericentromeric regions, we analysed deep exome capture data from a final panel of 815 Hordeum vulgare (barley) cultivars, landraces and wild barleys, sampled from across their eco-geographical ranges. We defined and compared variant data across the pericentromeric and non-pericentromeric regions, observing a clear partitioning of diversity both within and between chromosomes and germplasm groups. Dramatically reduced diversity was found in the pericentromeres of both cultivars and landraces when compared with wild barley. We observed a mixture of completely and partially differentiated single-nucleotide polymorphisms (SNPs) between domesticated and wild gene pools, suggesting that domesticated gene pools were derived from multiple wild ancestors. Patterns of genome-wide linkage disequilibrium, haplotype block size and number, and variant frequency within blocks showed clear contrasts among individual chromosomes and between cultivars and wild barleys. Although most cultivar chromosomes shared a single major pericentromeric haplotype, chromosome 7H clearly differentiated the two-row and six-row types associated with different geographical origins. Within the pericentromeric regions we identified 22 387 non-synonymous SNPs, 92 of which were fixed for alternative alleles in cultivar versus wild accessions. Surprisingly, only 29 SNPs found exclusively in the cultivars were predicted to be 'highly deleterious'. Overall, our data reveal an unconventional pericentromeric genetic landscape among distinct barley gene pools, with different evolutionary processes driving domestication and diversification.

The genome and gene editing system of sea barleygrass provide a novel platform for cereal domestication and stress tolerance studies.

The tribe Triticeae provides important staple cereal crops and contains elite wild species with wide genetic diversity and high tolerance to abiotic stresses. Sea barleygrass (Hordeum marinum Huds.), a wild Triticeae species, thrives in saline marshlands and is well known for its high tolerance to salinity and waterlogging. Here, a 3.82-Gb high-quality reference genome of sea barleygrass is assembled de novo, with 3.69 Gb (96.8%) of its sequences anchored onto seven chromosomes. In total, 41 045 high-confidence (HC) genes are annotated by homology, de novo prediction, and transcriptome analysis. Phylogenetics, non-synonymous/synonymous mutation ratios (Ka/Ks), and transcriptomic and functional analyses provide genetic evidence for the divergence in morphology and salt tolerance among sea barleygrass, barley, and wheat. The large variation in post-domestication genes (e.g. IPA1 and MOC1) may cause interspecies differences in plant morphology. The extremely high salt tolerance of sea barleygrass is mainly attributed to low Na+ uptake and root-to-shoot translocation, which are mainly controlled by SOS1, HKT, and NHX transporters. Agrobacterium-mediated transformation and CRISPR/Cas9-mediated gene editing systems were developed for sea barleygrass to promote its utilization for exploration and functional studies of hub genes and for the genetic improvement of cereal crops.

Identifying plant genes shaping microbiota composition in the barley rhizosphere.

A prerequisite to exploiting soil microbes for sustainable crop production is the identification of the plant genes shaping microbiota composition in the rhizosphere, the interface between roots and soil. Here, we use metagenomics information as an external quantitative phenotype to map the host genetic determinants of the rhizosphere microbiota in wild and domesticated genotypes of barley, the fourth most cultivated cereal globally. We identify a small number of loci with a major effect on the composition of rhizosphere communities. One of those, designated the QRMC-3HS, emerges as a major determinant of microbiota composition. We subject soil-grown sibling lines harbouring contrasting alleles at QRMC-3HS and hosting contrasting microbiotas to comparative root RNA-seq profiling. This allows us to identify three primary candidate genes, including a Nucleotide-Binding-Leucine-Rich-Repeat (NLR) gene in a region of structural variation of the barley genome. Our results provide insights into the footprint of crop improvement on the plant's capacity of shaping rhizosphere microbes.

BaRTv2: a highly resolved barley reference transcriptome for accurate transcript-specific RNA-seq quantification.

Accurate characterisation of splice junctions (SJs) as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data, and enables detailed investigations of transcriptional and post-transcriptional regulation. Using novel computational methods and a combination of PacBio Iso-seq and Illumina short-read sequences from 20 diverse tissues and conditions, we generated a comprehensive and highly resolved barley reference transcript dataset from the European 2-row spring barley cultivar Barke (BaRTv2.18). Stringent and thorough filtering was carried out to maintain the quality and accuracy of the SJs and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared with an earlier version, BaRTv1.0. The accuracy of transcript level quantification, SJs and transcript start and end sites have been validated extensively using parallel technologies and analysis, including high-resolution reverse transcriptase-polymerase chain reaction and 5'-RACE. BaRTv2.18 contains 39 434 genes and 148 260 transcripts, representing the most comprehensive and resolved reference transcriptome in barley to date. It provides an important and high-quality resource for advanced transcriptomic analyses, including both transcriptional and post-transcriptional regulation, with exceptional resolution and precision.

The value of genotype-specific reference for transcriptome analyses in barley.

It is increasingly apparent that although different genotypes within a species share "core" genes, they also contain variable numbers of "specific" genes and different structures of "core" genes that are only present in a subset of individuals. Using a common reference genome may thus lead to a loss of genotype-specific information in the assembled Reference Transcript Dataset (RTD) and the generation of erroneous, incomplete or misleading transcriptomics analysis results. In this study, we assembled genotype-specific RTD (sRTD) and common reference-based RTD (cRTD) from RNA-seq data of cultivated Barke and Morex barley, respectively. Our quantitative evaluation showed that the sRTD has a significantly higher diversity of transcripts and alternative splicing events, whereas the cRTD missed 40% of transcripts present in the sRTD and it only has ∼70% accurate transcript assemblies. We found that the sRTD is more accurate for transcript quantification as well as differential expression analysis. However, gene-level quantification is less affected, which may be a reasonable compromise when a high-quality genotype-specific reference is not available.

Barley (Hordeum Vulgare) Anther and Meiocyte RNA Sequencing: Mapping Sequencing Reads and Downstream Data Analyses.

RNA sequencing (RNA-seq) data is by now the most common method to study differential gene expression. Here we present a pipeline from RNA-seq generation to analysis with examples based on our own barley anther and meiocyte transcriptome. The bioinformatics pipeline will give everyone, from a beginner to a more experienced user, the possibility to analyze their datasets and identify significantly differentially expressed genes. It also allows differential alternative splicing analysis which will become increasingly common due to the high regulatory impact on the gene expression. We describe use of the Galaxy interface for RNA-seq read quantification and the 3D RNA-seq app for the downstream data analysis.

Grain dispersal mechanism in cereals arose from a genome duplication followed by changes in spatial expression of genes involved in pollen development.

KEY MESSAGE: Grain disarticulation in wild progenitor of wheat and barley evolved through a local duplication event followed by neo-functionalization resulting from changes in location of gene expression. One of the most critical events in the process of cereal domestication was the loss of the natural mode of grain dispersal. Grain dispersal in barley is controlled by two major genes, Btr1 and Btr2, which affect the thickness of cell walls around the disarticulation zone. The barley genome also encodes Btr1-like and Btr2-like genes, which have been shown to be the ancestral copies. While Btr and Btr-like genes are non-redundant, the biological function of Btr-like genes is unknown. We explored the potential biological role of the Btr-like genes by surveying their expression profile across 212 publicly available transcriptome datasets representing diverse organs, developmental stages and stress conditions. We found that Btr1-like and Btr2-like are expressed exclusively in immature anther samples throughout Prophase I of meiosis within the meiocyte. The similar and restricted expression profile of these two genes suggests they are involved in a common biological function. Further analysis revealed 141 genes co-expressed with Btr1-like and 122 genes co-expressed with Btr2-like, with 105 genes in common, supporting Btr-like genes involvement in a shared molecular pathway. We hypothesize that the Btr-like genes play a crucial role in pollen development by facilitating the formation of the callose wall around the meiocyte or in the secretion of callase by the tapetum. Our data suggest that Btr genes retained an ancestral function in cell wall modification and gained a new role in grain dispersal due to changes in their spatial expression becoming spike specific after gene duplication.

The proteome of developing barley anthers during meiotic prophase I.

Flowering plants reproduce sexually by combining a haploid male and female gametophyte during fertilization. Male gametophytes are localized in the anthers, each containing reproductive (meiocyte) and non-reproductive tissue necessary for anther development and maturation. Meiosis, where chromosomes pair and exchange their genetic material during a process called recombination, is one of the most important and sensitive stages in breeding, ensuring genetic diversity. Most anther development studies have focused on transcript variation, but very few have been correlated with protein abundance. Taking advantage of a recently published barley anther transcriptomic (BAnTr) dataset and a newly developed sensitive mass spectrometry-based approach to analyse the barley anther proteome, we conducted high-resolution mass spectrometry analysis of barley anthers, collected at six time points and representing their development from pre-meiosis to metaphase. Each time point was carefully staged using immunocytology, providing a robust and accurate staging mirroring our previous BAnTr dataset. We identified >6100 non-redundant proteins including 82 known and putative meiotic proteins. Although the protein abundance was relatively stable throughout prophase I, we were able to quantify the dynamic variation of 336 proteins. We present the first quantitative comparative proteomics study of barley anther development during meiotic prophase I when the important process of homologous recombination is taking place.