Search

ABSTRACT Faba bean (Vicia faba L.) is an important pulse crop with a wide range of agroecological adaptations. The development of genomic tools and a comprehensive catalog of extant genetic diversity are crucial for developing improved faba bean cultivars. The lack of a cost‐effective genotyping platform limits the characterization of large germplasm collections, understanding of genetic diversity across populations, and implementing breeder's tools like genomic selection. Genotyping‐by‐sequencing (GBS) offers high‐resolution genotyping for both model and crop plant species, even without a reference genome sequence. The genome fragments targeted by GBS depend substantially on the restriction enzyme (RE) used for the complexity reduction step. Species with complex genomic architecture require optimization of GBS with proper RE to realize the full potential of GBS. Here, we evaluated various REs in the GBS method and identified that the combination of ApeKI/MseI proved to be the most appropriate for faba bean based on the best library quality, a high number of genomic loci spread across chromosomes, and high enrichment loci associated with the gene space. With the new optimized protocol, we constructed a genetic map using a recombinant inbred line (RIL) population and identified a QTL for seed hilum color on Chromosome 1. In addition, we also genotyped a diversity panel and performed a genome‐wide association studies (GWAS) for important agronomic traits, including plant height (PH), flowering time (FT), and number of pods per plant (PPP). We identified six SNP markers significantly associated with these traits and listed potential candidate genes. The optimized faba bean‐specific GBS procedure will facilitate access to the untapped genetic diversity for genetic research and breeding and may facilitate functional genomics.

Abstract 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.

Abstract 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.

This work is an update and extension of the previously published article “Ultralong Oxford Nanopore Reads Enable the Development of a Reference-Grade Perennial Ryegrass Genome Assembly” by Frei et al. The published genome assembly of the doubled haploid perennial ryegrass (Lolium perenne L.) genotype Kyuss (Kyuss v1.0) marked a milestone for forage grass research and breeding. However, order and orientation errors may exist in the pseudo-chromosomes of Kyuss, since barley (Hordeum vulgare L.), which diverged 30 million years ago from perennial ryegrass, was used as the reference to scaffold Kyuss. To correct for structural errors possibly present in the published Kyuss assembly, we de novo assembled the genome again and generated 50-fold coverage high-throughput chromosome conformation capture (Hi-C) data to assist pseudo-chromosome construction. The resulting new chromosome-level assembly Kyuss v2.0 showed improved quality with high contiguity (contig N50 = 120 Mb), high completeness (total BUSCO score = 99%), high base-level accuracy (QV = 50), and correct pseudo-chromosome structure (validated by Hi-C contact map). This new assembly will serve as a better reference genome for Lolium spp. and greatly benefit the forage and turf grass research community.

Measurement(s) genome sequence • wheat stripe rust disease resistance • grain yield trait Technology Type(s) genotyping-by-sequencing • whole genome shotgun sequencing • Manual scoring in natural and artificially infected fields • Combine harvester Factor Type(s) genotype • location Sample Characteristic - Organism Triticum aestivum Sample Characteristic - Environment agricultural field Sample Characteristic - Location Germany

Genome-wide prediction is a powerful tool in breeding. Initial results suggest that genome-wide approaches are also promising for enhancing the use of the genebank material: predicting the performance of plant genetic resources can unlock their hidden potential and fill the information gap in genebanks across the world and, hence, underpin prebreeding programs. As a proof of concept, we evaluated the power of across-genebank prediction for extensive germplasm collections relying on historical data on flowering/heading date, plant height, and thousand kernel weight of 9,344 barley (Hordeum vulgare L.) plant genetic resources from the German Federal Ex situ Genebank for Agricultural and Horticultural Crops (IPK) and of 1,089 accessions from the International Center for Agriculture Research in the Dry Areas (ICARDA) genebank. Based on prediction abilities for each trait, three scenarios for predictive characterization were compared: 1) a benchmark scenario, where test and training sets only contain ICARDA accessions, 2) across-genebank predictions using IPK as training and ICARDA as test set, and 3) integrated genebank predictions that include IPK with 30% of ICARDA accessions as a training set to predict the rest of ICARDA accessions. Within the population of ICARDA accessions, prediction abilities were low to moderate, which was presumably caused by a limited number of accessions used to train the model. Interestingly, ICARDA prediction abilities were boosted up to ninefold by using training sets composed of IPK plus 30% of ICARDA accessions. Pervasive genotype × environment interactions (GEIs) can become a potential obstacle to train robust genome-wide prediction models across genebanks. This suggests that the potential adverse effect of GEI on prediction ability was counterbalanced by the augmented training set with certain connectivity to the test set. Therefore, across-genebank predictions hold the promise to improve the curation of the world’s genebank collections and contribute significantly to the long-term development of traditional genebanks toward biodigital resource centers.

Leaf rust is an economically significant disease of barley. Here the authors describe cloning of the barley Rph3 leaf rust resistance gene and reveal it encodes a predicted transmembrane protein that is expressed upon infection by Rph3-avirulent Puccinia hordei isolates.

Aegilops sharonensis is a wild diploid relative of wheat. Here, the authors assemble the genome of Ae. sharonensis and use the assembly as an aid to clone the Ae. sharonensis-derived stem rust resistance gene Sr62 in the allohexaploid genome of wheat.

Consideration of the properties of awns is important for the phenotypic description of grain crops. Awns have a number of important functions in grasses, including assimilation, mechanical protection, and seed dispersal and burial. An important feature of the awn is the presence or absence of barbs—tiny hook-like single-celled trichomes on the outer awn surface that can be visualized using microscopic imaging. There are, however, no suitable software tools for the automated analysis of these small, semi-transparent structures in a high-throughput manner. Furthermore, automated analysis of barbs using conventional methods of pattern detection and segmentation is hampered by high variability of their optical appearance including size, shape, and surface density. In this work, we present a software tool for automated detection and phenotyping of barbs in microscopic images of awns, which is based on a dedicated deep learning model (BarbNet). Our experimental results show that BarbNet is capable of detecting barb structures in different awn phenotypes with an average accuracy of 90%. Furthermore, we demonstrate that phenotypic traits derived from BarbNet-segmented images enable a quite robust categorization of 4 contrasting awn phenotypes with an accuracy of >85%. Based on the promising results of this work, we see that the proposed model has potential applications in the automation of barley awns sorting for plant developmental analysis.

Abstract Wheat has one of the largest and most repetitive genomes among major crop plants, containing over 85% transposable elements (TEs). TEs populate genomes much in the way that individuals populate ecosystems, diversifying into different lineages, sub‐families and sub‐populations. The recent availability of high‐quality, chromosome‐scale genome sequences from ten wheat lines enables a detailed analysis how TEs evolved in allohexaploid wheat, its diploids progenitors, and in various chromosomal haplotype segments. LTR retrotransposon families evolved into distinct sub‐populations and sub‐families that were active in waves lasting several hundred thousand years. Furthermore, It is shown that different retrotransposon sub‐families were active in the three wheat sub‐genomes, making them useful markers to study and date polyploidization events and chromosomal rearrangements. Additionally, haplotype‐specific TE sub‐families are used to characterize chromosomal introgressions in different wheat lines. Additionally, populations of non‐autonomous TEs co‐evolved over millions of years with their autonomous partners, leading to complex systems with multiple types of autonomous, semi‐autonomous and non‐autonomous elements. Phylogenetic and TE population analyses revealed the relationships between non‐autonomous elements and their mobilizing autonomous partners. TE population analysis provided insights into genome evolution of allohexaploid wheat and genetic diversity of species, and may have implication for future crop breeding.

Genetic pathogen control is an economical and sustainable alternative to the use of chemicals. In order to breed resistant varieties, information about potentially unused genetic resistance mechanisms is of high value. We phenotyped 8,316 genotypes of the winter wheat collection of the German Federal ex situ gene bank for Agricultural and Horticultural Crops, Germany, for resistance to powdery mildew (PM), Blumeria graminis f. sp. tritici, one of the most important biotrophic pathogens in wheat. To achieve this, we used a semi-automatic phenotyping facility to perform high-throughput detached leaf assays. This data set, combined with genotyping-by-sequencing (GBS) marker data, was used to perform a genome-wide association study (GWAS). Alleles of significantly associated markers were compared with SNP profiles of 171 widely grown wheat varieties in Germany to identify currently unexploited resistance conferring genes. We also used the Chinese Spring reference genome annotation and various domain prediction algorithms to perform a domain enrichment analysis and produced a list of candidate genes for further investigation. We identified 51 significantly associated regions. In most of these, the susceptible allele was fixed in the tested commonly grown wheat varieties. Eleven of these were located on chromosomes for which no resistance conferring genes have been previously reported. In addition to enrichment of leucine-rich repeats (LRR), we saw enrichment of several domain types so far not reported as relevant to PM resistance, thus, indicating potentially novel candidate genes for the disease resistance research and prebreeding in wheat.

Grasses have diverse inflorescence morphologies, but the underlying genetic mechanisms are unclear. Here, the authors report a TCP transcription factor COM1 affects cell growth through regulation of cell wall properties and promotes branch formation in non-Triticeae grasses but branch inhibition in barley (Triticeae).

Barley (Hordeum vulgare) is one of the most important crops worldwide and is also considered a research model for the large-genome small grain temperate cereals. Despite genomic resources improving all the time, they are limited for the cv. Golden Promise, the most efficient genotype for genetic transformation. We have developed a barley cv. Golden Promise reference assembly integrating Illumina paired-end reads, long mate-pair reads, Dovetail Chicago in vitro proximity ligation libraries and chromosome conformation capture sequencing (Hi-C) libraries into a contiguous reference assembly. The assembled genome of 7 chromosomes and 4.13Gb in size, has a super-scaffold N50 after Chicago libraries of 4.14Mb and contains only 2.2% gaps. Using BUSCO (benchmarking universal single copy orthologous genes) as evaluation the genome assembly contains 95.2% of complete and single copy genes from the plant database. A high-quality Golden Promise reference assembly will be useful and utilized by the whole barley research community but will prove particularly useful for CRISPR-Cas9 experiments.

Abstract Chromosome-scale genome sequence assemblies underpin pan-genomic studies. Recent genome assembly efforts in the large-genome Triticeae crops wheat and barley have relied on the commercial closed-source assembly algorithm DeNovoMagic. We present TRITEX, an open-source computational workflow that combines paired-end, mate-pair, 10X Genomics linked-read with chromosome conformation capture sequencing data to construct sequence scaffolds with megabase-scale contiguity ordered into chromosomal pseudomolecules. We evaluate the performance of TRITEX on publicly available sequence data of tetraploid wild emmer and hexaploid bread wheat, and construct an improved annotated reference genome sequence assembly of the barley cultivar Morex as a community resource.

Gene pairs resulting from whole genome duplication (WGD), so-called ohnologous genes, are retained if at least one member of the pair undergoes neo- or sub-functionalization. Phylogenetic analyses of the ohnologous genes ALBOSTRIANS (HvAST/HvCMF7) and ALBOSTRIANS-LIKE (HvASL/HvCMF3) of barley (Hordeum vulgare) revealed them as members of a subfamily of genes coding for CCT motif (CONSTANS, CONSTANS-LIKE and TIMING OF CAB1) proteins characterized by a single CCT domain and a putative N-terminal chloroplast transit peptide. Recently, we showed that HvCMF7 is needed for chloroplast ribosome biogenesis. Here we demonstrate that mutations in HvCMF3 lead to seedlings delayed in development. They exhibit a yellowish/light green – xantha – phenotype and successively develop pale green leaves. Compared to wild type, plastids of mutant seedlings show a decreased PSII efficiency, impaired processing and reduced amounts of ribosomal RNAs; they contain less thylakoids and grana with a higher number of more loosely stacked thylakoid membranes. Site-directed mutagenesis of HvCMF3 identified a previously unknown functional domain, which is highly conserved within this subfamily of CCT domain containing proteins. HvCMF3:GFP fusion constructs were localized to plastids and nucleus. Hvcmf3Hvcmf7 double mutants exhibited a xantha-albino or albino phenotype depending on the strength of molecular lesion of the HvCMF7 allele. The chloroplast ribosome deficiency is discussed as the primary observed defect of the Hvcmf3 mutants. Based on our observations, the genes HvCMF3 and HvCMF7 have similar but not identical functions in chloroplast development of barley supporting our hypothesis of neo-/sub-functionalization between both ohnologous genes.

Improvement of grain yield is the ultimate goal for wheat breeding under water-limited environments. In the present study, a high-density linkage map was developed by using genotyping-by-sequencing (GBS) of a recombinant inbred line (RIL) population derived from the cross between Iranian landrace #49 and cultivar Yecora Rojo. The population was evaluated in three locations in Iran during two years under irrigated and water deficit conditions for the agronomic traits grain yield (GY), plant height (PH), spike number per square meter (SM), 1000 kernel weight (TKW), grain number per spike (GNS), spike length (SL), biomass (BIO) and harvest index (HI). A linkage map was constructed using 5831 SNPs assigned to 21 chromosomes, spanning 3642.14 cM of the hexaploid wheat genome with an average marker density of 0.62 (markers/cM). In total, 85 QTLs were identified on 19 chromosomes (all except 5D and 6D) explaining 6.06–19.25% of the traits phenotypic variance. We could identify 20 novel QTLs explaining 8.87–19.18% of phenotypic variance on chromosomes 1A, 1B, 1D, 2B, 3A, 3B, 6A, 6B and 7A. For 35 out of 85 mapped QTLs functionally annotated genes were identified which could be related to a potential role in drought stress.

AbstractBarley (Hordeum vulgare

The Arabidopsis gene Chloroplast Import Apparatus 2 (CIA2) encodes a transcription factor that positively affects the activity of nuclear genes for chloroplast ribosomal proteins and chloroplast protein import machineries. CIA2-like (CIL) is the paralogous gene of CIA2. We generated a cil mutant by site-directed mutagenesis and compared it with cia2 and cia2cil double mutant. Phenotype of the cil mutant did not differ from the wild type under our growth conditions, except faster growth and earlier time to flowering. Compared to cia2, the cia2cil mutant showed more impaired chloroplast functions and reduced amounts of plastid ribosomal RNAs. In silico analyses predict for CIA2 and CIL a C-terminal CCT domain and an N-terminal chloroplast transit peptide (cTP). Chloroplast (and potentially nuclear) localization was previously shown for HvCMF3 and HvCMF7, the homologs of CIA2 and CIL in barley. We observed nuclear localization of CIL after transient expression in Arabidopsis protoplasts. Surprisingly, transformation of cia2 with HvCMF3, HvCMF7, or with a truncated CIA2 lacking the predicted cTP could partially rescue the pale-green phenotype of cia2. These data are discussed with respect to potentially overlapping functions between CIA2, CIL, and their barley homologs and to the function of the putative cTPs of CIA2 and CIL.

The origin of Tibetan barley (qingke) has been a controversial issue for many years. Here, the authors conduct population genomics study to support that qingke is derived from eastern domesticated barley instead of Tibetan wild barley and suggest southern Tibetan Plateau as its introduction route.

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.

Genebanks harbor a large treasure trove of untapped plant genetic diversity. A growing world population and a changing climate require an increase in the production and development of stress resistant plant cultivars while decreasing the acreage. These requirements for improved plant cultivars can be supported by the broader exploitation of plant genetic resources (PGR) as inputs for genomics-assisted breeding. To support this process we have developed BRIDGE, a data warehouse and exploratory data analysis tool for genebank genomics of barley (Hordeum vulgare L.). Using efficient technologies for data storage, data transfer and web development, we facilitate access to digital genebank resources of barley by prioritizing the interactive and visual analysis of integrated genotypic and phenotypic data. The underlying data resulted from a barley genebank genomics study cataloging sequence and morphological data of 22,626 barley accessions, mainly from the German Federal ex situ genebank. BRIDGE consists of interactively coupled modules to visualize integrated, curated and quality checked data, such as variation data, results of dimensionality reduction and genome wide association studies (GWAS), phenotyping results, passport data as well as the geographic distribution of germplasm samples. The core component is a manager for custom collections of germplasm. A search module to find and select germplasm by passport and phenotypic attributes is included as well as modules to export genotypic data in gzip-compressed variant call format (VCF) files and phenotypic data in MIAPPE-compliant ISA-Tab files. BRIDGE is accessible at the following URL: https://bridge.ipk-gatersleben.de.

Powdery mildew caused by Blumeria graminis f. sp. hordei (Bgh) is one of the main foliar diseases in barley (Hordeum vulgare L.; Hv). Naturally occurring resistance genes used in barley breeding are a cost effective and environmentally sustainable strategy to minimize the impact of pathogens, however, the primary gene pool of H. vulgare contains limited diversity owing to recent domestication bottlenecks. To ensure durable resistance against this pathogen, more genes are required that could be unraveled by investigation of secondary barley gene-pool. A large set of Hordeum bulbosum (Hb) introgression lines (ILs) harboring a diverse set of desirable resistance traits have been developed and are being routinely used as source of novel diversity in gene mapping studies. Nevertheless, this strategy is often compromised by a lack of recombination between the introgressed fragment and the orthologous chromosome of the barley genome. In this study, we fine-mapped a Hb gene conferring resistance to barley powdery mildew. The initial genotyping of two Hb ILs mapping populations with differently sized 2HS introgressions revealed severely reduced interspecific recombination in the region of the introgressed segment, preventing precise localization of the gene. To overcome this difficulty, we developed an alternative strategy, exploiting intraspecific recombination by crossing two Hv/Hb ILs with collinear Hb introgressions, one of which carries a powdery mildew resistance gene, while the other doesn’t. The intraspecific recombination rate in the Hb-introgressed fragment of 2HS was approximately 20 times higher than it was in the initial simple ILs mapping populations. Using high-throughput genotyping-by-sequencing (GBS), we allocated the resistance gene to a 1.4 Mb interval, based on an estimate using the Hv genome as reference, in populations of only 103 and 146 individuals, respectively, similar to what is expected at this locus in barley. The most likely candidate resistance gene within this interval is part of the coiled-coil nucleotide-binding-site leucine-rich-repeat (CC-NBS-LLR) gene family, which is over-represented among genes conferring strong dominant resistance to pathogens. The reported strategy can be applied as a general strategic approach for identifying genes underlying traits of interest in crop wild relatives.

Genebank genomics promises to unlock valuable diversity for plant breeding but first, one key question is which marker system is most suitable to fingerprint entire genebank collections. Using wheat as model species, we tested for the presence of an ascertainment bias and investigated its impact on estimates of genetic diversity and prediction ability obtained using three marker platforms: simple sequence repeat (SSR), genotyping-by-sequencing (GBS), and array-based SNP markers. We used a panel of 378 winter wheat genotypes including 190 elite lines and 188 plant genetic resources (PGR), which were phenotyped in multi-environmental trials for grain yield and plant height. We observed an ascertainment bias for the array-based SNP markers, which led to an underestimation of the molecular diversity within the population of PGR. In contrast, the marker system played only a minor role for the overall picture of the population structure and precision of genome-wide predictions. Interestingly, we found that rare markers contributed substantially to the prediction ability. This combined with the expectation that valuable novel diversity is most likely rare suggests that markers with minor allele frequency deserve careful consideration in the design of a pre-breeding program.

Abstract Understanding how crop plants evolved from their wild relatives and spread around the world can inform about the origins of agriculture. Here, we review how the rapid development of genomic resources and tools has made it possible to conduct genetic mapping and population genetic studies to unravel the molecular underpinnings of domestication and crop evolution in diverse crop species. We propose three future avenues for the study of crop evolution: establishment of high-quality reference genomes for crops and their wild relatives; genomic characterization of germplasm collections; and the adoption of novel methodologies such as archaeogenetics, epigenomics, and genome editing.

Abstract Background The large and highly repetitive genomes of the cultivated species Hordeum vulgare (barley), Triticum aestivum (wheat), and Secale cereale (rye) belonging to the Triticeae tribe of grasses appear to be particularly rich in gene-like sequences including partial duplicates. Most of them have been classified as putative pseudogenes. In this study we employ transient and stable gene silencing- and over-expression systems in barley to study the function of HvARM1 (for H. vulgare Armadillo 1), a partial gene duplicate of the U-box/armadillo-repeat E3 ligase HvPUB15 (for H. vulgare Plant U-Box 15). Results The partial ARM1 gene is derived from a gene-duplication event in a common ancestor of the Triticeae and contributes to quantitative host as well as nonhost resistance to the biotrophic powdery mildew fungus Blumeria graminis. In barley, allelic variants of HvARM1 but not of HvPUB15 are significantly associated with levels of powdery mildew infection. Both HvPUB15 and HvARM1 proteins interact in yeast and plant cells with the susceptibility-related, plastid-localized barley homologs of THF1 (for Thylakoid formation 1) and of ClpS1 (for Clp-protease adaptor S1) of Arabidopsis thaliana. A genome-wide scan for partial gene duplicates reveals further events in barley resulting in stress-regulated, potentially neo-functionalized, genes. Conclusion The results suggest neo-functionalization of the partial gene copy HvARM1 increases resistance against powdery mildew infection. It further links plastid function with susceptibility to biotrophic pathogen attack. These findings shed new light on a novel mechanism to employ partial duplication of protein-protein interaction domains to facilitate the expansion of immune signaling networks.

Abstract Background While transposable elements (TEs) comprise the bulk of plant genomic DNA, how they contribute to genome structure and organization is still poorly understood. Especially in large genomes where TEs make the majority of genomic DNA, it is still unclear whether TEs target specific chromosomal regions or whether they simply accumulate where they are best tolerated. Results Here, we present an analysis of the repetitive fraction of the 5100 Mb barley genome, the largest angiosperm genome to have a near-complete sequence assembly. Genes make only about 2% of the genome, while over 80% is derived from TEs. The TE fraction is composed of at least 350 different families. However, 50% of the genome is comprised of only 15 high-copy TE families, while all other TE families are present in moderate or low copy numbers. We found that the barley genome is highly compartmentalized with different types of TEs occupying different chromosomal “niches”, such as distal, interstitial, or proximal regions of chromosome arms. Furthermore, gene space represents its own distinct genomic compartment that is enriched in small non-autonomous DNA transposons, suggesting that these TEs specifically target promoters and downstream regions. Furthermore, their presence in gene promoters is associated with decreased methylation levels. Conclusions Our data show that TEs are major determinants of overall chromosome structure. We hypothesize that many of the the various chromosomal distribution patterns are the result of TE families targeting specific niches, rather than them accumulating where they have the least deleterious effects.

Abstract Background Short-culm mutants have been widely used in breeding programs to increase lodging resistance. In barley (Hordeum vulgare L.), several hundreds of short-culm mutants have been isolated over the years. The objective of the present study was to identify the Brachytic1 (Brh1) semi-dwarfing gene and to test its effect on yield and malting quality. Results Double-haploid lines generated through a cross between a brh1.a mutant and the European elite malting cultivar Quench, showed good malting quality but a decrease in yield. Especially the activities of the starch degrading enzymes β-amylase and free limit dextrinase were high. A syntenic approach comparing markers in barley to those in rice (Oryza sativa L.), sorghum (Sorghum bicolor Moench) and brachypodium (Brachypodium distachyon P. Beauv) helped us to identify Brh1 as an orthologue of rice D1 encoding the Gα subunit of a heterotrimeric G protein. We demonstrated that Brh1 is allelic to Ari-m. Sixteen different mutant alleles were described at the DNA level. Conclusions Mutants in the Brh1 locus are deficient in the Gα subunit of a heterotrimeric G protein, which shows that heterotrimeric G proteins are important regulators of culm length in barley. Mutant alleles do not have any major negative effects on malting quality.

Markers linked to agronomic traits are of the prerequisite for molecular breeding. Genotyping-by-sequencing (GBS) data enables to detect small polymorphisms including single nucleotide polymorphisms (SNPs) and short insertions or deletions (InDels) that can be used, for instance, for marker-assisted selection, population genetics, and genome-wide association studies (GWAS). Here, we aim at detecting large chromosomal modifications in barley and wheat based on GBS data. These modifications could be duplications, deletions, substitutions including introgressions as well as alterations of DNA methylation. We demonstrate that GBS coverage analysis is capable to detect Hordeum vulgare/Hordeum bulbosum introgression lines. Furthermore, we identify large chromosomal modifications in barley and wheat collections. Hence, large chromosomal modifications, including introgressions and copy number variations (CNV), can be detected easily and can be used as markers in research and breeding without additional wet-lab experiments.

We compared the performance of two commonly used genotyping platforms, genotyping-by-sequencing (GBS) and single nucleotide polymorphism-arrays (SNP), to investigate the extent and pattern of genetic variation within a collection of 1,000 diverse barley genotypes selected from the German Federal ex situ GenBank hosted at IPK Gatersleben. Each platform revealed equivalent numbers of robust bi-allelic SNPs (39,733 and 37,930 SNPs for the 50K SNP-array and GBS datasets respectively). A small overlap of 464 SNPs was common to both platforms, indicating that the methodologies we used selectively access informative polymorphism in different portions of the barley genome. Approximately half of the GBS dataset was comprised of SNPs with minor allele frequencies (MAFs) below 1%, illustrating the power of GBS to detect rare alleles in diverse germplasm collections. While desired for certain applications, the highly robust calling of alleles at the same SNPs across multiple populations is an advantage of the SNP-array, allowing direct comparisons of data from related or unrelated studies. Overall MAFs and diversity statistics (π) were higher for the SNP-array data, potentially reflecting the conscious removal of markers with a low MAF in the ascertainment population. A comparison of similarity matrices revealed a positive correlation between both approaches, supporting the validity of using either for entire GenBank characterization. To explore the potential of each dataset for focused genetic analyses we explored the outcomes of their use in genome-wide association scans for row type, growth habit and non-adhering hull, and discriminant analysis of principal components for the drivers of sub-population differentiation. Interpretation of the results from both types of analysis yielded broadly similar conclusions indicating that choice of platform used for such analyses should be determined by the research question being asked, group preferences and their capabilities to extract and interpret the different types of output data easily and quickly. Access to the requisite infrastructure for running, processing, analyzing, querying, storing, and displaying either datatype is an additional consideration. Our investigations reveal that for barley the cost per genotyping assay is less for SNP-arrays than GBS, which translates to a cost per informative datapoint being significantly lower for the SNP-array.

Powdery mildew caused by Blumeria graminis f. sp. hordei is a foliar disease with highly negative impact on yield and grain quality in barley. Thus, breeding for powdery mildew resistance is an important goal and requires constantly the discovery of new sources of natural resistance. Here, we report the high resolution genetic and physical mapping of a dominant race-specific powdery mildew resistance locus, originating from an Ethiopian spring barley accession ‘HOR2573,’ conferring resistance to several modern mildew isolates. High-resolution genetic mapping narrowed down the interval containing the resistance locus to a physical span of 850 kb. Four candidate genes with homology to known disease resistance gene families were identified. The mapped resistance locus coincides with a previously reported resistance locus from Hordeum laevigatum, suggesting allelism at the same locus in two different barley lines. Therefore, we named the newly mapped resistance locus from HOR2573 as MlLa-H. The reported co-segregating and flanking markers may provide new tools for marker-assisted selection of this resistance locus in barley breeding.

Chromosome conformation capture sequencing (Hi-C) is a powerful method to comprehensively interrogate the three-dimensional positioning of chromatin in the nucleus. The development of Hi-C can be traced back to successive increases in the resolution and throughput of chromosome conformation capture (3C) (Dekker et al., 2002). The basic workflow of 3C consists of (i) fixation of intact chromatin, usually by formaldehyde, (ii) cutting the fixed chromatin with a restriction enzyme, (iii) religation of sticky ends under diluted conditions to favor ligations between cross-linked fragments or those between random fragments and (iv) quantifying the number of ligations events between pairs of genomic loci (de Wit and de Laat, 2012). In the original 3C protocol, ligation frequency was measured by amplification of selected ligation junctions corresponding to a small number of genomic loci (‘one versus one’) through semi-quantitative PCR (Dekker et al., 2002). The chromosome conformation capture-on-chip (4C) and chromosome conformation capture carbon copy (5C) technologies then extended 3C to count ligation events in a ‘one versus many’ or ‘many versus many’ manner, respectively. Hi-C (Lieberman-Aiden et al., 2009) finally combined 3C with next-generation sequencing (Metzker, 2010). Here, before religation sticky ends are filled in with biotin-labeled nucleotide analogs to enrich for fragments with a ligation junction in a later step. The Hi-C libraries are then subjected to high-throughput sequencing and the resultant reads mapped to a reference genome, allowing the determination of contact probabilities in a ‘many versus many’ way with a resolution that is limited only by the distribution of restriction sites and the read depth. The first application of Hi-C was the elucidation of global chromatin folding principles in the human genome (Lieberman-Aiden et al., 2009). Similar efforts have since been carried out in other eukaryotic model species such as yeast (Duan et al., 2010), Drosophila (Sexton et al., 2012) and Arabidopsis (Grob et al., 2014; Wang et al., 2015; Liu et al., 2016). Other uses of Hi-C include the study of chromatin looping at high-resolution (Rao et al., 2014; Liu et al., 2016), of chromatin reorganization along the cell cycle (Naumova et al., 2013) and of differences in chromatin organization in mutant individuals (Feng et al., 2014). The tethered conformation capture protocol (TCC) (Kalhor et al., 2011) described here is a variant of the original Hi-C method (Lieberman-Aiden et al., 2009) and was adapted to Triticeae.

L., a wild grass and close relative of cultivated barley ( L.), gained importance in plant breeding as inducer of haploid plants in crosses with barley and also as a genetic resource for introgression of disease resistance/tolerance genes into cultivated barley. Genetic mapping of genes introgressed from is a prerequisite for their efficient utilization in barley breeding, but often hindered due to repressed recombination. The mechanism underlying the reduced frequency or lack of meiotic recombination between . and . chromatin in introgressed segments is not understood. It may be explained by lack of genome collinearity or other structural differences between both genomes. In the present study, two F mapping populations of were analyzed by genotyping-by-sequencing (GBS) and four dense genetic maps containing 1449, 996, 720, and 943 SNP markers, respectively, revealed overall a high degree of collinearity for all seven homeologous linkage groups of and . The patterns of distribution of recombination along chromosomes differed between barley and , indicating organizational differences between both genomes.

Race-nonspecific and durable resistance of plant genotypes to major pathogens is highly relevant for yield stability and sustainable crop production but difficult to handle in practice due to its polygenic inheritance by quantitative trait loci (QTL). As far as the underlying genes are concerned, very little is currently known in the most important crop plants such as the cereals. Here, we integrated publicly available data for barley (Hordeum vulgare subsp. vulgare) in order to detect the most important genomic regions for QTL-mediated resistance to a number of fungal pathogens and localize specific functional groups of genes within these regions. This identified 20 meta-QTL, including eight hot spots for resistance to multiple diseases that were distributed over all chromosomes. At least one meta-QTL region for resistance to the powdery mildew fungus Blumeria graminis was found to be co-linear between barley and wheat, suggesting partial evolutionary conservation. Large-scale genetic mapping revealed that functional groups of barley genes involved in secretory processes and cell-wall reinforcement were significantly over-represented within QTL for resistance to powdery mildew. Overall, the results demonstrate added value resulting from large-scale genetic and genomic data integration and may inform genomic-selection procedures for race-nonspecific and durable disease resistance in barley.

Second-generation sequencing now provides the potential for low-cost generation of whole-genome sequences. However, for large-genome organisms with high repetitive DNA content, genome-wide short read sequence assembly is currently impossible, with accurate ordering and localization of genes still relying heavily on integration with physical and genetic maps. To facilitate this process, we have used Agilent microarrays to simultaneously address thousands of gene sequences to individual BAC clones and contiguous sequences that form part of an emerging physical map of the large and currently unsequenced 5.3-Gb barley genome. The approach represents a cost-effective, highly parallel alternative to traditional addressing methods. By coupling the gene-to-BAC address data with gene-based molecular markers, thousands of BACs can be anchored directly to the genetic map, thereby generating a framework for orientating and ordering genes, and providing direct links to phenotypic traits.

The C-value paradox remains incompletely resolved after >40 yr and is exemplified by 2,350-fold variation in genome sizes of flowering plants. The carnivorous Lentibulariaceae genus , displaying a 25-fold range of genome sizes, is a promising subject to study mechanisms and consequences of evolutionary genome size variation. Applying genomic, phylogenetic, and cytogenetic approaches, we uncovered bidirectional genome size evolution within the genus . The Steyerm. genome (86 Mbp) has probably shrunk by retroelement silencing and deletion-biased double-strand break (DSB) repair, from an ancestral size of 400 to 800 Mbp to become one of the smallest among flowering plants. The Stapf genome has expanded by whole-genome duplication (WGD) and retrotransposition to 1550 Mbp. became allotetraploid after the split from the clade ∼29 Ma. A. St.-Hil. (179 Mbp), a close relative of , proved to be a recent (auto)tetraploid. Our analyses suggest a common ancestor of the genus a with an intermediate 1C value (400–800 Mbp) and subsequent rapid genome size evolution in opposite directions. Many abundant repeats of the larger genome are absent in the smaller, casting doubt on their functionality for the organism, while recurrent WGD seems to safeguard against the loss of essential elements in the face of genome shrinkage. We cannot identify any consistent differences in habitat or life strategy that correlate with genome size changes, raising the possibility that these changes may be selectively neutral.

Recent advances in next-generation sequencing techniques allow the detection of a large number of SNPs and their use in a high throughput manner. However, Cleaved Amplified Polymorphic Sequences (CAPSs) still play a significant role as complement to other high throughput methods for SNP genotyping. Therefore, new methods focusing on the acceleration of this type of markers are highly desirable. The combination of the classical CAPS technique and a M13-tailed primer multiplexing assay was used to develop an agarose gel free protocol for the analysis of SNPs via restriction enzyme digestion. PCR products were fluorescence labeled with a universal M13 primer and subsequently digested with the appropriate restriction endonuclease. After mixing differently labeled products, they were detected on a capillary electrophoresis system. This method allows the cost-effective genotyping of several SNPs in a multiplexed manner at an overall low cost in a short period of time. Additionally, this method could be efficiently combined with the simultaneous detection of SSRs at the same electrophoresis run resulting in a procedure well suited for marker-based selection procedures, genotyping of mapping populations and the assay of genetic diversity.

Plant Sequence Capture is used for targeted resequencing of whole exomes (all exons of a genome) of complex genomes e.g. barley and its relatives (Mascher et al., 2013). Sequencing and computing costs are significantly reduced since only the greatly enriched and gene-coding part of the barley genome is targeted, that corresponds to only 1-2% of the entire genome. Thus, applications such as genetic diversity studies and the isolation of single genes (“cloning-by-sequencing”) are greatly facilitated. Here, a protocol is provided describing the construction of shotgun DNA libraries from genomic barley DNA for sequencing on the Illumina HiSeq/MiSeq systems. The shotgun DNA sequencing libraries are hybridized to an oligonucleotide pool (Exome Library) encompassing the whole exome of barley. The Exome Library is provided as a liquid array containing biotinylated probes (Roche/NimbleGen). Subsequently, genomic shotgun DNA fragments hybridized to the Exome Library are affinity-purified using streptavidin coated magnetic beads. The captured library is PCR-amplified and sequenced using high-throughput short read sequencing-by-synthesis.

Barley is a model plant in genomic studies of Triticeae species. However, barley’s large genome size and high repetitive sequence content complicate the whole-genome sequencing. The majority of the barley genome is composed of transposable elements (TEs). In this study, TE repeat junctions (RJs) were used to develop a large-scale molecular marker platform, as a prerequisite to genome assembly. A total of 10.22 Gb of barley nonassembled 454 sequencing data were screened with RJPrimers pipeline. In total, 9,881,561 TE junctions were identified. From detected RJs, 400,538 polymerase chain reaction (PCR)-based RJ markers (RJMs) were designed across the genome, with an average of 39 markers/Mb. The utility of designed markers was tested using a random subset of RJMs. Over 94% of the markers successfully amplified amplicons, among which ∼90% were genome specific. In addition to marker design, identified RJs were utilized to detect 1190,885 TEs across the genome. In gene-poor regions of the genome elements comprised the majority of TEs (∼65%), while in gene-rich regions , , and were the main transposons, each representing an average ∼23% of total TEs. The numerous RJ primer pairs developed in this study will be a valuable resource for barley genomic studies including genomic selection, fine mapping, and genome assembly. In addition, the results of this study show that characterizing RJs provides insight into TE composition of species without a sequenced genome but for which short-read sequence data is available.

Using Roche/454 technology, we sequenced the chloroplast genomes of 12 Triticeae species, including bread wheat, barley and rye, as well as the diploid progenitors and relatives of bread wheat Triticum urartu, Aegilops speltoides and Ae. tauschii. Two wild tetraploid taxa, Ae. cylindrica and Ae. geniculata, were also included. Additionally, we incorporated wild Einkorn wheat Triticum boeoticum and its domesticated form T. monococcum and two Hordeum spontaneum (wild barley) genotypes. Chloroplast genomes were used for overall sequence comparison, phylogenetic analysis and dating of divergence times. We estimate that barley diverged from rye and wheat approximately 8-9 million years ago (MYA). The genome donors of hexaploid wheat diverged between 2.1-2.9 MYA, while rye diverged from Triticum aestivum approximately 3-4 MYA, more recently than previously estimated. Interestingly, the A genome taxa T. boeoticum and T. urartu were estimated to have diverged approximately 570,000 years ago. As these two have a reproductive barrier, the divergence time estimate also provides an upper limit for the time required for the formation of a species boundary between the two. Furthermore, we conclusively show that the chloroplast genome of hexaploid wheat was contributed by the B genome donor and that this unknown species diverged from Ae. speltoides about 980,000 years ago. Additionally, sequence alignments identified a translocation of a chloroplast segment to the nuclear genome which is specific to the rye/wheat lineage. We propose the presented phylogeny and divergence time estimates as a reference framework for future studies on Triticeae.

The rapid development of next-generation sequencing platforms has enabled the use of sequencing for routine genotyping across a range of genetics studies and breeding applications. Genotyping-by-sequencing (GBS), a low-cost, reduced representation sequencing method, is becoming a common approach for whole-genome marker profiling in many species. With quickly developing sequencing technologies, adapting current GBS methodologies to new platforms will leverage these advancements for future studies. To test new semiconductor sequencing platforms for GBS, we genotyped a barley recombinant inbred line (RIL) population. Based on a previous GBS approach, we designed bar code and adapter sets for the Ion Torrent platforms. Four sets of 24-plex libraries were constructed consisting of 94 RILs and the two parents and sequenced on two Ion platforms. In parallel, a 96-plex library of the same RILs was sequenced on the Illumina HiSeq 2000. We applied two different computational pipelines to analyze sequencing data; the reference-independent TASSEL pipeline and a reference-based pipeline using SAMtools. Sequence contigs positioned on the integrated physical and genetic map were used for read mapping and variant calling. We found high agreement in genotype calls between the different platforms and high concordance between genetic and reference-based marker order. There was, however, paucity in the number of SNP that were jointly discovered by the different pipelines indicating a strong effect of alignment and filtering parameters on SNP discovery. We show the utility of the current barley genome assembly as a framework for developing very low-cost genetic maps, facilitating high resolution genetic mapping and negating the need for developing de novo genetic maps for future studies in barley. Through demonstration of GBS on semiconductor sequencing platforms, we conclude that the GBS approach is amenable to a range of platforms and can easily be modified as new sequencing technologies, analysis tools and genomic resources develop.

Numerous CCT domain genes are known to control flowering in plants. They belong to the CONSTANS-like (COL) and PREUDORESPONSE REGULATOR (PRR) gene families, which in addition to a CCT domain possess B-box or response-regulator domains, respectively. Ghd7 is the most recently identified COL gene to have a proven role in the control of flowering time in the Poaceae. However, as it lacks B-box domains, its inclusion within the COL gene family, technically, is incorrect. Here, we show Ghd7 belongs to a larger family of previously uncharacterized Poaceae genes which possess just a single CCT domain, termed here CCT MOTIF FAMILY (CMF) genes. We molecularly describe the CMF (and related COL and PRR) gene families in four sequenced Poaceae species, as well as in the draft genome assembly of barley (Hordeum vulgare). Genetic mapping of the ten barley CMF genes identified, as well as twelve previously unmapped HvCOL and HvPRR genes, finds the majority map to colinear positions relative to their Poaceae orthologues. Combined inter-/intra-species comparative and phylogenetic analysis of CMF, COL and PRR gene families indicates they evolved prior to the monocot/dicot divergence ∼200 mya, with Poaceae CMF evolution described as the interplay between whole genome duplication in the ancestral cereal, and subsequent clade-specific mutation, deletion and duplication events. Given the proven role of CMF genes in the modulation of cereals flowering, the molecular, phylogenetic and comparative analysis of the Poaceae CMF, COL and PRR gene families presented here provides the foundation from which functional investigation can be undertaken.