TRITEX: chromosome-scale sequence assembly of Triticeae genomes with open-source tools

Genome Assembly of the Canadian two-row Malting Barley cultivar AAC Synergy

Barley (Hordeum vulgare L.) is one of the most important global crops. The six-row barley cultivar Morex reference genome has been used by the barley research community worldwide. However, this reference genome can have limitations when used for genomic and genetic diversity analysis studies, gene discovery, and marker development when working in two-row germplasm that is more common to Canadian barley. Here we assembled, for the first time, the genome sequence of a Canadian two-row malting barley, cultivar AAC Synergy. We applied deep Illumina paired-end reads, long mate-pair reads, PacBio sequences, 10X chromium linked read libraries, and chromosome conformation capture sequencing (Hi-C) to generate a contiguous assembly. The genome assembled from super-scaffolds had a size of 4.85 Gb, N50 of 2.32 Mb, and an estimated 93.9% of complete genes from a plant database (BUSCO, benchmarking universal single-copy orthologous genes). After removal of small scaffolds (< 300 Kb), the assembly was arranged into pseudomolecules of 4.14 Gb in size with seven chromosomes plus unanchored scaffolds. The completeness and annotation of the assembly were assessed by comparing it with the updated version of six-row Morex and recently released two-row Golden Promise genome assemblies.

Building pan-genome infrastructures for crop plants and their use in association genetics

Pan-genomic studies aim at representing the entire sequence diversity within a species to provide useful resources for evolutionary studies, functional genomics and breeding of cultivated plants. Cost reductions in high-throughput sequencing and advances in sequence assembly algorithms have made it possible to create multiple reference genomes along with a catalogue of all forms of genetic variations in plant species with large and complex or polyploid genomes. In this review, we summarize the current approaches to building pan-genomes as an in silico representation of plant sequence diversity and outline relevant methods for their effective utilization in linking structural with phenotypic variation. We propose as future research avenues (i) transcriptomic and epigenomic studies across multiple reference genomes and (ii) the development of user-friendly and feature-rich pan-genome browsers.

Improved genome assembly and pan-genome provide key insights into eggplant domestication and breeding

Eggplant (Solanum melongena L.) is an important horticultural crop and one of the most widely grown vegetables from the Solanaceae family. It was domesticated from a wild, prickly progenitor carrying small, round, non-anthocyanic fruits. We obtained a novel, highly contiguous genome assembly of the eggplant '67/3' reference line, by Hi-C retrofitting of a previously released short read- and optical mapping-based assembly. The sizes of the 12 chromosomes and the fraction of anchored genes in the improved assembly were comparable to those of a chromosome-level assembly. We resequenced 23 accessions of S. melongena representative of the worldwide phenotypic, geographic, and genetic diversity of the species, and one each from the closely related species Solanum insanum and Solanum incanum. The eggplant pan-genome contained approximately 51.5 additional megabases and 816 additional genes compared with the reference genome, while the pan-plastome showed little genetic variation. We identified 53 selective sweeps related to fruit color, prickliness, and fruit shape in the nuclear genome, highlighting selection leading to the emergence of present-day S. melongena cultivars from its wild ancestors. Candidate genes underlying the selective sweeps included a MYBL1 repressor and CHALCONE ISOMERASE (for fruit color), homologs of Arabidopsis GLABRA1 and GLABROUS INFLORESCENCE STEMS2 (for prickliness), and orthologs of tomato FW2.2, OVATE, LOCULE NUMBER/WUSCHEL, SUPPRESSOR OF OVATE, and CELL SIZE REGULATOR (for fruit size/shape), further suggesting that selection for the latter trait relied on a common set of orthologous genes in tomato and eggplant.

Identification and spatio-temporal expression analysis of barley genes that encode putative modular xylanolytic enzymes

Arabinoxylans are cell wall polysaccharides whose re-modelling and degradation during plant development are mediated by several classes of xylanolytic enzymes. Here, we present the identification and new annotation of twelve putative (1,4)-β-xylanase and six β-xylosidase genes, and their spatio-temporal expression patterns during vegetative and reproductive growth of barley (Hordeum vulgare cv. Navigator). The encoded xylanase proteins are all predicted to contain a conserved carbohydrate-binding module (CBM) and a catalytic glycoside hydrolase (GH) 10 domain. Additional domains in some xylanases define three discrete phylogenetic clades: one clade contains proteins with an additional N-terminal signal sequence, while another clade contains proteins with multiple CBMs. Homology modelling revealed that all fifteen xylanases likely contain a third domain, a β-sandwich folded from two non-contiguous sequence segments that bracket the catalytic GH domain, which may explain why the full length protein is required for correct folding of the active enzyme. Similarly, predicted xylosidase proteins share a highly conserved domain structure, each with an N-terminal signal peptide, a split GH 3 domain, and a C-terminal fibronectin-like domain. Several genes appear to be ubiquitously expressed during barley growth and development, while four newly annotated xylanase and xylosidase genes are expressed at extremely high levels, which may be of broader interest for industrial applications where cell wall degradation is necessary.

Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly

Until recently, achieving a reference-quality genome sequence for bread wheat was long thought beyond the limits of genome sequencing and assembly technology, primarily due to the large genome size and > 80% repetitive sequence content. The release of the chromosome scale 14.5-Gb IWGSC RefSeq v1.0 genome sequence of bread wheat cv. Chinese Spring (CS) was, therefore, a milestone. Here, we used a direct label and stain (DLS) optical map of the CS genome together with a prior nick, label, repair and stain (NLRS) optical map, and sequence contigs assembled with Pacific Biosciences long reads, to refine the v1.0 assembly. Inconsistencies between the sequence and maps were reconciled and gaps were closed. Gap filling and anchoring of 279 unplaced scaffolds increased the total length of pseudomolecules by 168 Mb (excluding Ns). Positions and orientations were corrected for 233 and 354 scaffolds, respectively, representing 10% of the genome sequence. The accuracy of the remaining 90% of the assembly was validated. As a result of the increased contiguity, the numbers of transposable elements (TEs) and intact TEs have increased in IWGSC RefSeq v2.1 compared with v1.0. In total, 98% of the gene models identified in v1.0 were mapped onto this new assembly through development of a dedicated approach implemented in the MAGAAT pipeline. The numbers of high-confidence genes on pseudomolecules have increased from 105 319 to 105 534. The reconciled assembly enhances the utility of the sequence for genetic mapping, comparative genomics, gene annotation and isolation, and more general studies on the biology of wheat.

Analysis of barley mutants ert-c.1 and ert-d.7 reveals two loci with additive effect on plant architecture

Both mutant ert-c.1 and ert-d.7 carry T2-T3 translocations in the Ert-c gene. Principal coordinate analyses revealed the translocation types and translocation breakpoints. Mutant ert-d.7 is an Ert-c Ert-d double mutant. Mutations in the Ert-c and Ert-d loci are among the most common barley mutations affecting plant architecture. The mutants have various degrees of erect and compact spikes, often accompanied with short and stiff culms. In the current study, complementation tests, linkage mapping, principal coordinate analyses and fine mapping were conducted. We conclude that the original ert-d.7 mutant does not only carry an ert-d mutation but also an ert-c mutation. Combined, mutations in Ert-c and Ert-d cause a pyramid-dense spike phenotype, whereas mutations in only Ert-c or Ert-d give a pyramid and dense phenotype, respectively. Associations between the Ert-c gene and T2-T3 translocations were detected in both mutant ert-c.1 and ert-d.7. Different genetic association patterns indicate different translocation breakpoints in these two mutants. Principal coordinate analysis based on genetic distance and screening of recombinants from all four ends of polymorphic regions was an efficient way to narrow down the region of interest in translocation-involved populations. The Ert-c gene was mapped to the marker interval of 2_0801to1_0224 on 3HL near the centromere. The results illuminate a complex connection between two single genes having additive effects on barley spike architecture and will facilitate the identification of the Ert-c and Ert-d genes.

Improvements in the sequencing and assembly of plant genomes

Advances in DNA sequencing have made it easier to sequence and assemble plant genomes. Here, we extend an earlier study, and compare recent methods for long read sequencing and assembly. Updated Oxford Nanopore Technology software improved assemblies. Using more accurate sequences produced by repeated sequencing of the same molecule (Pacific Biosciences HiFi) resulted in less fragmented assembly of sequencing reads. Using data for increased genome coverage resulted in longer contigs, but reduced total assembly length and improved genome completeness. The original model species, Macadamia jansenii, was also compared with three other Macadamia species, as well as avocado (Persea americana) and jojoba (Simmondsia chinensis). In these angiosperms, increasing sequence data volumes caused a linear increase in contig size, decreased assembly length and further improved already high completeness. Differences in genome size and sequence complexity influenced the success of assembly. Advances in long read sequencing technology continue to improve plant genome sequencing and assembly. However, results were improved by greater genome coverage, with the amount needed to achieve a particular level of assembly being species dependent.

Genome-wide identification of bZIP transcription factor genes related to starch synthesis in barley (Hordeum vulgare L.)

The basic leucine zipper (bZIP) family of genes encode transcription factors that play key roles in plant growth and development. In this study, a total of 92 HvbZIP genes were identified and compared with previous studies using recently released barley genome data. Two novel genes were characterized in this study, and some misannotated and duplicated genes from previous studies have been corrected. Phylogenetic analysis results showed that 92 HvbZIP genes were classified into 10 groups and three unknown groups. The gene structure and motif distribution of the three unknown groups implied that the genes of the three groups may be functionally different. Expression profiling indicated that the HvbZIP genes exhibited different patterns of spatial and temporal expression. Using qRT-PCR, more than 10 HvbZIP genes were identified with expression patterns similar to those of starch synthase genes in barley. Yeast one-hybrid analysis revealed that two of the HvbZIP genes exhibited in vitro binding activity to the promoter of HvAGP-S. The two HvbZIP genes may be candidate genes for further study to explore the mechanism by which they regulate the synthesis of barley starch.

Genomic Prediction of Grain Yield in a Barley MAGIC Population Modeling Genotype per Environment Interaction

Multi-parent Advanced Generation Inter-crosses (MAGIC) lines have mosaic genomes that are generated shuffling the genetic material of the founder parents following pre-defined crossing schemes. In cereal crops, these experimental populations have been extensively used to investigate the genetic bases of several traits and dissect the genetic bases of epistasis. In plants, genomic prediction models are usually fitted using either diverse panels of mostly unrelated accessions or individuals of biparental families and several empirical analyses have been conducted to evaluate the predictive ability of models fitted to these populations using different traits. In this paper, we constructed, genotyped and evaluated a barley MAGIC population of 352 individuals developed with a diverse set of eight founder parents showing contrasting phenotypes for grain yield. We combined phenotypic and genotypic information of this MAGIC population to fit several genomic prediction models which were cross-validated to conduct empirical analyses aimed at examining the predictive ability of these models varying the sizes of training populations. Moreover, several methods to optimize the composition of the training population were also applied to this MAGIC population and cross-validated to estimate the resulting predictive ability. Finally, extensive phenotypic data generated in field trials organized across an ample range of water regimes and climatic conditions in the Mediterranean were used to fit and cross-validate multi-environment genomic prediction models including G×E interaction, using both genomic best linear unbiased prediction and reproducing kernel Hilbert space along with a non-linear Gaussian Kernel. Overall, our empirical analyses showed that genomic prediction models trained with a limited number of MAGIC lines can be used to predict grain yield with values of predictive ability that vary from 0.25 to 0.60 and that beyond QTL mapping and analysis of epistatic effects, MAGIC population might be used to successfully fit genomic prediction models. We concluded that for grain yield, the single-environment genomic prediction models examined in this study are equivalent in terms of predictive ability while, in general, multi-environment models that explicitly split marker effects in main and environmental-specific effects outperform simpler multi-environment models.

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.

Aluminum or Low pH - Which Is the Bigger Enemy of Barley? Transcriptome Analysis of Barley Root Meristem Under Al and Low pH Stress

Aluminum (Al) toxicity is considered to be the most harmful abiotic stress in acidic soils that today comprise more than 50% of the world’s arable lands. Barley belongs to a group of crops that are most sensitive to Al in low pH soils. We present the RNA-seq analysis of root meristems of barley seedlings grown in hydroponics at optimal pH (6.0), low pH (4.0), and low pH with Al (10 μM of bioavailable Al³⁺ ions). Two independent experiments were conducted: with short-term (24 h) and long-term (7 days) Al treatment. In the short-term experiment, more genes were differentially expressed (DEGs) between root meristems grown at pH = 6.0 and pH = 4.0, than between those grown at pH = 4.0 with and without Al treatment. The genes upregulated by low pH were associated mainly with response to oxidative stress, cell wall organization, and iron ion binding. Among genes upregulated by Al, overrepresented were those related to response to stress condition and calcium ion binding. In the long-term experiment, the number of DEGs between hydroponics at pH = 4.0 and 6.0 were lower than in the short-term experiment, which suggests that plants partially adapted to the low pH. Interestingly, 7 days Al treatment caused massive changes in the transcriptome profile. Over 4,000 genes were upregulated and almost 2,000 genes were downregulated by long-term Al stress. These DEGs were related to stress response, cell wall development and metal ion transport. Based on our results we can assume that both, Al³⁺ ions and low pH are harmful to barley plants. Additionally, we phenotyped the root system of barley seedlings grown in the same hydroponic conditions for 7 days at pH = 6.0, pH = 4.0, and pH = 4.0 with Al. The results correspond to transcriptomic data and show that low pH itself is a stress factor that causes a significant reduction of root growth and the addition of aluminum further increases this reduction. It should be noted that in acidic arable lands, plants are exposed simultaneously to both of these stresses. The presented transcriptome analysis may help to find potential targets for breeding barley plants that are more tolerant to such conditions.

Horizontally Acquired nrDNAs Persist in Low Amounts in Host Hordeum Genomes and Evolve Independently of Native nrDNA

Nuclear ribosomal DNA (nrDNA) has displayed extraordinary dynamics during the evolution of plant species. However, the patterns and evolutionary significance of nrDNA array expansion or contraction are still relatively unknown. Moreover, only little is known of the fate of minority nrDNA copies acquired between species via horizontal transfer. The barley genus Hordeum (Poaceae) represents a good model for such a study, as species of section Stenostachys acquired nrDNA via horizontal transfer from at least five different panicoid genera, causing long-term co-existence of native (Hordeum-like) and non-native (panicoid) nrDNAs. Using quantitative PCR, we investigated copy number variation (CNV) of nrDNA in the diploid representatives of the genus Hordeum. We estimated the copy number of the foreign, as well as of the native ITS types (ribotypes), and followed the pattern of their CNV in relation to the genus' phylogeny, species' genomes size and the number of nrDNA loci. For the native ribotype, we encountered an almost 19-fold variation in the mean copy number among the taxa analysed, ranging from 1689 copies (per 2C content) in H. patagonicum subsp. mustersii to 31342 copies in H. murinum subsp. glaucum. The copy numbers did not correlate with any of the genus' phylogeny, the species' genome size or the number of nrDNA loci. The CNV was high within the recognised groups (up to 13.2 × in the American I-genome species) as well as between accessions of the same species (up to 4×). Foreign ribotypes represent only a small fraction of the total number of nrDNA copies. Their copy numbers ranged from single units to tens and rarely hundreds of copies. They amounted, on average, to between 0.1% (Setaria ribotype) and 1.9% (Euclasta ribotype) of total nrDNA. None of the foreign ribotypes showed significant differences with respect to phylogenetic groups recognised within the sect. Stenostachys. Overall, no correlation was found between copy numbers of native and foreign nrDNAs suggesting the sequestration and independent evolution of native and non-native nrDNA arrays. Therefore, foreign nrDNA in Hordeum likely poses a dead-end by-product of horizontal gene transfer events.

Identification and characterization of HAK/KUP/KT potassium transporter gene family in barley and their expression under abiotic stress

BACKGROUND

HAK/KUP/KT (High-affinity K⁺ transporters/K⁺ uptake permeases/K⁺ transporters) is the largest potassium transporter family in plants, and plays pivotal roles in K⁺ uptake and transport, as well as biotic and abiotic stress responses. However, our understanding of the gene family in barley (Hordeum vulgare L.) is quite limited.

RESULTS

In the present study, we identified 27 barley HAK/KUP/KT genes (hereafter called HvHAKs) through a genome-wide analysis. These HvHAKs were unevenly distributed on seven chromosomes, and could be phylogenetically classified into four clusters. All HvHAK protein sequences possessed the conserved motifs and domains. However, the substantial difference existed among HAK members in cis-acting elements and tissue expression patterns. Wheat had the most orthologous genes to barley HAKs, followed by Brachypodium distachyon, rice and maize. In addition, six barley HAK genes were selected to investigate their expression profiling in response to three abiotic stresses by qRT-PCR, and their expression levels were all up-regulated under salt, hyperosmotic and potassium deficiency treatments.

CONCLUSION

Twenty seven HAK genes (HvHAKs) were identified in barley, and they differ in tissue expression patterns and responses to salt stress, drought stress and potassium deficiency.

FINDER: an automated software package to annotate eukaryotic genes from RNA-Seq data and associated protein sequences

BACKGROUND

Gene annotation in eukaryotes is a non-trivial task that requires meticulous analysis of accumulated transcript data. Challenges include transcriptionally active regions of the genome that contain overlapping genes, genes that produce numerous transcripts, transposable elements and numerous diverse sequence repeats. Currently available gene annotation software applications depend on pre-constructed full-length gene sequence assemblies which are not guaranteed to be error-free. The origins of these sequences are often uncertain, making it difficult to identify and rectify errors in them. This hinders the creation of an accurate and holistic representation of the transcriptomic landscape across multiple tissue types and experimental conditions. Therefore, to gauge the extent of diversity in gene structures, a comprehensive analysis of genome-wide expression data is imperative.

RESULTS

We present FINDER, a fully automated computational tool that optimizes the entire process of annotating genes and transcript structures. Unlike current state-of-the-art pipelines, FINDER automates the RNA-Seq pre-processing step by working directly with raw sequence reads and optimizes gene prediction from BRAKER2 by supplementing these reads with associated proteins. The FINDER pipeline (1) reports transcripts and recognizes genes that are expressed under specific conditions, (2) generates all possible alternatively spliced transcripts from expressed RNA-Seq data, (3) analyzes read coverage patterns to modify existing transcript models and create new ones, and (4) scores genes as high- or low-confidence based on the available evidence across multiple datasets. We demonstrate the ability of FINDER to automatically annotate a diverse pool of genomes from eight species.

CONCLUSIONS

FINDER takes a completely automated approach to annotate genes directly from raw expression data. It is capable of processing eukaryotic genomes of all sizes and requires no manual supervision-ideal for bench researchers with limited experience in handling computational tools.

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 discovery of G-quadruplexes in barley

G-quadruplexes (G4s) are four-stranded nucleic acid structures with closely spaced guanine bases forming square planar G-quartets. Aberrant formation of G4 structures has been associated with genomic instability. However, most plant species are lacking comprehensive studies of G4 motifs. In this study, genome-wide identification of G4 motifs in barley was performed, followed by a comparison of genomic distribution and molecular functions to other monocot species, such as wheat, maize, and rice. Similar to the reports on human and some plants like wheat, G4 motifs peaked around the 5′ untranslated region (5′ UTR), the first coding domain sequence, and the first intron start sites on antisense strands. Our comparative analyses in human, Arabidopsis, maize, rice, and sorghum demonstrated that the peak points could be erroneously merged into a single peak when large window sizes are used. We also showed that the G4 distributions around genic regions are relatively similar in the species studied, except in the case of Arabidopsis. G4 containing genes in monocots showed conserved molecular functions for transcription initiation and hydrolase activity. Additionally, we provided examples of imperfect G4 motifs.

Episodes of gene flow and selection during the evolutionary history of domesticated barley

Background

Barley is one of the founder crops of Neolithic agriculture and is among the most-grown cereals today. The only trait that universally differentiates the cultivated and wild subspecies is ‘non-brittleness’ of the rachis (the stem of the inflorescence), which facilitates harvesting of the crop. Other phenotypic differences appear to result from facultative or regional selective pressures. The population structure resulting from these regional events has been interpreted as evidence for multiple domestications or a mosaic ancestry involving genetic interaction between multiple wild or proto-domesticated lineages. However, each of the three mutations that confer non-brittleness originated in the western Fertile Crescent, arguing against multiregional origins for the crop.

Results

We examined exome data for 310 wild, cultivated and hybrid/feral barley accessions and showed that cultivated barley is structured into six genetically-defined groups that display admixture, resulting at least in part from two or more significant passages of gene flow with distinct wild populations. The six groups are descended from a single founding population that emerged in the western Fertile Crescent. Only a few loci were universally targeted by selection, the identity of these suggesting that changes in seedling emergence and pathogen resistance could represent crucial domestication switches. Subsequent selection operated on a regional basis and strongly contributed to differentiation of the genetic groups.

Conclusions

Identification of genetically-defined groups provides clarity to our understanding of the population history of cultivated barley. Inference of population splits and mixtures together with analysis of selection sweeps indicate descent from a single founding population, which emerged in the western Fertile Crescent. This founding population underwent relatively little genetic selection, those changes that did occur affecting traits involved in seedling emergence and pathogen resistance, indicating that these phenotypes should be considered as ‘domestication traits’. During its expansion out of the western Fertile Crescent, the crop underwent regional episodes of gene flow and selection, giving rise to a modern genetic signature that has been interpreted as evidence for multiple domestications, but which we show can be rationalized with a single origin.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07511-7.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

A highly differentiated region of wheat chromosome 7AL encodes a Pm1a immune receptor that recognizes its corresponding AvrPm1a effector from Blumeria graminis

Pm1a, the first powdery mildew resistance gene described in wheat, is part of a complex resistance (R) gene cluster located in a distal region of chromosome 7AL that has suppressed genetic recombination. A nucleotide-binding, leucine-rich repeat (NLR) immune receptor gene was isolated using mutagenesis and R gene enrichment sequencing (MutRenSeq). Stable transformation confirmed Pm1a identity which induced a strong resistance phenotype in transgenic plants upon challenge with avirulent Blumeria graminis (wheat powdery mildew) pathogens. A high-density genetic map of a B. graminis family segregating for Pm1a avirulence combined with pathogen genome resequencing and RNA sequencing (RNAseq) identified AvrPm1a effector gene candidates. In planta expression identified an effector, with an N terminal Y/FxC motif, that induced a strong hypersensitive response when co-expressed with Pm1a in Nicotiana benthamiana. Single chromosome enrichment sequencing (ChromSeq) and assembly of chromosome 7A suggested that suppressed recombination around the Pm1a region was due to a rearrangement involving chromosomes 7A, 7B and 7D. The cloning of Pm1a and its identification in a highly rearranged region of chromosome 7A provides insight into the role of chromosomal rearrangements in the evolution of this complex resistance cluster.

High-resolution mapping of Rym14Hb, a wild relative resistance gene to barley yellow mosaic disease

We mapped the Rym14^Hb resistance locus to barley yellow mosaic disease in a 2Mbp interval. The co-segregating markers will be instrumental for marker-assisted selection in barley breeding. Barley yellow mosaic disease is caused by Barley yellow mosaic virus and Barley mild mosaic virus and leads to severe yield losses in barley (Hordeum vulgare) in Central Europe and East-Asia. Several resistance loci are used in barley breeding. However, cases of resistance-breaking viral strains are known, raising concerns about the durability of those genes. Rym14^Hb is a dominant major resistance gene on chromosome 6HS, originating from barley's secondary genepool wild relative Hordeum bulbosum. As such, the resistance mechanism may represent a case of non-host resistance, which could enhance its durability. A susceptible barley variety and a resistant H. bulbosum introgression line were crossed to produce a large F₂ mapping population (n = 7500), to compensate for a ten-fold reduction in recombination rate compared to intraspecific barley crosses. After high-throughput genotyping, the Rym14^Hb locus was assigned to a 2Mbp telomeric interval on chromosome 6HS. The co-segregating markers developed in this study can be used for marker-assisted introgression of this locus into barley elite germplasm with a minimum of linkage drag.

A Panicum-derived chromosomal segment captured by Hordeum a few million years ago preserves a set of stress-related genes

Intra-specific variability is a cornerstone of evolutionary success of species. Acquiring genetic material from distant sources is an important adaptive mechanism in bacteria, but it can also play a role in eukaryotes. In this paper, we investigate the nature and evolution of a chromosomal segment of panicoid (Poaceae, Panicoideae) origin occurring in the nuclear genomes of species of the barley genus Hordeum (Pooideae). The segment, spanning over 440 kb in the Asian Hordeum bogdanii and 219 kb in the South American Hordeum pubiflorum, resides on a pair of nucleolar organizer region (NOR)-bearing chromosomes. Conserved synteny and micro-collinearity of the segment in both species indicate a common origin of the segment, which was acquired before the split of the respective barley lineages 5-1.7 million years ago. A major part of the foreign DNA consists of several approximately 68 kb long repeated blocks containing five stress-related protein-coding genes and transposable elements (TEs). Whereas outside these repeats, the locus was invaded by multiple TEs from the host genome, the repeated blocks are rather intact and appear to be preserved. The protein-coding genes remained partly functional, as indicated by conserved reading frames, a low amount of non-synonymous mutations, and expression of mRNA. A screen across Hordeum species targeting the panicoid protein-coding genes revealed the presence of the genes in all species of the section Stenostachys. In summary, our study shows that grass genomes can contain large genomic segments obtained from distantly related species. These segments usually remain undetected, but they may play an important role in the evolution and adaptation of species.

An Induced Mutation in HvRECQL4 Increases the Overall Recombination and Restores Fertility in a Barley HvMLH3 Mutant Background

Plant breeding relies on the meiotic recombination or crossing over to generate the new combinations of the alleles along and among the chromosomes. However, crossing over is constrained in the crops such as barley by a combination of the low frequency and biased distribution. In this study, we attempted to identify the genes that limit the recombination by performing a suppressor screen for the restoration of fertility to the semi-fertile barley mutant desynaptic10 (des10), carrying a mutation in the barley ortholog of MutL-Homolog 3 (HvMLH3), a member of the MutL-homolog (MLH) family of DNA mismatch repair genes. des10 mutants exhibit reduced recombination and fewer chiasmata, resulting in the loss of obligate crossovers (COs) leading to chromosome mis-segregation. We identified several candidate suppressor lines and confirmed their restored fertility in an Hvmlh3 background in the subsequent generations. We focus on one of the candidate suppressor lines, SuppLine2099, which showed the most complete restoration of fertility. We characterized this line by using a target-sequence enrichment and sequencing (TENSEQ) capture array representing barley orthologs of 46 meiotic genes. We found that SuppLine2099 contained a C/T change in the anti-CO gene RecQ-like helicase 4 (RECQL4) resulting in the substitution of a non-polar glycine to a polar aspartic acid (G700D) amino acid in the conserved helicase domain. Single nucleotide polymorphism (SNP) genotyping of F₃ populations revealed a significant increase in the recombination frequency in lines with Hvrecql4 in the Hvmlh3 background that was associated with the restoration of fertility. The genotyping also indicated that there was nearly double the recombination levels in homozygous Hvrecql4 lines compared to the wild type (WT). However, we did not observe any significant change in the distribution of CO events. Our results confirm the anti-CO role of RECQL4 in a large genome cereal and establish the possibility of testing the utility of increasing recombination in the context of traditional crop improvement.

Re-examination of the APETALA2/Ethylene-Responsive Factor Gene Family in Barley (Hordeum vulgare L.) Indicates a Role in the Regulation of Starch Synthesis

The APETALA2/Ethylene-Responsive factor (AP2/ERF) gene family is a large plant-specific transcription factor family, which plays important roles in regulating plant growth and development. A role in starch synthesis is among the multiple functions of this family of transcription factors. Barley (Hordeum vulgare L.) is one of the most important cereals for starch production. However, there are limited data on the contribution of AP2 transcription factors in barley. In this study, we used the recently published barley genome database (Morex) to identify 185 genes of the HvAP2/ERF family. Compared with previous work, we identified 64 new genes in the HvAP2/ERF gene family and corrected some previously misannotated and duplicated genes. After phylogenetic analysis, HvAP2/ERF genes were classified into four subfamilies and 18 subgroups. Expression profiling showed different patterns of spatial and temporal expression for HvAP2/ERF genes. Most of the 12 HvAP2/ERF genes analyzed using quantitative reverse transcription-polymerase chain reaction had similar expression patterns when compared with those of starch synthase genes in barley, except for HvAP2-18 and HvERF-73. HvAP2-18 is homologous to OsRSR1, which negatively regulates the synthesis of rice starch. Luciferase reporter gene, and yeast one-hybrid assays showed that HvAP2-18 bound the promoter of AGP-S and SBE1 in vitro. Thus, HvAP2-18 might be an interesting candidate gene to further explore the mechanisms involved in the regulation of starch synthesis in barley.

Identification, Characterization, and Expression Profile Analysis of the mTERF Gene Family and Its Role in the Response to Abiotic Stress in Barley (Hordeum vulgare L.)

Plant mitochondrial transcription termination factor (mTERF) family regulates organellar gene expression (OGE) and is functionally characterized in diverse species. However, limited data are available about its functions in the agriculturally important cereal barley (Hordeum vulgare L.). In this study, we identified 60 mTERFs in the barley genome (HvmTERFs) through a comprehensive search against the most updated barley reference genome, Morex V2. Then, phylogenetic analysis categorized these genes into nine subfamilies, with approximately half of the HvmTERFs belonging to subfamily IX. Members within the same subfamily generally possessed conserved motif composition and exon-intron structure. Both segmental and tandem duplication contributed to the expansion of HvmTERFs, and the duplicated gene pairs were subjected to strong purifying selection. Expression analysis suggested that many HvmTERFs may play important roles in barley development (e.g., seedlings, leaves, and developing inflorescences) and abiotic stresses (e.g., cold, salt, and metal ion), and HvmTERF21 and HvmTERF23 were significant induced by various abiotic stresses and/or phytohormone treatment. Finally, the nucleotide diversity was decreased by only 4.5% for HvmTERFs during the process of barley domestication. Collectively, this is the first report to characterize HvmTERFs, which will not only provide important insights into further evolutionary studies but also contribute to a better understanding of the potential functions of HvmTERFs and ultimately will be useful in future gene functional studies.

Chromosome-scale genome assembly provides insights into rye biology, evolution and agronomic potential

Rye (Secale cereale L.) is an exceptionally climate-resilient cereal crop, used extensively to produce improved wheat varieties via introgressive hybridization and possessing the entire repertoire of genes necessary to enable hybrid breeding. Rye is allogamous and only recently domesticated, thus giving cultivated ryes access to a diverse and exploitable wild gene pool. To further enhance the agronomic potential of rye, we produced a chromosome-scale annotated assembly of the 7.9-gigabase rye genome and extensively validated its quality by using a suite of molecular genetic resources. We demonstrate applications of this resource with a broad range of investigations. We present findings on cultivated rye's incomplete genetic isolation from wild relatives, mechanisms of genome structural evolution, pathogen resistance, low-temperature tolerance, fertility control systems for hybrid breeding and the yield benefits of rye-wheat introgressions.

Recognition and defence of plant-infecting fungal pathogens

Attempted infections of plants with fungi result in diverse outcomes ranging from symptom-less resistance to severe disease and even death of infected plants. The deleterious effect on crop yield have led to intense focus on the cellular and molecular mechanisms that explain the difference between resistance and susceptibility. This research has uncovered plant resistance or susceptibility genes that explain either dominant or recessive inheritance of plant resistance with many of them coding for receptors that recognize pathogen invasion. Approaches based on cell biology and phytochemistry have contributed to identifying factors that halt an invading fungal pathogen from further invasion into or between plant cells. Plant chemical defence compounds, antifungal proteins and structural reinforcement of cell walls appear to slow down fungal growth or even prevent fungal penetration in resistant plants. Additionally, the hypersensitive response, in which a few cells undergo a strong local immune reaction, including programmed cell death at the site of infection, stops in particular biotrophic fungi from spreading into surrounding tissue. In this review, we give a general overview of plant recognition and defence of fungal parasites tracing back to the early 20th century with a special focus on Triticeae and on the progress that was made in the last 30 years.

Genomic Dissection of Peduncle Morphology in Barley through Nested Association Mapping

Straw biomass and stability are crucial for stable yields. Moreover, straw harbors the potential to serve as a valuable raw material for bio-economic processes. The peduncle is the top part of the last shoot internode and carries the spike. This study investigates the genetic control of barley peduncle morphology. Therefore, 1411 BC1S3 lines of the nested association mapping (NAM) population “Halle Exotic Barley 25” (HEB-25), generated by crossing the spring barley elite cultivar Barke with an assortment of 25 exotic barley accessions, were used. Applying 50k Illumina Infinium iSelect SNP genotyping yielded new insights and a better understanding of the quantitative trait loci (QTL) involved in controlling the peduncle diameter traits, we found the total thickness of peduncle tissues and the area of the peduncle cross-section. We identified three major QTL regions on chromosomes 2H and 3H mainly impacting the traits. Remarkably, the exotic allele at the QTL on chromosome 3H improved all three traits investigated in this work. Introgressing this QTL in elite cultivars might facilitate to adjust peduncle morphology for improved plant stability or enlarged straw biomass production independent of flowering time and without detrimental effects on grain yield.

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.

A genetic analysis of the resistance in barley to Soil-borne wheat mosaic virus

Soil-borne wheat mosaic virus (SBWMV), a ubiquitous pathogen commonly encountered in temperate regions of the Northern hemisphere, can damage a number of economically important cereal crops, notably wheat and barley. Given that the plasmodiophorid cercozoan Polymyxa graminis, which acts as the vector of SBWMV, can survive in the soil for many decades, the only feasible control measure is the deployment of resistant cultivars. Here, a quantitative trait locus (QTL) approach was taken to characterize the genetic basis of the SBWMV resistance exhibited by the barley cultivar Haruna Nijo. The analysis revealed that between 33% and 41% of the variation for the measure chosen to represent resistance was under the control of a gene(s) mapping to a region at the distal end of the short arm of chromosome 2H. In contrast to most of the genes known to encode resistance to soil-borne mosaic viruses, the allele specifying resistance was dominant over those present in a susceptible genotype.

The barley pan-genome reveals the hidden legacy of mutation breeding

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Identification of polycomb repressive complex 1 and 2 core components in hexaploid bread wheat

BACKGROUND

Polycomb repressive complexes 1 and 2 play important roles in epigenetic gene regulation by posttranslationally modifying specific histone residues. Polycomb repressive complex 2 is responsible for the trimethylation of lysine 27 on histone H3; Polycomb repressive complex 1 catalyzes the monoubiquitination of histone H2A at lysine 119. Both complexes have been thoroughly studied in Arabidopsis, but the evolution of polycomb group gene families in monocots, particularly those with complex allopolyploid origins, is unknown.

RESULTS

Here, we present the in silico identification of the Polycomb repressive complex 1 and 2 (PRC2, PRC1) subunits in allohexaploid bread wheat, the reconstruction of their evolutionary history and a transcriptional analysis over a series of 33 developmental stages. We identified four main subunits of PRC2 [E(z), Su(z), FIE and MSI] and three main subunits of PRC1 (Pc, Psc and Sce) and determined their chromosomal locations. We found that most of the genes coding for subunit proteins are present as paralogs in bread wheat. Using bread wheat RNA-seq data from different tissues and developmental stages throughout plant ontogenesis revealed variable transcriptional activity for individual paralogs. Phylogenetic analysis showed a high level of protein conservation among temperate cereals.

CONCLUSIONS

The identification and chromosomal location of the Polycomb repressive complex 1 and 2 core components in bread wheat may enable a deeper understanding of developmental processes, including vernalization, in commonly grown winter wheat.

Fonio millet genome unlocks African orphan crop diversity for agriculture in a changing climate

Sustainable food production in the context of climate change necessitates diversification of agriculture and a more efficient utilization of plant genetic resources. Fonio millet (Digitaria exilis) is an orphan African cereal crop with a great potential for dryland agriculture. Here, we establish high-quality genomic resources to facilitate fonio improvement through molecular breeding. These include a chromosome-scale reference assembly and deep re-sequencing of 183 cultivated and wild Digitaria accessions, enabling insights into genetic diversity, population structure, and domestication. Fonio diversity is shaped by climatic, geographic, and ethnolinguistic factors. Two genes associated with seed size and shattering showed signatures of selection. Most known domestication genes from other cereal models however have not experienced strong selection in fonio, providing direct targets to rapidly improve this crop for agriculture in hot and dry environments.

Fonio millet genome unlocks African orphan crop diversity for agriculture in a changing climate

Sustainable food production in the context of climate change necessitates diversification of agriculture and a more efficient utilization of plant genetic resources. Fonio millet (Digitaria exilis) is an orphan African cereal crop with a great potential for dryland agriculture. Here, we establish high-quality genomic resources to facilitate fonio improvement through molecular breeding. These include a chromosome-scale reference assembly and deep re-sequencing of 183 cultivated and wild Digitaria accessions, enabling insights into genetic diversity, population structure, and domestication. Fonio diversity is shaped by climatic, geographic, and ethnolinguistic factors. Two genes associated with seed size and shattering showed signatures of selection. Most known domestication genes from other cereal models however have not experienced strong selection in fonio, providing direct targets to rapidly improve this crop for agriculture in hot and dry environments.

History and future perspectives of barley genomics

Barley (Hordeum vulgare), one of the most widely cultivated cereal crops, possesses a large genome of 5.1 Gbp. Through various international collaborations, the genome has recently been sequenced and assembled at the chromosome-scale by exploiting available genetic and genomic resources. Many wild and cultivated barley accessions have been collected and preserved around the world. These accessions are crucial to obtain diverse natural and induced barley variants. The barley bioresource project aims to investigate the diversity of this crop based on purified seed and DNA samples of a large number of collected accessions. The long-term goal of this project is to analyse the genome sequences of major barley accessions worldwide. In view of technical limitations, a strategy has been employed to establish the exome structure of a selected number of accessions and to perform high-quality chromosome-scale assembly of the genomes of several major representative accessions. For the future project, an efficient annotation pipeline is essential for establishing the function of genomes and genes as well as for using this information for sequence-based digital barley breeding. In this article, the author reviews the existing barley resources along with their applications and discuss possible future directions of research in barley genomics.

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.

A Genome Assembly of the Barley 'Transformation Reference' Cultivar Golden Promise

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.