BARLEX - the Barley Draft Genome Explorer

Non-cell-autonomous signaling associated with barley ALOG1 specifies spikelet meristem determinacy

Inflorescence architecture and crop productivity are often tightly coupled in our major cereal crops. However, the underlying genetic mechanisms controlling cereal inflorescence development remain poorly understood. Here, we identified recessive alleles of barley (Hordeum vulgare L.) HvALOG1 (Arabidopsis thaliana LSH1 and Oryza G1) that produce non-canonical extra spikelets and fused glumes abaxially to the central spikelet from the upper-mid portion until the tip of the inflorescence. Notably, we found that HvALOG1 exhibits a boundary-specific expression pattern that specifically excludes reproductive meristems, implying the involvement of previously proposed localized signaling centers for branch regulation. Importantly, during early spikelet formation, non-cell-autonomous signals associated with HvALOG1 expression may specify spikelet meristem determinacy, while boundary formation of floret organs appears to be coordinated in a cell-autonomous manner. Moreover, barley ALOG family members synergistically modulate inflorescence morphology, with HvALOG1 predominantly governing meristem maintenance and floral organ development. We further propose that spatiotemporal redundancies of expressed HvALOG members specifically in the basal inflorescence may be accountable for proper patterning of spikelet formation in mutant plants. Our research offers new perspectives on regulatory signaling roles of ALOG transcription factors during the development of reproductive meristems in cereal inflorescences.

Genome-Wide Identification and Expression Analysis of the Stearoyl-Acyl Carrier Protein Δ9 Desaturase Gene Family under Abiotic Stress in Barley

Stearoyl-acyl carrier protein (ACP) Δ9 desaturase (SAD) is a critical fatty acid dehydrogenase in plants, playing a prominent role in regulating the synthesis of unsaturated fatty acids (UFAs) and having a significant impact on plant growth and development. In this study, we conducted a comprehensive genomic analysis of the SAD family in barley (Hordeum vulgare L.), identifying 14 HvSADs with the FA_desaturase_2 domain, which were divided into four subgroups based on sequence composition and phylogenetic analysis, with members of the same subgroup possessing similar genes and motif structures. Gene replication analysis suggested that tandem and segmental duplication may be the major reasons for the expansion of the SAD family in barley. The promoters of HvSADs contained various cis-regulatory elements (CREs) related to light, abscisic acid (ABA), and methyl jasmonate (MeJA). In addition, expression analysis indicated that HvSADs exhibit multiple tissue expression patterns in barley as well as different response characteristics under three abiotic stresses: salt, drought, and cold. Briefly, this evolutionary and expression analysis of HvSADs provides insight into the biological functions of barley, supporting a comprehensive analysis of the regulatory mechanisms of oil biosynthesis and metabolism in plants under abiotic stress.

HvSL1 and HvMADS16 promote stamen identity to restrict multiple ovary formation in barley

Correct floral development is the result of a sophisticated balance of molecular cues. Floral mutants provide insight into the main genetic determinants that integrate these cues, as well as providing opportunities to assess functional variation across species. In this study, we characterize the barley (Hordeum vulgare) multiovary mutants mov2.g and mov1, and propose causative gene sequences: a C2H2 zinc-finger gene HvSL1 and a B-class gene HvMADS16, respectively. In the absence of HvSL1, florets lack stamens but exhibit functional supernumerary carpels, resulting in multiple grains per floret. Deletion of HvMADS16 in mov1 causes homeotic conversion of lodicules and stamens into bract-like organs and carpels that contain non-functional ovules. Based on developmental, genetic, and molecular data, we propose a model by which stamen specification in barley is defined by HvSL1 acting upstream of HvMADS16. The present work identifies strong conservation of stamen formation pathways with other cereals, but also reveals intriguing species-specific differences. The findings lay the foundation for a better understanding of floral architecture in Triticeae, a key target for crop improvement.

Identifying developmental QTL alleles with favorable effect on grain yield components under late-terminal drought in spring barley MAGIC population

Barley is the fourth most cultivated cereal worldwide, and drought is a major cause of its yield loss by negatively affecting its development. Hence, better understanding developmental mechanisms that control complex polygenic yield-related traits under drought is essential to uncover favorable yield regulators. This study evaluated seven above-ground yield-related traits under well-watered (WW) and late-terminal drought (TD) treatment using 534 spring barley multiparent advanced generation intercross double haploid (DH) lines. The analysis of quantitative trait loci (QTL) for WW, TD, marker by treatment interaction, and drought stress tolerance identified 69, 64, 25, and 25 loci, respectively, for seven traits from which 15 loci were common for at least three traits and 17 were shared by TD and drought stress tolerance. Evaluation of allelic effects for a QTL revealed varying effect of parental alleles. Results showed prominent QTL located on major flowering time gene Ppd-H1 with favorable effects for grain weight under TD when flowering time was not significantly affected, suggesting that this gene might be linked with increasing grain weight by ways other than timing of flowering under late-terminal drought stress. Furthermore, a desirable novel QTL allele was identified on chromosome 5H for grain number under TD nearby sucrose transporter gene HvSUT2. The findings indicated that spring barley multiparent advanced generation intercross population can provide insights to improve yield under complex condition of drought.

Genetic resources and precise gene editing for targeted improvement of barley abiotic stress tolerance

Abiotic stresses, predominately drought, heat, salinity, cold, and waterlogging, adversely affect cereal crops. They limit barley production worldwide and cause huge economic losses. In barley, functional genes under various stresses have been identified over the years and genetic improvement to stress tolerance has taken a new turn with the introduction of modern gene-editing platforms. In particular, clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) is a robust and versatile tool for precise mutation creation and trait improvement. In this review, we highlight the stress-affected regions and the corresponding economic losses among the main barley producers. We collate about 150 key genes associated with stress tolerance and combine them into a single physical map for potential breeding practices. We also overview the applications of precise base editing, prime editing, and multiplexing technologies for targeted trait modification, and discuss current challenges including high-throughput mutant genotyping and genotype dependency in genetic transformation to promote commercial breeding. The listed genes counteract key stresses such as drought, salinity, and nutrient deficiency, and the potential application of the respective gene-editing technologies will provide insight into barley improvement for climate resilience.

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.

Barley's inability to germinate after submergence depends on hypoxia-induced secondary dormancy

Global climate change has dramatically increased flooding events, which have a strong impact on crop production. Barley is one of the most important cereals and its cultivation includes a broad range of different environments. We tested the capacity to germinate of a large barley panel after a short period of submergence followed by a recovery phase. We demonstrated that sensitive barley varieties activate underwater secondary dormancy because of a lower permeability to oxygen dissolved in water. In sensitive barley accessions, secondary dormancy is removed by nitric oxide donors. Our genome wide association study results uncovered a laccase gene located in a region of significant marker-trait association that is differently regulated during grain development and plays a key role in this process. We believe that our findings will help to improve the genetics of barley thereby increasing the capacity of seeds to germinate after a short period of flooding.

MicroRNA172b-5p/trehalose-6-phosphate synthase module stimulates trehalose synthesis and microRNA172b-3p/AP2-like module accelerates flowering in barley upon drought stress

MicroRNAs (miRNAs) are major regulators of gene expression during plant development under normal and stress conditions. In this study, we analyzed the expression of 150 conserved miRNAs during drought stress applied to barley ready to flower. The dynamics of miRNAs expression was also observed after rewatering. Target messenger RNA (mRNAs) were experimentally identified for all but two analyzed miRNAs, and 41 of the targets were not reported before. Drought stress applied to barley induced accelerated flowering coordinated by a pair of two differently expressed miRNAs originating from a single precursor: hvu-miR172b-3p and hvu-miR172b-5p. Increased expression of miRNA172b-3p during drought leads to the downregulation of four APETALA2(AP2)-like genes by their mRNA cleavage. In parallel, the downregulation of the miRNA172b-5p level results in an increased level of a newly identified target, trehalose-6-phosphate synthase, a key enzyme in the trehalose biosynthesis pathway. Therefore, drought-treated plants have higher trehalose content, a known osmoprotectant, whose level is rapidly dropping after watering. In addition, trehalose-6-phosphate, an intermediate of the trehalose synthesis pathway, is known to induce flowering. The hvu-miRNA172b-5p/trehalose-6-phosphate synthase and hvu-miRNA172b-3p/AP2-like create a module leading to osmoprotection and accelerated flowering induction during drought.

Genome-wide association study in two-row spring barley landraces identifies QTL associated with plantlets root system architecture traits in well-watered and osmotic stress conditions

Water availability is undoubtedly one of the most important environmental factors affecting crop production. Drought causes a gradual deprivation of water in the soil from top to deep layers and can occur at diverse stages of plant development. Roots are the first organs that perceive water deficit in soil and their adaptive development contributes to drought adaptation. Domestication has contributed to a bottleneck in genetic diversity. Wild species or landraces represent a pool of genetic diversity that has not been exploited yet in breeding program. In this study, we used a collection of 230 two-row spring barley landraces to detect phenotypic variation in root system plasticity in response to drought and to identify new quantitative trait loci (QTL) involved in root system architecture under diverse growth conditions. For this purpose, young seedlings grown for 21 days in pouches under control and osmotic-stress conditions were phenotyped and genotyped using the barley 50k iSelect SNP array, and genome-wide association studies (GWAS) were conducted using three different GWAS methods (MLM GAPIT, FarmCPU, and BLINK) to detect genotype/phenotype associations. In total, 276 significant marker-trait associations (MTAs; p-value (FDR)< 0.05) were identified for root (14 and 12 traits under osmotic-stress and control conditions, respectively) and for three shoot traits under both conditions. In total, 52 QTL (multi-trait or identified by at least two different GWAS approaches) were investigated to identify genes representing promising candidates with a role in root development and adaptation to drought stress.

HvGST4 enhances tolerance to multiple abiotic stresses in barley: Evidence from integrated meta-analysis to functional verification

Extreme weather events have become more frequent, increasing crop yield fluctuations in many regions and thus the risk to global food security. Breeding crop cultivars with improved tolerance to a combination of abiotic stresses is an effective solution to counter the adverse impact of climate change. The ever-increasing genomic data and analytical tools provide unprecedented opportunities to mine genes with tolerance to multiple abiotic stresses through bioinformatics analysis. We undertook an integrated meta-analysis using 260 transcriptome data of barley related to drought, salt, heat, cold, and waterlogging stresses. A total of 223 shared differentially expressed genes (DEGs) were identified in response to five abiotic stresses, and significantly enriched in 'glutathione metabolism' and 'monoterpenoid biosynthesis' pathways. Using weighted gene co-expression network analysis (WGCNA), we further identified 15 hub genes (e.g., MYB, WRKY, NADH, and GST4) and selected the GST4 gene for functional validation. HvGST4 overexpression in Arabidopsis thaliana enhanced the tolerance to multiple abiotic stresses, likely through increasing the content of glutathione to scavenge reactive oxygen species and alleviate cell membrane peroxidation. Furthermore, we showed that virus-induced gene silencing (VIGS) of HvGST4 in barley leaves exacerbated cell membrane peroxidation under five abiotic stresses, reducing tolerance to multiple abiotic stress. Our study provides a new solution for identifying genes with tolerance to multiple abiotic stresses based on meta-analysis, which could contribute to breeding new varieties adapted genetically to adverse environmental conditions.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

The TCP transcription factor HvTB2 heterodimerizes with VRS5 and controls spike architecture in barley

Understanding the molecular network, including protein-protein interactions, of VRS5 provide new routes towards the identification of other key regulators of plant architecture in barley. The TCP transcriptional regulator TEOSINTE BRANCHED 1 (TB1) is a key regulator of plant architecture. In barley, an important cereal crop, HvTB1 (also referred to as VULGARE SIX-ROWED spike (VRS) 5), inhibits the outgrowth of side shoots, or tillers, and grains. Despite its key role in barley development, there is limited knowledge on the molecular network that is utilized by VRS5. In this work, we performed protein-protein interaction studies of VRS5. Our analysis shows that VRS5 potentially interacts with a diverse set of proteins, including other class II TCP's, NF-Y TF, but also chromatin remodelers. Zooming in on the interaction capacity of VRS5 with other TCP TFs shows that VRS5 preferably interacts with other class II TCP TFs in the TB1 clade. Induced mutagenesis through CRISPR-Cas of one of the putative VRS5 interactors, HvTB2 (also referred to as COMPOSITUM 1 and BRANCHED AND INDETERMINATE SPIKELET 1), resulted in plants that have lost their characteristic unbranched spike architecture. More specifically, hvtb2 mutants exhibited branches arising at the main spike, suggesting that HvTB2 acts as inhibitor of branching. Our protein-protein interaction studies of VRS5 resulted in the identification of HvTB2 as putative interactor of VRS5, another key regulator of spike architecture in barley. The study presented here provides a first step to underpin the protein-protein interactome of VRS5 and to identify other, yet unknown, key regulators of barley plant architecture.

Evaluating Variation in Germination and Growth of Landraces of Barley (Hordeum vulgare L.) Under Salinity Stress

Ongoing climate change is resulting in increasing areas of salinity affected soils, rising saline groundwater and droughts resulting in irrigation with brackish water. This leads to increased salinity stress in crops that are already grown on marginal agricultural lands, such as barley. Tolerance to salinity stress is limited in the elite barley cultivar pools, but landraces of barley hold potential sources of tolerance due to their continuous selection on marginal lands. This study analyzed 140 heritage cultivars and landrace lines of barley, including 37 Scottish Bere lines that were selected from coastal regions, to screen for tolerance to salinity stress. Tolerance to salinity stress was screened by looking at the germination speed and the early root growth during germination, and the pre-maturity biomass accumulation during early growth stages. Results showed that most lines increased germination time, and decreased shoot biomass and early root growth with greater salinity stress. Elite cultivars showed increased response to the salinity, compared to the landrace lines. Individual Bere and landrace lines showed little to no effect of increased salinity in one or more experiments, one line showed high salinity tolerance in all experiments-Bere 49 A 27 Shetland. A Genome Wide Association Screening identified a number of genomic regions associated with increased tolerance to salinity stress. Two chromosomal regions were found, one associated with shoot biomass on 5HL, and another associated with early root growth, in each of the salinities, on 3HS. Within these regions a number of promising candidate genes were identified. Further analysis of these new regions and candidate genes should be undertaken, along with field trials, to identify targets for future breeding for salinity tolerance.

Recombination landscape divergence between populations is marked by larger low-recombining regions in domesticated rye

The genomic landscape of recombination plays an essential role in evolution. Patterns of recombination are highly variable along chromosomes, between sexes, individuals, populations, and species. In many eukaryotes, recombination rates are elevated in sub-telomeric regions and drastically reduced near centromeres, resulting in large low-recombining (LR) regions. The processes of recombination are influenced by genetic factors, such as different alleles of genes involved in meiosis and chromatin structure, as well as external environmental stimuli like temperature and overall stress. In this work, we focused on the genomic landscapes of recombination in a collection of 916 rye (Secale cereale) individuals. By analysing population structure among individuals of different domestication status and geographic origin, we detected high levels of admixture, reflecting the reproductive biology of a self-incompatible, wind-pollinating grass species. We then analysed patterns of recombination in overlapping subpopulations, which revealed substantial variation in the physical size of LR regions, with a tendency for larger LR regions in domesticated subpopulations. Genome-wide association scans (GWAS) for LR region size revealed a major quantitative-trait-locus (QTL) at which, among 18 annotated genes, an ortholog of histone H4 acetyltransferase ESA1 was located. Rye individuals belonging to domesticated subpopulations showed increased synaptonemal complex length, but no difference in crossover frequency, indicating that only the recombination landscape is different. Furthermore, the genomic region harbouring rye ScESA1 showed moderate patterns of selection in domesticated subpopulations, suggesting that larger LR regions were indirectly selected for during domestication to achieve more homogeneous populations for agricultural use.

A Genome-Wide Association Study Identifying Single-Nucleotide Polymorphisms for Iron and Zinc Biofortification in a Worldwide Barley Collection

Micronutrient deficiency affects half of the world’s population, mostly in developing countries. Severe health issues such as anemia and inadequate growth in children below five years of age and pregnant women have been linked to mineral deficiencies (mostly zinc and iron). Improving the mineral content in staple crops, also known as mineral biofortification, remains the best approach to address mineral malnutrition. Barley is a staple crop in some parts of the world and is a healthy choice since it contains β-glucan, a high dietary protein. Barley mineral biofortification, especially with zinc and iron, can be beneficial since barley easily adapts to marginalized areas and requires less input than other frequently consumed cereals. In this study, we analyzed zinc and iron content in 496 barley samples. The samples were genotyped with an Illumina 50 K SNP chip. Genome-wide association studies (GWAS) identified 62 SNPs and 68 SNPs (p < 0.001) associated with iron and zinc content in grains, respectively. After a Bonferroni correction (p < 0.005), there were 12 SNPs (single-nucleotide polymorphism) associated with Zn and 6 for iron. SNP annotations revealed proteins involved in membrane transport, Zn and Fe binding, linked to nutrient remobilization in grains. These results can be used to develop biofortified barley via marker-assisted selection (MAS), which could alleviate mineral malnutrition.

Analyses of MADS-box Genes Suggest HvMADS56 to Regulate Lateral Spikelet Development in Barley

MADS-box transcription factors are crucial regulators of inflorescence and flower development in plants. Therefore, the recent interest in this family has received much attention in plant breeding programs due to their impact on plant development and inflorescence architecture. The aim of this study was to investigate the role of HvMADS-box genes in lateral spikelet development in barley (Hordeum vulgare L.). A set of 30 spike-contrasting barley lines were phenotypically and genotypically investigated under controlled conditions. We detected clear variations in the spike and spikelet development during the developmental stages among the tested lines. The lateral florets in the deficiens and semi-deficiens lines were more reduced than in two-rowed cultivars except cv. Kristina. Interestingly, cv. Kristina, int-h.43 and int-i.39 exhibited the same behavior as def.5, def.6, semi-def.1, semi-def.8 regarding development and showed reduced lateral florets size. In HOR1555, HOR7191 and HOR7041, the lateral florets continued their development, eventually setting seeds. In contrast, lateral florets in two-rowed barley stopped differentiating after the awn primordia stage giving rise to lateral floret sterility. At harvest, the lines tested showed large variation for all central and lateral spikelet-related traits. Phylogenetic analysis showed that more than half of the 108 MADS-box genes identified are highly conserved and are expressed in different barley tissues. Re-sequence analysis of a subset of these genes showed clear polymorphism in either SNPs or in/del. Variation in HvMADS56 correlated with altered lateral spikelet morphology. This suggests that HvMADS56 plays an important role in lateral spikelet development in barley.

Cold Conditioned: Discovery of Novel Alleles for Low-Temperature Tolerance in the Vavilov Barley Collection

Climate changes leading to higher summer temperatures can adversely affect cool season crops like spring barley. In the Upper Midwest region of the United States, one option for escaping this stress factor is to plant winter or facultative type cultivars in the autumn and then harvest in early summer before the onset of high-temperature stress. However, the major challenge in breeding such cultivars is incorporating sufficient winter hardiness to survive the extremely low temperatures that commonly occur in this production region. To broaden the genetic base for winter hardiness in the University of Minnesota breeding program, 2,214 accessions from the N. I. Vavilov Institute of Plant Industry (VIR) were evaluated for winter survival (WS) in St. Paul, Minnesota. From this field trial, 267 (>12%) accessions survived [designated as the VIR-low-temperature tolerant (LTT) panel] and were subsequently evaluated for WS across six northern and central Great Plains states. The VIR-LTT panel was genotyped with the Illumina 9K SNP chip, and then a genome-wide association study was performed on seven WS datasets. Twelve significant associations for WS were identified, including the previously reported frost resistance gene FR-H2 as well as several novel ones. Multi-allelic haplotype analysis revealed the most favorable alleles for WS in the VIR-LTT panel as well as another recently studied panel (CAP-LTT). Seventy-eight accessions from the VIR-LTT panel exhibited a high and consistent level of WS and select ones are being used in winter barley breeding programs in the United States and in a multiparent population.

Identification of novel QTL contributing to barley yellow mosaic resistance in wild barley (Hordeum vulgare spp. spontaneum)

BACKGROUND

Barley yellow mosaic disease (BYMD) caused by Barley yellow mosaic virus (BaYMV) and Barley mild mosaic virus (BaMMV) seriously threatens the production of winter barley. Cultivating and promoting varieties that carry disease-resistant genes is one of the most powerful ways to minimize the disease's effect on yield. However, as the BYMD virus mutates rapidly, resistance conferred by the two cloned R genes to the virus had been overcome by new virus strains. There is an urgent need for novel resistance genes in barley that convey sustainable resistance to newly emerging virus strains causing BYMD.

RESULTS

A doubled haploid (DH) population derived from a cross of SRY01 (BYMD resistant wild barley) and Gairdner (BYMD susceptible barley cultivar) was used to explore for QTL of resistance to BYMD in barley. A total of six quantitative trait loci (qRYM-1H, qRYM-2Ha, qRYM-2Hb, qRYM-3H, qRYM-5H, and qRYM-7H) related to BYMD resistance were detected, which were located on chromosomes 1H, 2H, 3H, 5H, and 7H. Both qRYM-1H and qRYM-2Ha were detected in all environments. qRYM-1H was found to be overlapped with rym7, a known R gene to the disease, whereas qRYM-2Ha is a novel QTL on chromosome 2H originated from SRY01, explaining phenotypic variation from 9.8 to 17.8%. The closely linked InDel markers for qRYM-2Ha were developed which could be used for marker-assisted selection in barley breeding. qRYM-2Hb and qRYM-3H were stable QTL for specific resistance to Yancheng and Yangzhou virus strains, respectively. qRYM-5H and qRYM-7H identified in Yangzhou were originated from Gairdner.

CONCLUSIONS

Our work is focusing on a virus disease (barley yellow mosaic) of barley. It is the first report on BYMD-resistant QTL from wild barley accessions. One novel major QTL (qRYM-2Ha) for the resistance was detected. The consistently detected new genes will potentially serve as novel sources for achieving pre-breeding barley materials with resistance to BYMD.

Genetic and molecular characterization of determinant of six-rowed spike of barley carrying vrs1.a4

Decisive role of reduced vrs1 transcript abundance in six-rowed spike of barley carrying vrs1.a4 was genetically proved and its potential causes were preliminarily analyzed. Six-rowed spike 1 (vrs1) is the major determinant of the six-rowed spike phenotype of barley (Hordeum vulgare L.). Alleles of Vrs1 have been extensively investigated. Allele vrs1.a4 in six-rowed barley is unique in that it has the same coding sequence as Vrs1.b4 in two-rowed barley. The determinant of row-type in vrs1.a4 carriers has not been experimentally identified. Here, we identified Vrs1.b4 in two-rowed accessions and vrs1.a4 in six-rowed accessions from the Qinghai-Tibet Plateau at high frequency. Genetic analyses revealed a single nuclear gene accounting for row-type alteration in these accessions. Physical mapping identified a 0.08-cM (~ 554-kb) target interval on chromosome 2H, wherein Vrs1 was the most likely candidate gene. Further analysis of Vrs1 expression in offspring of the mapping populations or different Vrs1.b4 and vrs1.a4 lines confirmed that downregulated expression of vrs1.a4 causes six-rowed spike. Regulatory sequence analysis found a single 'TA' dinucleotide deletion in vrs1.a4 carriers within a 'TA' tandem-repeat-enriched region ~ 1 kb upstream of the coding region. DNA methylation levels did not correspond to the expression difference and therefore did not affect Vrs1 expression. More evidence is needed to verify the causal link between the 'TA' deletion and the downregulated Vrs1 expression and hence the six-rowed spike phenotype.

The Jacalin-Related Lectin HvHorcH Is Involved in the Physiological Response of Barley Roots to Salt Stress

Salt stress tolerance of crop plants is a trait with increasing value for future food production. In an attempt to identify proteins that participate in the salt stress response of barley, we have used a cDNA library from salt-stressed seedling roots of the relatively salt-stress-tolerant cv. Morex for the transfection of a salt-stress-sensitive yeast strain (Saccharomyces cerevisiae YSH818 Δhog1 mutant). From the retrieved cDNA sequences conferring salt tolerance to the yeast mutant, eleven contained the coding sequence of a jacalin-related lectin (JRL) that shows homology to the previously identified JRL horcolin from barley coleoptiles that we therefore named the gene HvHorcH. The detection of HvHorcH protein in root extracellular fluid suggests a secretion under stress conditions. Furthermore, HvHorcH exhibited specificity towards mannose. Protein abundance of HvHorcH in roots of salt-sensitive or salt-tolerant barley cultivars were not trait-specific to salinity treatment, but protein levels increased in response to the treatment, particularly in the root tip. Expression of HvHorcH in Arabidopsis thaliana root tips increased salt tolerance. Hence, we conclude that this protein is involved in the adaptation of plants to salinity.

Barley Viridis-k links an evolutionarily conserved C-type ferredoxin to chlorophyll biosynthesis

Ferredoxins are single-electron carrier proteins involved in various cellular reactions. In chloroplasts, the most abundant ferredoxin accepts electrons from photosystem I and shuttles electrons via ferredoxin NADP+ oxidoreductase to generate NADPH or directly to ferredoxin dependent enzymes. In addition, plants contain other isoforms of ferredoxins. Two of these, named FdC1 and FdC2 in Arabidopsis thaliana, have C-terminal extensions and functions that are poorly understood. Here we identified disruption of the orthologous FdC2 gene in barley (Hordeum vulgare L.) mutants at the Viridis-k locus; these mutants are deficient in the aerobic cyclase reaction of chlorophyll biosynthesis. The magnesium-protoporphyrin IX monomethyl ester cyclase is one of the least characterized enzymes of the chlorophyll biosynthetic pathway and its electron donor has long been sought. Agroinfiltrations showed that the viridis-k phenotype could be complemented in vivo by Viridis-k but not by canonical ferredoxin. VirK could drive the cyclase reaction in vitro and analysis of cyclase mutants showed that in vivo accumulation of VirK is dependent on cyclase enzyme levels. The chlorophyll deficient phenotype of viridis-k mutants suggests that VirK plays an essential role in chlorophyll biosynthesis that cannot be replaced by other ferredoxins, thus assigning a specific function to this isoform of C-type ferredoxins.

Description and functional analysis of the transcriptome from malting barley

The present study aimed to establish an early model of the malting barley transcriptome, which describes the expression of genes and their ontologies, identify the period during malting with the largest dynamic shift in gene expression for future investigation, and to determine the expression patterns of all starch degrading enzyme genes relevant to the malting and brewing industry. Large dynamic increases in gene expression occurred early in malting with differential expressed genes enriched for cell wall and starch hydrolases amongst many malting related categories. Twenty-five of forty starch degrading enzyme genes were differentially expressed in the malting barley transcriptome including eleven α-amylase genes, six β-amylase genes, three α-glucosidase genes, and all five starch debranching enzyme genes. Four new or novel α-amylase genes, one β-amylase gene (Bmy3), three α-glucosidase genes, and two isoamylase genes had appreciable expression that requires further exploration into their potential relevance to the malting and brewing industry.

Genome-Wide Characterization of WRKY Transcription Factors Revealed Gene Duplication and Diversification in Populations of Wild to Domesticated Barley

The WRKY transcription factors (WRKYs) are known for their crucial roles in biotic and abiotic stress responses, and developmental and physiological processes. In barley, early studies revealed their importance, whereas their diversity at the population scale remains hardly estimated. In this study, 98 HsWRKYs and 103 HvWRKYs have been identified from the reference genome of wild and cultivated barley, respectively. The tandem duplication and segmental duplication events from the cultivated barley were observed. By taking advantage of early released exome-captured sequencing datasets in 90 wild barley accessions and 137 landraces, the diversity analysis uncovered synonymous and non-synonymous variants instead of loss-of-function mutations that had occurred at all WRKYs. For majority of WRKYs, the haplotype and nucleotide diversity both decreased in cultivated barley relative to the wild population. Five WRKYs were detected to have undergone selection, among which haplotypes of WRKY9 were enriched, correlating with the geographic collection sites. Collectively, profiting from the state-of-the-art barley genomic resources, this work represented the characterization and diversity of barley WRKY transcription factors, shedding light on future deciphering of their roles in barley domestication and adaptation.

ATP-Dependent Clp Protease Subunit C1, HvClpC1, Is a Strong Candidate Gene for Barley Variegation Mutant luteostrians as Revealed by Genetic Mapping and Genomic Re-sequencing

Implementation of next-generation sequencing in forward genetic screens greatly accelerated gene discovery in species with larger genomes, including many crop plants. In barley, extensive mutant collections are available, however, the causative mutations for many of the genes remains largely unknown. Here we demonstrate how a combination of low-resolution genetic mapping, whole-genome resequencing and comparative functional analyses provides a promising path toward candidate identification of genes involved in plastid biology and/or photosynthesis, even if genes are located in recombination poor regions of the genome. As a proof of concept, we simulated the prediction of a candidate gene for the recently cloned variegation mutant albostrians (HvAST/HvCMF7) and adopted the approach for suggesting HvClpC1 as candidate gene for the yellow-green variegation mutant luteostrians.

Genome-wide association study reveals the genetic complexity of fructan accumulation patterns in barley grain

We profiled the grain oligosaccharide content of 154 two-row spring barley genotypes and quantified 27 compounds, mainly inulin- and neoseries-type fructans, showing differential abundance. Clustering revealed two profile groups where the 'high' set contained greater amounts of sugar monomers, sucrose, and overall fructans, but lower fructosylraffinose. A genome-wide association study (GWAS) identified a significant association for the variability of two fructan types: neoseries-DP7 and inulin-DP9, which showed increased strength when applying a novel compound ratio-GWAS approach. Gene models within this region included three known fructan biosynthesis genes (fructan:fructan 1-fructosyltransferase, sucrose:sucrose 1-fructosyltransferase, and sucrose:fructan 6-fructosyltransferase). Two other genes in this region, 6(G)-fructosyltransferase and vacuolar invertase1, have not previously been linked to fructan biosynthesis and showed expression patterns distinct from those of the other three genes, including exclusive expression of 6(G)-fructosyltransferase in outer grain tissues at the storage phase. From exome capture data, several single nucleotide polymorphisms related to inulin- and neoseries-type fructan variability were identified in fructan:fructan 1-fructosyltransferase and 6(G)-fructosyltransferase genes. Co-expression analyses uncovered potential regulators of fructan biosynthesis including transcription factors. Our results provide the first scientific evidence for the distinct biosynthesis of neoseries-type fructans during barley grain maturation and reveal novel gene candidates likely to be involved in the differential biosynthesis of various types of fructan in barley.

Combined effects of a glycine-rich RNA-binding protein and a NAC transcription factor extend grain fill duration and improve malt barley agronomic performance

Two key barley genes independently control anthesis and senescence timing, enabling the manipulation of grain fill duration, grain size/plumpness, and grain protein concentration. Plant developmental processes such as flowering and senescence have direct effects on cereal yield and quality. Previous work highlighted the importance of two tightly linked genes encoding a glycine-rich RNA-binding protein (HvGR-RBP1) and a NAC transcription factor (HvNAM1), controlling barley anthesis timing, senescence, and percent grain protein. Varieties that differ in HvGR-RBP1 expression, 'Karl'(low) and 'Lewis'(high), also differ in sequence 1 KB upstream of translation start site, including an ~ 400 bp G rich insertion in the 5'-flanking region of the 'Karl' allele, which could disrupt gene expression. To improve malt quality, the (low-grain protein, delayed-senescence) 'Karl' HvNAM1 allele was introgressed into Montana germplasm. After several seasons of selection, the resulting germplasm was screened for the allelic combinations of HvGR-RBP1 and HvNAM1, finding lines combining 'Karl' alleles for both genes (-/-), lines combining 'Lewis' (functional, expressed) HvGR-RBP1 with 'Karl' HvNAM1 alleles ( ±), and lines combining 'Lewis' alleles for both genes (+ / +). Field experiments indicate that the functional ('Lewis,'  +) HvGR-RBP1 allele is associated with earlier anthesis and with slightly shorter plants, while the 'Karl' (-) HvNAM1 allele delays maturation. Genotypes carrying the ± allele combination therefore had a significantly (3 days) extended grain fill duration, leading to a higher percentage of plump kernels, slightly enhanced test weight, and lower grain protein concentration when compared to the other allele combinations. Overall, our data suggest an important function for HvGR-RBP1 in the control of barley reproductive development and set the stage for a more detailed functional analysis of this gene.

Natural variation in meiotic recombination rate shapes introgression patterns in intraspecific hybrids between wild and domesticated barley

Meiotic recombination rates vary considerably between species, populations and individuals. The genetic exchange between homologous chromosomes plays a major role in evolution by breaking linkage between advantageous and deleterious alleles in the case of introgressions. Identifying recombination rate modifiers is thus of both fundamental and practical interest to understand and utilize variation in meiotic recombination rates. We investigated recombination rate variation in a large intraspecific hybrid population (named HEB-25) derived from a cross between domesticated barley and 25 wild barley accessions. We observed quantitative variation in total crossover number with a maximum of a 1.4-fold difference between subpopulations and increased recombination rates across pericentromeric regions. The meiosis-specific α-kleisin cohesin subunit REC8 was identified as a candidate gene influencing crossover number and patterning. Furthermore, we quantified wild barley introgression patterns and revealed how local and genome-wide recombination rate variation shapes patterns of introgression. The identification of allelic variation in REC8 in combination with the observed changes in crossover patterning suggest a difference in how chromatin loops are tethered to the chromosome axis, resulting in reduced crossover suppression across pericentromeric regions. Local and genome-wide recombination rate variation is shaping patterns of introgressions and thereby directly influences the consequences of linkage drag.

BarleyVarDB: a database of barley genomic variation

Barley (Hordeum vulgare L.) is one of the first domesticated grain crops and represents the fourth most important cereal source for human and animal consumption. BarleyVarDB is a database of barley genomic variation. It can be publicly accessible through the website at http://146.118.64.11/BarleyVar. This database mainly provides three sets of information. First, there are 57 754 224 single nuclear polymorphisms (SNPs) and 3 600 663 insertions or deletions (InDels) included in BarleyVarDB, which were identified from high-coverage whole genome sequencing of 21 barley germplasm, including 8 wild barley accessions from 3 barley evolutionary original centers and 13 barley landraces from different continents. Second, it uses the latest barley genome reference and its annotation information publicly accessible, which has been achieved by the International Barley Genome Sequencing Consortium (IBSC). Third, 522 212 whole genome-wide microsatellites/simple sequence repeats (SSRs) were also included in this database, which were identified in the reference barley pseudo-molecular genome sequence. Additionally, several useful web-based applications are provided including JBrowse, BLAST and Primer3. Users can design PCR primers to asses polymorphic variants deposited in this database and use a user-friendly interface for accessing the barley reference genome. We envisage that the BarleyVarDB will benefit the barley genetic research community by providing access to all publicly available barley genomic variation information and barley reference genome as well as providing them with an ultra-high density of SNP and InDel markers for molecular breeding and identification of functional genes with important agronomic traits in barley. Database URL: http://146.118.64.11/BarleyVar.

Dissecting the Genetic Basis of Lateral and Central Spikelet Development and Grain Traits in Intermedium-Spike Barley (Hordeum vulgare Convar. Intermedium)

Barley (Hordeum vulgare L.) is one of the major grain crops worldwide and considered as a model plant for temperate cereals. One of the barley row-type groups, named intermedium-barley, was used in our previous study where we reported that other genetic loci rather than vrs1 and Int-c could play a role in lateral spikelet development and even in setting grains. To continue this work, we used phenotypic and genotypic data of 254 intermedium-spike barley accessions aimed at dissecting the genetic basis of development and grain traits of lateral and central spikelet using genome wide association (GWAS) analysis. After genotypic data filtering, 8,653 single-nucleotide polymorphism (SNPs) were used for GWAS analysis. A total of 169 significant associations were identified and we focused only on the subset of associations that exceeded the p < 10⁻⁴ threshold. Thirty-three highly significant marker-trait-associations (MTAs), represented in 28 different SNPs on all seven chromosomes for the central and/or lateral spikelet traits; such as kernel length, width, area, weight, unfilled spikelet and 1000-kernel weight, were detected. Highly significant associated markers were anchored physically using barley genome sequencing to identify candidate genes to either contain the SNPs or the closest gene to the SNP position. The results showed that 12 MTAs were specific for lateral spikelet traits, nine MTAs were specific for central spikelet traits and seven MTAs for both central and lateral traits. All together, the GWAS and candidate gene results support our hypothesis that lateral spikelet development could be regulated by loci different from those regulating central spikelet development.

Genome-wide diversity analysis of TCP transcription factors revealed cases of selection from wild to cultivated barley

Plant-specific TEOSINTE BRANCHED 1/CYCLOIDEA/PROLIFERATING CELL FACTORS 1/2 (TCP) transcription factors have known roles in inflorescence architecture. In barley, there are two family members INTERMEDIUM-C (INT-c/HvTB1-1) and COMPOSITUM 1 (COM1/HvTCP24) which are involved in the manipulation of spike architecture, whereas the participation of TCP family genes in selection from wild (Hordeum vulgare subsp. spontaneum, Hs) to cultivated barley (Hordeum vulgare subsp. vulgare, Hv) remains poorly investigated. Here, by conducting a genome-wide survey for TCP-like sequences in publicly-released datasets, 22 HsTCP and 20 HvTCP genes encoded for mature proteins were identified and assigned into two classes (I and II) based on their functional domains and the phylogenetic analysis. Each counterpart of the orthologous gene in wild and cultivated barley usually represented a similarity on the transcriptional profile across the tissues. The diversity analysis of TCPs in 90 wild barley accessions and 137 landraces with geographically-referenced passport information revealed the detectable selection at three loci including INT-c/HvTB1-1, HvPCF2, and HvPCF8. Especially, the HvPCF8 haplotypes in cultivated barley were found correlating with their geographical collection sites. There was no difference observed in either transactivation activity in yeast or subcellular localization in Nicotiana benthamiana among these haplotypes. Nevertheless, the genome-wide diversity analysis of barley TCP genes in wild and cultivated populations provided insight for future functional characterization in plant development such as spike architecture.

Targeted mutation of barley (1,3;1,4)-β-glucan synthases reveals complex relationships between the storage and cell wall polysaccharide content

Barley (Hordeum vulgare L) grain is comparatively rich in (1,3;1,4)-β-glucan, a source of fermentable dietary fibre that protects against various human health conditions. However, low grain (1,3;1,4)-β-glucan content is preferred for brewing and distilling. We took a reverse genetics approach, using CRISPR/Cas9 to generate mutations in members of the Cellulose synthase-like (Csl) gene superfamily that encode known (HvCslF6 and HvCslH1) and putative (HvCslF3 and HvCslF9) (1,3;1,4)-β-glucan synthases. Resultant mutations ranged from single amino acid (aa) substitutions to frameshift mutations causing premature stop codons, and led to specific differences in grain morphology, composition and (1,3;1,4)-β-glucan content. (1,3;1,4)-β-Glucan was absent in the grain of cslf6 knockout lines, whereas cslf9 knockout lines had similar (1,3;1,4)-β-glucan content to wild-type (WT). However, cslf9 mutants showed changes in the abundance of other cell-wall-related monosaccharides compared with WT. Thousand grain weight (TGW), grain length, width and surface area were altered in cslf6 knockouts, and to a lesser extent TGW in cslf9 knockouts. cslf3 and cslh1 mutants had no effect on grain (1,3;1,4)-β-glucan content. Our data indicate that multiple members of the CslF/H family fulfil important functions during grain development but, with the exception of HvCslF6, do not impact the abundance of (1,3;1,4)-β-glucan in mature grain.

Characterization of Genetic Diversity and Genome-Wide Association Mapping of Three Agronomic Traits in Qingke Barley (Hordeum Vulgare L.) in the Qinghai-Tibet Plateau

Barley (Hordeum vulgare L.) is one of the most important cereal crops worldwide. In the Qinghai-Tibet Plateau, six-rowed hulless (or naked) barley, called “qingke” in Chinese or “nas” in Tibetan, is produced mainly in Tibet. The complexity of the environment in the Qinghai-Tibet Plateau has provided unique opportunities for research on the breeding and adaptability of qingke barley. However, the genetic architecture of many important agronomic traits for qingke barley remains elusive. Heading date (HD), plant height (PH), and spike length (SL) are three prominent agronomic traits in barley. Here, we used genome-wide association (GWAS) mapping and GWAS with eigenvector decomposition (EigenGWAS) to detect quantitative trait loci (QTL) and selective signatures for HD, PH, and SL in a collection of 308 qingke barley accessions. The accessions were genotyped using a newly-developed, proprietary genotyping-by-sequencing (tGBS) technology, that yielded 14,970 high quality single nucleotide polymorphisms (SNPs). We found that the number of SNPs was higher in the varieties than in the landraces, which suggested that Tibetan varieties and varieties in the Tibetan area may have originated from different landraces in different areas. We have identified 62 QTLs associated with three important traits, and the observed phenotypic variation is well-explained by the identified QTLs. We mapped 114 known genes that include, but are not limited to, vernalization, and photoperiod genes. We found that 83.87% of the identified QTLs are located in the non-coding regulatory regions of annotated barley genes. Forty-eight of the QTLs are first reported here, 28 QTLs have pleotropic effects, and three QTL are located in the regions of the well-characterized genes HvVRN1, HvVRN3, and PpD-H2. EigenGWAS analysis revealed that multiple heading-date-related loci bear signatures of selection. Our results confirm that the barley panel used in this study is highly diverse, and showed a great promise for identifying the genetic basis of adaptive traits. This study should increase our understanding of complex traits in qingke barley, and should facilitate genome-assisted breeding for qingke barley improvement.

Fine mapping qGL2H, a major locus controlling grain length in barley (Hordeum vulgare L.)

A major grain length QTL on chromosome 2H was fine mapped to a 140.9 Kb region containing three genes. Increasing yield is an important target for barley breeding programs. One approach to increase yield is by enhancing individual grain weights through the regulation of grain size. Fine mapping major grain size-related quantitative trait loci is necessary for future marker-assisted selection strategies, yet studies of this nature are limited in barley. In the present study, we utilised a doubled haploid population derived from two Australian malt barley varieties, Vlamingh and Buloke, coupled with extensive genotypic and phenotypic data from three independent environments. A major grain length locus identified on chromosome 2H designated qGL2H was fine mapped to a 140.9 Kb interval. qGL2H was able to account for 25.4% of the phenotypic variation for grain length and 10.2% for grain yield. Underlying qGL2H were three high-confidence predicted genes. One of these genes encodes a MYB transcription factor and represents a promising candidate for further genetic research.

Barley's Second Spring as A Model Organism for Chloroplast Research

Barley (Hordeum vulgare) has been widely used as a model crop for studying molecular and physiological processes such as chloroplast development and photosynthesis. During the second half of the 20th century, mutants such as albostrians led to the discovery of the nuclear-encoded, plastid-localized RNA polymerase and the retrograde (chloroplast-to-nucleus) signalling communication pathway, while chlorina-f2 and xantha mutants helped to shed light on the chlorophyll biosynthetic pathway, on the light-harvesting proteins and on the organization of the photosynthetic apparatus. However, during the last 30 years, a large fraction of chloroplast research has switched to the more “user-friendly” model species Arabidopsis thaliana, the first plant species whose genome was sequenced and published at the end of 2000. Despite its many advantages, Arabidopsis has some important limitations compared to barley, including the lack of a real canopy and the absence of the proplastid-to-chloroplast developmental gradient across the leaf blade. These features, together with the availability of large collections of natural genetic diversity and mutant populations for barley, a complete genome assembly and protocols for genetic transformation and gene editing, have relaunched barley as an ideal model species for chloroplast research. In this review, we provide an update on the genomics tools now available for barley, and review the biotechnological strategies reported to increase photosynthesis efficiency in model species, which deserve to be validated in barley.

Genome-wide characterization and expression analyses of the auxin/indole-3-acetic acid (Aux/IAA) gene family in barley (Hordeum vulgare L.)

Aux/IAA genes are early auxin-responsive genes and essential for auxin signaling transduction. There is little information about Aux/IAAs in the agriculturally important cereal, barley. Using in silico method, we identified and subsequently characterized 36 Aux/IAAs from the barley genome. Based on their genomic sequences and the phylogenic relationship with Arabidopsis and rice Aux/IAA, the 36 HvIAAs were categorized into two major groups and 14 subgroups. The indication of the presence or absence of these domains for the biological functions and acting mechanisms was discussed. The cis-element distributions in HvIAA promoters suggests that the HvIAAs expressions may not only regulated by auxin (the presence of AuxREs and TGA-element) but also by other hormones and developmental and environmental cues. We then studied the HvIAAs expression in response to NAA (1-Naphthaleneacetic acid) using quantitative real-time PCR (qRT-PCR). Like the promoter analysis, only 14 HvIAAs were upregulated by NAA over two-fold at 4 h. HvIAAs were clustered into three groups based on the spatiotemporal expression data. We confirmed by qRT-PCR that most HvIAAs, especially HvIAA3, HvIAA7, HvIAA8, HvIAA18, HvIAA24 and HvIAA34, are expressed in the developing barley spike compared within seedling, suggesting their roles in regulating spike development. Taken together, our data provide a foundation for further revealing the biological function of these HvIAAs.

Supernumerary B chromosomes of Aegilops speltoides undergo precise elimination in roots early in embryo development

Not necessarily all cells of an organism contain the same genome. Some eukaryotes exhibit dramatic differences between cells of different organs, resulting from programmed elimination of chromosomes or their fragments. Here, we present a detailed analysis of programmed B chromosome elimination in plants. Using goatgrass Aegilops speltoides as a model, we demonstrate that the elimination of B chromosomes is a strictly controlled and highly efficient root-specific process. At the onset of embryo differentiation B chromosomes undergo elimination in proto-root cells. Independent of centromere activity, B chromosomes demonstrate nondisjunction of chromatids and lagging in anaphase, leading to micronucleation. Chromatin structure and DNA replication differ between micronuclei and primary nuclei and degradation of micronucleated DNA is the final step of B chromosome elimination. This process might allow root tissues to survive the detrimental expression, or overexpression of B chromosome-located root-specific genes with paralogs located on standard chromosomes.

B chromosomes are supernumerary chromosomes exhibiting dramatic differences between different organs in same species. Here, the authors show programmed B chromosome elimination in goatgrass starts at the onset of embryo differentiation by nondisjunction of chromatids, anaphase lagging, and ends with the degradation of micronucleated DNA.

Genome-Wide Analysis of the GRAS Gene Family in Barley (Hordeum vulgare L.)

The GRAS (named after first three identified proteins within this family, GAI, RGA, and SCR) family contains plant-specific genes encoding transcriptional regulators that play a key role in gibberellin (GA) signaling, which regulates plant growth and development. Even though GRAS genes have been characterized in some plant species, little research is known about the GRAS genes in barley (Hordeum vulgare L.). In this study, we observed 62 GRAS members from barley genome, which were grouped into 12 subgroups by using phylogenomic analysis together with the GRAS genes from Arabidopsis (Arabidopsis thaliana), maize (Zea mays), and rice (Oryza sativa). Chromosome localization and gene structure analysis suggested that duplication events and abundant presence of intronless genes might account for the massive expansion of GRAS gene family in barley. The analysis of RNA-seq data indicates the expression pattern of GRAS genes in various tissues at different stages in barley. Noteworthy, our qRT-PCR analysis showed the expression of 18 candidate GRAS genes abundantly in the developing inflorescence, indicating their potential roles in the barley inflorescence development and reproduction. Collectively, our evolutionary and expression analysis of GRAS family are useful for future functional characterization of GA signaling in barley and agricultural improvement.

Kmasker plants - a tool for assessing complex sequence space in plant species

Many plant genomes display high levels of repetitive sequences. The assembly of these complex genomes using short high-throughput sequence reads is still a challenging task. Underestimation or disregard of repeat complexity in these datasets can easily misguide downstream analysis. Detection of repetitive regions by k-mer counting methods has proved to be reliable. Easy-to-use applications utilizing k-mer counting are in high demand, especially in the domain of plants. We present Kmasker plants, a tool that uses k-mer count information as an assistant throughout the analytical workflow of genome data that is provided as a command-line and web-based solution. Beside its core competence to screen and mask repetitive sequences, we have integrated features that enable comparative studies between different cultivars or closely related species and methods that estimate target specificity of guide RNAs for application of site-directed mutagenesis using Cas9 endonuclease. In addition, we have set up a web service for Kmasker plants that maintains pre-computed indices for 10 of the economically most important cultivated plants. Source code for Kmasker plants has been made publically available at https://github.com/tschmutzer/kmasker. The web service is accessible at https://kmasker.ipk-gatersleben.de.

Bymovirus-induced yellow mosaic diseases in barley and wheat: viruses, genetic resistances and functional aspects

Bymovirus-induced yellow mosaic diseases seriously threaten global production of autumn-sown barley and wheat, which are two of the presently most important crops around the world. Under natural field conditions, the diseases are caused by infection of soil-borne plasmodiophorid Polymyxa graminis-transmitted bymoviruses of the genus Bymovirus of the family Potyviridae. Focusing on barley and wheat, this article summarizes the achievements on taxonomy, geography and host specificity of these disease-conferring viruses, as well as the genetics of resistance in barley, wheat and wild relatives. Moreover, based on recent progress of barley and wheat genomics, germplasm resources and large-scale sequencing, the exploration and isolation of corresponding resistant genes from wheat and barley as well as relatives, no matter what a large and complicated genome is present, are becoming feasible and are discussed. Furthermore, the foreseen advances on cloning of the resistance or susceptibility-encoding genes, which will provide the possibility to explore the functional interaction between host plants and soil-borne viral pathogens, are discussed as well as the benefits for marker-assisted resistance breeding in barley and wheat.

Salinity tolerance in barley during germination- homologs and potential genes

Salinity affects more than 6% of the world's total land area, causing massive losses in crop yield. Salinity inhibits plant growth and development through osmotic and ionic stresses; however, some plants exhibit adaptations through osmotic regulation, exclusion, and translocation of accumulated Na+ or Cl-. Currently, there are no practical, economically viable methods for managing salinity, so the best practice is to grow crops with improved tolerance. Germination is the stage in a plant's life cycle most adversely affected by salinity. Barley, the fourth most important cereal crop in the world, has outstanding salinity tolerance, relative to other cereal crops. Here, we review the genetics of salinity tolerance in barley during germination by summarizing reported quantitative trait loci (QTLs) and functional genes. The homologs of candidate genes for salinity tolerance in Arabidopsis, soybean, maize, wheat, and rice have been blasted and mapped on the barley reference genome. The genetic diversity of three reported functional gene families for salt tolerance during barley germination, namely dehydration-responsive element-binding (DREB) protein, somatic embryogenesis receptor-like kinase and aquaporin genes, is discussed. While all three gene families show great diversity in most plant species, the DREB gene family is more diverse in barley than in wheat and rice. Further to this review, a convenient method for screening for salinity tolerance at germination is needed, and the mechanisms of action of the genes involved in salt tolerance need to be identified, validated, and transferred to commercial cultivars for field production in saline soil.

Effect of epistasis and environment on flowering time in barley reveals a novel flowering-delaying QTL allele

Flowering time is a complex trait and has a key role in crop yield and adaptation to environmental stressors such as heat and drought. This study aimed to better understand the interconnected dynamics of epistasis and environment and look for novel regulators. We investigated 534 spring barley MAGIC DH lines for flowering time at various environments. Analysis of quantitative trait loci (QTLs), epistatic interactions, QTL × environment (Q×E) interactions, and epistasis × environment (E×E) interactions were performed with single SNP and haplotype approaches. In total, 18 QTLs and 2420 epistatic interactions were detected, including intervals harboring major genes such as Ppd-H1, Vrn-H1, Vrn-H3, and denso/sdw1. Epistatic interactions found in field and semi-controlled conditions were distinctive. Q×E and E×E interactions revealed that temperature influenced flowering time by triggering different interactions between known and newly detected regulators. A novel flowering-delaying QTL allele was identified on chromosome 1H (named 'HvHeading') and was shown to be engaged in epistatic and environment interactions. Results suggest that investigating epistasis, environment, and their interactions, rather than only single QTLs, is an effective approach for detecting novel regulators. We assume that barley can adapt flowering time to the environment via alternative routes within the pathway.

Perspectives on Low Temperature Tolerance and Vernalization Sensitivity in Barley: Prospects for Facultative Growth Habit

One option to achieving greater resiliency for barley production in the face of climate change is to explore the potential of winter and facultative growth habits: for both types, low temperature tolerance (LTT) and vernalization sensitivity are key traits. Sensitivity to short-day photoperiod is a desirable attribute for facultative types. In order to broaden our understanding of the genetics of these phenotypes, we mapped quantitative trait loci (QTLs) and identified candidate genes using a genome-wide association studies (GWAS) panel composed of 882 barley accessions that was genotyped with the Illumina 9K single-nucleotide polymorphism (SNP) chip. Fifteen loci including 5 known and 10 novel QTL/genes were identified for LTT-assessed as winter survival in 10 field tests and mapped using a GWAS meta-analysis. FR-H1, FR-H2, and FR-H3 were major drivers of LTT, and candidate genes were identified for FR-H3. The principal determinants of vernalization sensitivity were VRN-H1, VRN-H2, and PPD-H1. VRN-H2 deletions conferred insensitive or intermediate sensitivity to vernalization. A subset of accessions with maximum LTT were identified as a resource for allele mining and further characterization. Facultative types comprised a small portion of the GWAS panel but may be useful for developing germplasm with this growth habit.

Barley grain (1,3;1,4)-β-glucan content: effects of transcript and sequence variation in genes encoding the corresponding synthase and endohydrolase enzymes

The composition of plant cell walls is important in determining cereal end uses. Unlike other widely consumed cereal grains barley is comparatively rich in (1,3;1,4)-β-glucan, a source of dietary fibre. Previous work showed Cellulose synthase-like genes synthesise (1,3;1,4)-β-glucan in several tissues. HvCslF6 encodes a grain (1,3;1,4)-β-glucan synthase, whereas the function of HvCslF9 is unknown. Here, the relationship between mRNA levels of HvCslF6, HvCslF9, HvGlbI (1,3;1,4)-β-glucan endohydrolase, and (1,3;1,4)-β-glucan content was studied in developing grains of four barley cultivars. HvCslF6 was differentially expressed during mid (8-15 DPA) and late (38 DPA) grain development stages while HvCslF9 transcript was only clearly detected at 8-10 DPA. A peak of HvGlbI expression was detected at 15 DPA. Differences in transcript abundance across the three genes could partially explain variation in grain (1,3;1,4)-β-glucan content in these genotypes. Remarkably narrow sequence variation was found within the HvCslF6 promoter and coding sequence and does not explain variation in (1,3;1,4)-β-glucan content. Our data emphasise the genotype-dependent accumulation of (1,3;1,4)-β-glucan during barley grain development and a role for the balance between hydrolysis and synthesis in determining (1,3;1,4)-β-glucan content, and suggests that other regulatory sequences or proteins are likely to be involved in this trait in developing grain.

Barley grain (1,3;1,4)-β-glucan content: effects of transcript and sequence variation in genes encoding the corresponding synthase and endohydrolase enzymes

The composition of plant cell walls is important in determining cereal end uses. Unlike other widely consumed cereal grains barley is comparatively rich in (1,3;1,4)-β-glucan, a source of dietary fibre. Previous work showed Cellulose synthase-like genes synthesise (1,3;1,4)-β-glucan in several tissues. HvCslF6 encodes a grain (1,3;1,4)-β-glucan synthase, whereas the function of HvCslF9 is unknown. Here, the relationship between mRNA levels of HvCslF6, HvCslF9, HvGlbI (1,3;1,4)-β-glucan endohydrolase, and (1,3;1,4)-β-glucan content was studied in developing grains of four barley cultivars. HvCslF6 was differentially expressed during mid (8-15 DPA) and late (38 DPA) grain development stages while HvCslF9 transcript was only clearly detected at 8-10 DPA. A peak of HvGlbI expression was detected at 15 DPA. Differences in transcript abundance across the three genes could partially explain variation in grain (1,3;1,4)-β-glucan content in these genotypes. Remarkably narrow sequence variation was found within the HvCslF6 promoter and coding sequence and does not explain variation in (1,3;1,4)-β-glucan content. Our data emphasise the genotype-dependent accumulation of (1,3;1,4)-β-glucan during barley grain development and a role for the balance between hydrolysis and synthesis in determining (1,3;1,4)-β-glucan content, and suggests that other regulatory sequences or proteins are likely to be involved in this trait in developing grain.

Non-Invasive Phenotyping Reveals Genomic Regions Involved in Pre-Anthesis Drought Tolerance and Recovery in Spring Barley

With ongoing climate change, drought events are becoming more frequent and will affect biomass formation when occurring during pre-flowering stages. We explored growth over time under such a drought scenario, via non-invasive imaging and revealed the underlying key genetic factors in spring barley. By comparing with well-watered conditions investigated in an earlier study and including information on timing, QTL could be classified as constitutive, drought or recovery-adaptive. Drought-adaptive QTL were found in the vicinity of genes involved in dehydration tolerance such as dehydrins (Dhn4, Dhn7, Dhn8, and Dhn9) and aquaporins (e.g. HvPIP1;5, HvPIP2;7, and HvTIP2;1). The influence of phenology on biomass formation increased under drought. Accordingly, the main QTL during recovery was the region of HvPPD-H1. The most important constitutive QTL for late biomass was located in the vicinity of HvDIM, while the main locus for seedling biomass was the HvWAXY region. The disappearance of QTL marked the genetic architecture of tiller number. The most important constitutive QTL was located on 6HS in the region of 1-FEH. Stage and tolerance specific QTL might provide opportunities for genetic manipulation to stabilize biomass and tiller number under drought conditions and thereby also grain yield.

Development of a Multiparent Population for Genetic Mapping and Allele Discovery in Six-Row Barley

Germplasm collections hold valuable allelic diversity for crop improvement and genetic mapping of complex traits. To gain access to the genetic diversity within the USDA National Small Grain Collection (NSGC), we developed the Barley Recombinant Inbred Diverse Germplasm Population (BRIDG6), a six-row spring barley multiparent population (MPP) with 88 cultivated accessions crossed to a common parent (Rasmusson). The parents were randomly selected from a core subset of the NSGC that represents the genetic diversity of landrace and breeding accessions. In total, we generated 6160 F5 recombinant inbred lines (RILs), with an average of 69 and a range of 37-168 RILs per family, that were genotyped with 7773 SNPs, with an average of 3889 SNPs segregating per family. We detected 23 quantitative trait loci (QTL) associated with flowering time with five QTL found coincident with previously described flowering time genes. A major QTL was detected near the flowering time gene, HvPpd-H1 which affects photoperiod. Haplotype-based analysis of HvPpd-H1 identified private alleles to families of Asian origin conferring both positive and negative effects, providing the first observation of flowering time-related alleles private to Asian accessions. We evaluated several subsampling strategies to determine the effect of sample size on the power of QTL detection, and found that, for flowering time in barley, a sample size >50 families or 3000 individuals results in the highest power for QTL detection. This MPP will be useful for uncovering large and small effect QTL for traits of interest, and identifying and utilizing valuable alleles from the NSGC for barley improvement.

Never the Two Shall Mix: Robust Indel Markers to Ensure the Fidelity of Two Pivotal and Closely-Related Accessions of Brachypodium distachyon

Brachypodium distachyon is an established model for monocotyledonous plants. Numerous markers intended for gene discovery and population genetics have been designed. However to date, very few indel markers with larger and easily scored length polymorphism differences, that distinguish between the two morphologically similar and highly utilized B. distachyon accessions, Bd21, the reference genome accession, and Bd21-3, the transformation-optimal accession, are publically available. In this study, 22 indel markers were designed and utilized to produce length polymorphism differences of 150 bp or more, for easy discrimination between Bd21 and Bd21-3. When tested on four other B. distachyon accessions, one case of multiallelism was observed. It was also shown that the markers could be used to determine homozygosity and heterozygosity at specific loci in a Bd21 x Bd3-1 F2 population. The work done in this study allows researchers to maintain the fidelity of Bd21 and Bd21-3 stocks for both transgenic and nontransgenic studies. It also provides markers that can be utilized in conjunction with others already available for further research on population genetics, gene discovery and gene characterization, all of which are necessary for the relevance of B. distachyon as a model species.

Genome-wide association mapping in a diverse spring barley collection reveals the presence of QTL hotspots and candidate genes for root and shoot architecture traits at seedling stage

BACKGROUND

Adaptation to drought-prone environments requires robust root architecture. Genotypes with a more vigorous root system have the potential to better adapt to soils with limited moisture content. However, root architecture is complex at both, phenotypic and genetic level. Customized mapping panels in combination with efficient screenings methods can resolve the underlying genetic factors of root traits.

RESULTS

A mapping panel of 233 spring barley genotypes was evaluated for root and shoot architecture traits under non-stress and osmotic stress. A genome-wide association study elucidated 65 involved genomic regions. Among them were 34 root-specific loci, eleven hotspots with associations to up to eight traits and twelve stress-specific loci. A list of candidate genes was established based on educated guess. Selected genes were tested for associated polymorphisms. By this, 14 genes were identified as promising candidates, ten remained suggestive and 15 were rejected. The data support the important role of flowering time genes, including HvPpd-H1, HvCry2, HvCO4 and HvPRR73. Moreover, seven root-related genes, HERK2, HvARF04, HvEXPB1, PIN5, PIN7, PME5 and WOX5 are confirmed as promising candidates. For the QTL with the highest allelic effect for root thickness and plant biomass a homologue of the Arabidopsis Trx-m3 was revealed as the most promising candidate.

CONCLUSIONS

This study provides a catalogue of hotspots for seedling growth, root and stress-specific genomic regions along with candidate genes for future potential incorporation in breeding attempts for enhanced yield potential, particularly in drought-prone environments. Root architecture is under polygenic control. The co-localization of well-known major genes for barley development and flowering time with QTL hotspots highlights their importance for seedling growth. Association analysis revealed the involvement of HvPpd-H1 in the development of the root system. The co-localization of root QTL with HERK2, HvARF04, HvEXPB1, PIN5, PIN7, PME5 and WOX5 represents a starting point to explore the roles of these genes in barley. Accordingly, the genes HvHOX2, HsfA2b, HvHAK2, and Dhn9, known to be involved in abiotic stress response, were located within stress-specific QTL regions and await future validation.

Fine mapping of Rha2 in barley reveals candidate genes for resistance against cereal cyst nematode

The cereal cyst nematode resistance locus Rha2 was mapped to a 978 kbp region on the long arm of barley chromosome 2H. Three candidate genes are discussed. The cereal cyst nematode (CCN) Heterodera avenae is a soil-borne obligate parasite that can cause severe damage to cereals. This research involved fine mapping of Rha2, a CCN resistance locus on chromosome 2H of barley. Rha2 was previously mapped relative to restriction fragment length polymorphisms (RFLPs) in two mapping populations. Anchoring of flanking RFLP clone sequences to the barley genome assembly defined an interval of 5077 kbp. Genotyping-by-sequencing of resistant and susceptible materials led to the discovery of potentially useful single nucleotide polymorphisms (SNPs). Assays were designed for these SNPs and applied to mapping populations. This narrowed the region of interest to 3966 kbp. Further fine mapping was pursued by crossing and backcrossing the resistant cultivar Sloop SA to its susceptible ancestor Sloop. Evaluation of F₂ progeny confirmed that the resistance segregates as a single dominant gene. Genotyping of 9003 BC₂F₂ progeny identified recombinants. Evaluation of recombinant BC₂F₃ progeny narrowed the region of interest to 978 kbp. Two of the SNPs within this region proved to be diagnostic of CCN resistance across a wide range of barley germplasm. Fluorescence-based and gel-based assays were developed for these SNPs for use in marker-assisted selection. Within the candidate region of the reference genome, there are nine high-confidence predicted genes. Three of these, one that encodes RAR1 (a cysteine- and histidine-rich domain-containing protein), one that is predicted to encode an acetylglutamate kinase and one that is predicted to encode a tonoplast intrinsic protein, are discussed as candidate genes for CCN resistance.

A tiered approach to genome-wide association analysis for the adherence of hulls to the caryopsis of barley seeds reveals footprints of selection

BACKGROUND

Seeds of domesticated barley are grouped into two distinct types, which differ in morphology. Caryopses covered by adaxial (palea) and abaxial (lemma) hulls that tightly adhere to the pericarp at maturity give rise to hulled seeds whereas caryopses without adhering hulls give rise to naked seeds. The naked caryopsis character is an essential trait regarding the end use of barley.

RESULTS

To uncover the genetic basis of the trait, a genome-wide association study (GWAS) has been performed in a panel comprising 222 2-rowed and 303 6-rowed spring barley landrace accessions. In addition to the well-described Nud locus on chromosome 7H, three novel loci showed strong associations with the trait: the first locus on 2H was specifically detected in 6-rowed accessions, the second locus on 3H was found in 2-rowed accessions from Eurasia and the third locus on 6H was revealed in 6-rowed accessions from Ethiopia. PCR analysis of naked accessions also confirmed the absence of a 17 kb region harboring the Nud gene on chromosome 7H for all but one naked accession. The latter was characterized by a slightly variant phenotype of the caryopsis.

CONCLUSION

Our findings provide evidence of the pervasiveness of the 17 kb deletion in spring barley from different geographic regions and at the same time reveal genomic footprints of selection in naked barley, which follow both geographic and morphological patterns.